[강릉원주대 대기환경과학과] 2014년 2학기 전선 구름물리학 과제물

 정보

• 업무명    : 2014년 2학기 전선 구름물리학 과제물

• 작성자    : 이상호

• 작성일    : 2019-09-13

• 설   명     :

• 수정이력 :

  • 2020-03-27 : 서식 수정

 

 내용

[특징]

• 구름물리학 과제물에 대한 이해를 돕기위해 작성

 

[기능]

• 쾰러 곡선 (Kohler Curve)를 이용한 가시화

 

[활용 자료]

• 인터넷 검색

 

[과제물 작성 요령]

• 과제물 작성 계획 세우기

  • 과제물 작성을 위하여 계획 및 준비, 글쓰기, 제출과 같이 3 단계로 구성

  • 1 단계는 우선 계획 및 준비로서 공고 확인, 제출 기간 확인, 과제물 문제를 파악

• 2단계의 경우 글쓰기로서 자료 검색 수집 및 개요 작성 그리고 본문 쓰기 및 수정

• 끝으로 제출 마감일뿐만 아니라 파일 용량 검사 및 온라인시스템 제출

• 과제물 문제 파악 요령

  • 과제물 문제에서 교수님의 지시사항을 정확히 간파하여 작성

  • 특히 제출 파일의 양식/용량 및 스타일 형식 그리고 작성 분량에 유의

• 나만의 글쓰기 전략

  • 앞서 작성 계획에서와 같이 단계별로 수행이 필연

  • 추가로 각 과제물마다 한글 또는 워드 양식에서 각자만의 스타일을 바탕으로 글의 내용 작성

  • 또한 각종 인터넷 검색 및 도서를 참조하여 내용을 제시할 뿐만 아니라 그에 따른 판단의 근거에 대한 명확한 인용 및 출처가 정확하게 표시

• 학술 검색을 위한 유용한 방법

  • 인터넷은 1950대에 컴퓨터의 개발과 더불어 너무 많은 정보의 바다에서 살고 있기 때문에 자신에게 필요한 정보 발굴 능력 요구 

  • 인터넷 검색이 익숙하지 않을 경우 네이버 검색을 추천하고 익숙하면 구글 검색이 유용

  • 그러나 일반적으로 네이버 또는 구글로 검색할 경우 다수의 블로그에서 중복되는 정보를 확인할 수 있고 이러한 정보들은 대부분 복사 붙여넣기로 인해 남용되고 있는 실정

  • 따라서 이 보다는 전 세계 검색엔진 2위를 차지한 Yahoo! JAPAN (https://www.yahoo.co.jp)에서 검색하여 구글 번역할 경우 원하는 정보를 쉽게 얻음

  • 연구 보고서 또는 SCI/SCIE 논문을 보기 위해서 대학교에서 저널을 구독하는 것이 일반적이나 모든 저널을 구독하기에는 상당한 경제적 부담이 되기 때문에 SCI-HUB (https://sci-hub.tw/)을 통해 보는 것을 추천

 

[사용법]

• 과제물 서식을 바탕으로 작성

 

[사용 OS]

• Windows 10

 

[사용 언어]

• Fortran v90

• 한글 v2018

 

 과제물

• NaCl 농도에 따른 가시화

그림. NaCl 농도에 따른 가시화.

 

• NaCl과 (NH4)SO4에 따른 가시화

그림. NaCl과 (NH4)SO4에 따른 가시화.

 

• NaCl과 (NH4)SO4 농도에 따른 가시화

그림. NaCl과 (NH4)SO4 농도에 따른 가시화.

 

[전체]

      real,parameter :: C=7.5*(10.**-2.) , Rv=482. , Pl=1000. , Mna=22.99 , Mcl=35.45 , Mh=1.01 , Mo=16. , Mn=14. , Ms=32.
      real :: T,Mv,i2,r,pi,max,j,a,b,S,Rc,Sc,dd
      real,dimension(5) :: bb 
      real,dimension(3) :: cc
      real,dimension(2) :: i,Mss
    
      open(1, file='1.txt')
      open(2, file='2.txt')
      open(3, file='3.txt')
      open(4, file='4.txt')

      i(1) = 2.
      i(2) = 3.
      bb(1) = 10.**-17.
      bb(2) = 10.**-16.
      bb(3) = 10.**-15.
      bb(4) = 10.**-18.
      bb(5) = 10.**-19.
      pi = atan(1.0d0)*4.0d0
      Mss(1) = Mna+Mcl
      Mss(2) = ((Mn+(Mh*4.))*2.)+(Ms+((Mo)*4.))
      Mv = Mh + Mh + Mo 
      T = 273.+5.
    

      !===========================================
      !  Figure 1
      !===========================================
      do k = 1, 3
         max = 0.1
         do j = 0.01, 10., 0.001
            r = real(j)
            a = (2.*C)/(Pl*Rv*T)*(10.**6.)
            b = (3.*i(1)*bb(k)*Mv)/(4.*pi*Pl*Mss(1))*(10.**15.)
            S = 1.+(a/r)-(b/(r**3.))
            Rc = sqrt((3.*b)/a)
            Sc = 1.+sqrt((4.*(a**3.))/(27.*b))
            if (max.lt.S) max=S
            if (S.eq.max) cc(k)=r
            ! write(1,*) r,S,Rc,max
         enddo 
         ! write(*,*) Rc,cc(k),Sc,max
      enddo

      
      !===========================================
      !  Figure 2
      !===========================================
      do k = 1, 2
         max = 0.1
         do j = 0.01, 10., 0.001
            r = real(j)
            a = (2.*C)/(Pl*Rv*T)*(10.**6.)
            b = (3.*i(k)*bb(2)*Mv)/(4.*pi*Pl*Mss(k))*(10.**15.)
            S = 1.+(a/r)-(b/(r**3.))
            Rc = sqrt((3.*b)/a)
            Sc = 1.+sqrt((4.*(a**3.))/(27.*b))
            if (max.lt.S) max=S
            if (S.eq.max) cc(k)=r
            ! write(2,*) i(k),r,S
         enddo 
         ! write(*,*) Rc,cc(k),Sc,max
      enddo

      !===========================================
      !  Figure 3a
      !===========================================
      max = 0.1
      do j = 0.01, 10., 0.001
         r = real(j)
         a = (2.*C)/(Pl*Rv*T)*(10.**6.)
         b = (3.*i(1)*bb(4)*Mv)/(4.*pi*Pl*Mss(1))*(10.**18.)
         S = 1.+(a/r)-(b/(r**3.))
         Rc = sqrt((3.*b)/a)
         Sc = 1.+sqrt((4.*(a**3.))/(27.*b))
         if (max.lt.S) max=S
         if (S.eq.max) dd=r
         ! write(3,*) r,S
      enddo 
      ! write(*,*) dd, Rc, max, Sc

      !===========================================
      !  Figure 3b
      !===========================================
      max = 0.1
      do j = 0.01,10.,0.001
         r = real(j)
         a = (2.*C)/(Pl*Rv*T)*(10.**6.)
         b = (3.*i(2)*bb(5)*Mv)/(4.*pi*Pl*Mss(2))*(10.**18.)
         S = 1.+(a/r)-(b/(r**3.))
         Rc = sqrt((3.*b)/a)
         Sc = 1.+sqrt((4.*(a**3.))/(27.*b))
         ! write(4,*) r,S
         if (max.lt.S) max=S
         if (S.eq.max) dd=r
      enddo 
      ! write(*,*) max,Sc,dd,Rc

      end

 

[GihHub GIST]

	real,parameter :: C=7.5*(10.**-2.) , Rv=482. , Pl=1000. , Mna=22.99 , Mcl=35.45 , Mh=1.01 , Mo=16. , Mn=14. , Ms=32.
	real :: T,Mv,i2,r,pi,max,j,a,b,S,Rc,Sc,dd
	real,dimension(5) :: bb
	real,dimension(3) :: cc
	real,dimension(2) :: i,Mss
	open(1, file='1.txt')
	open(2, file='2.txt')
	open(3, file='3.txt')
	open(4, file='4.txt')
	i(1) = 2.
	i(2) = 3.
	bb(1) = 10.**-17.
	bb(2) = 10.**-16.
	bb(3) = 10.**-15.
	bb(4) = 10.**-18.
	bb(5) = 10.**-19.
	pi = atan(1.0d0)*4.0d0
	Mss(1) = Mna+Mcl
	Mss(2) = ((Mn+(Mh*4.))*2.)+(Ms+((Mo)*4.))
	Mv = Mh + Mh + Mo
	T = 273.+5.
	!===========================================
	! Figure 1
	!===========================================
	do k = 1, 3
	max = 0.1
	do j = 0.01, 10., 0.001
	r = real(j)
	a = (2.*C)/(Pl*Rv*T)*(10.**6.)
	b = (3.*i(1)*bb(k)*Mv)/(4.*pi*Pl*Mss(1))*(10.**15.)
	S = 1.+(a/r)-(b/(r**3.))
	Rc = sqrt((3.*b)/a)
	Sc = 1.+sqrt((4.*(a**3.))/(27.*b))
	if (max.lt.S) max=S
	if (S.eq.max) cc(k)=r
	! write(1,*) r,S,Rc,max
	enddo
	! write(*,*) Rc,cc(k),Sc,max
	enddo
	!===========================================
	! Figure 2
	!===========================================
	do k = 1, 2
	max = 0.1
	do j = 0.01, 10., 0.001
	r = real(j)
	a = (2.*C)/(Pl*Rv*T)*(10.**6.)
	b = (3.*i(k)*bb(2)*Mv)/(4.*pi*Pl*Mss(k))*(10.**15.)
	S = 1.+(a/r)-(b/(r**3.))
	Rc = sqrt((3.*b)/a)
	Sc = 1.+sqrt((4.*(a**3.))/(27.*b))
	if (max.lt.S) max=S
	if (S.eq.max) cc(k)=r
	! write(2,*) i(k),r,S
	enddo
	! write(*,*) Rc,cc(k),Sc,max
	enddo
	!===========================================
	! Figure 3a
	!===========================================
	max = 0.1
	do j = 0.01, 10., 0.001
	r = real(j)
	a = (2.*C)/(Pl*Rv*T)*(10.**6.)
	b = (3.*i(1)*bb(4)*Mv)/(4.*pi*Pl*Mss(1))*(10.**18.)
	S = 1.+(a/r)-(b/(r**3.))
	Rc = sqrt((3.*b)/a)
	Sc = 1.+sqrt((4.*(a**3.))/(27.*b))
	if (max.lt.S) max=S
	if (S.eq.max) dd=r
	! write(3,*) r,S
	enddo
	! write(*,*) dd, Rc, max, Sc
	!===========================================
	! Figure 3b
	!===========================================
	max = 0.1
	do j = 0.01,10.,0.001
	r = real(j)
	a = (2.*C)/(Pl*Rv*T)*(10.**6.)
	b = (3.*i(2)*bb(5)*Mv)/(4.*pi*Pl*Mss(2))*(10.**18.)
	S = 1.+(a/r)-(b/(r**3.))
	Rc = sqrt((3.*b)/a)
	Sc = 1.+sqrt((4.*(a**3.))/(27.*b))
	! write(4,*) r,S
	if (max.lt.S) max=S
	if (S.eq.max) dd=r
	enddo
	! write(*,*) max,Sc,dd,Rc
	end

view raw Visualization_Using_Kohler_Curve.f90 hosted with ❤ by GitHub

 

 관련 자료

퀼러 곡선를 이용한 자료 처리 및 가시화

[강릉원주대 대기환경과학과] 2014년 2학기 전선 구름물리학 과제물 from Lee Sang-Ho

 

 참고 문헌

[논문]

• 없음

[보고서]

• 없음

[URL]

• 없음

 

 문의사항

[기상학/프로그래밍 언어]

• sangho.lee.1990@gmail.com

[해양학/천문학/빅데이터]

• saimang0804@gmail.com