Changeset 1648 for trunk/LMDZ.TITAN/libf/phytitan/su_gases.F90

Timestamp:

Jan 19, 2017, 2:46:53 PM (8 years ago)

Author:

jvatant

Message:

Modifications to custom radiative transfer to Titan
+ Enables an altitude dependant gfrac for CIA computations

-> many radical changes in su_gases and co ..
-> read vertical CH4 profile with call_profilgases
-> Now you need a 'profile.def' that I will add in the deftank

+ Added interpolate CIA routines for CH4
+ Added temporary mean aerosol profile opacity routine (disr_haze)

File:

• trunk/LMDZ.TITAN/libf/phytitan/su_gases.F90 (modified) (3 diffs)

trunk/LMDZ.TITAN/libf/phytitan/su_gases.F90

@@ -1 @@
-subroutine su_gases
+subroutine su_gases(nlayer,tracer)
 
   use gases_h
@@ -5 +5 @@
   implicit none
 
-  integer igas, ierr, count
-
+  integer,intent(in) :: nlayer
+  logical,intent(in) :: tracer
+  
+  integer igas, ierr
+  
   !==================================================================
   !     
@@ -17 +20 @@
   !     R. Wordsworth (2011)
   !     Allocatable arrays by A. Spiga (2011)
-  !     
+  !     Titan's version : J. Vatant d'Ollone (2017)
+  !              * force igas_X labels and def gnom for N2,CH4,H2
+  !              * get rid of variable species ( enrichment dimension will be set in corrk routines )
+  !
   !==================================================================
 
 !$OMP MASTER
-  ! load gas names from file 'gases.def'
+
+
+  ngasmx   = 3 ! N2, CH4, H2
+
+
+  ! load reference level and reference molar fractions from file 'gases.def'
   open(90,file='gases.def',status='old',form='formatted',iostat=ierr)
+  
   if (ierr.eq.0) then
      write(*,*) "sugases.F90: reading file gases.def"
-     read(90,*)
-     read(90,*,iostat=ierr) ngasmx
+     read(90,*) ! header
+
+     ! We allocate gfrac and we set gas molar fractions for the reference level only.
+     ! This will be useful for the cpp_mu and rayleigh routines
+     ! Other gas molar fractions are now set in routine callprofilgases
+     
+     write(*,*) 'sugases.F90: allocating and reading gas molar fractions from reference level in gases.def...'
+     
+     if(.not.allocated(gfrac)) allocate(gfrac(ngasmx,nlayer))
+     
+     read(90,*,iostat=ierr) nivref     
      if (ierr.ne.0) then
-        write(*,*) "sugases.F90: error reading number of gases"
+        write(*,*) "sugases.F90: error reading reference level"
         write(*,*) "   (first line of gases.def) "
         call abort
      endif
-
-     print*,ngasmx, " gases found in gases.def. Allocating names and molar fractions..."
-
-     if (.not.allocated(gnom)) allocate(gnom(ngasmx))
+     
+     print*, "layer", nivref, "is reference level found in gases.def ..."
+     
      do igas=1,ngasmx
-        read(90,*,iostat=ierr) gnom(igas)
+        read(90,*,iostat=ierr) gfrac(igas,nivref)
         if (ierr.ne.0) then
-           write(*,*) 'sugases.F90: error reading gas names in gases.def...'
+           write(*,*) 'sugases.F90: error reading reference gas molar fractions in gases.def... aborting'
            call abort
         endif
      enddo                  !of do igas=1,ngasmx
 
-     vgas=0
-     if(.not.allocated(gfrac)) allocate(gfrac(ngasmx))
-     do igas=1,ngasmx
-        read(90,*,iostat=ierr) gfrac(igas)
-        if (ierr.ne.0) then
-           write(*,*) 'sugases.F90: error reading gas molar fractions in gases.def...'
-           call abort
-        endif
 
-        ! find variable gas (if any)
-        if(gfrac(igas).eq.-1.0)then 
-           if(vgas.eq.0)then
-              vgas=igas
-           else
-              print*,'You seem to be choosing two variable gases'
-              print*,'Check that gases.def is correct'
-              call abort
-           endif
-        endif
-
-     enddo                  !of do igas=1,ngasmx
+     ! We force gnom = (N2, CH4, H2) and igas_X for Titan
+     
+     igas_N2  = 1
+     igas_CH4 = 2
+     igas_H2  = 3
 
 
-     ! assign the 'igas_X' labels
-     count=0
-     do igas=1,ngasmx
-        if (trim(gnom(igas)).eq."H2_" .or. trim(gnom(igas)).eq."H2") then
-           igas_H2=igas
-           count=count+1
-        elseif (trim(gnom(igas)).eq."N2_" .or. trim(gnom(igas)).eq."N2") then
-           igas_N2=igas
-           count=count+1
-        elseif (trim(gnom(igas)).eq."CH4") then
-           igas_CH4=igas
-           count=count+1
-        elseif (trim(gnom(igas)).eq."C2H6") then
-           igas_C2H6=igas
-           count=count+1
-        elseif (trim(gnom(igas)).eq."C2H2") then
-           igas_C2H2=igas
-           count=count+1
-        endif
-     enddo
-
-     if(count.ne.ngasmx)then
-        print*,'Mismatch between ngas and number of recognised gases in sugas_corrk.F90.'
-        print*,'Either we haven`t managed to assign all the gases, or there are duplicates.'
-        print*,'Please try again.'
-     endif
-
+     if (.not.allocated(gnom)) allocate(gnom(ngasmx))
+     gnom(igas_N2)  = "N2_"
+     gnom(igas_CH4) = "CH4"
+     gnom(igas_H2)  = "H2_"
+     
+     
   else
      write(*,*) 'Cannot find required file "gases.def"'
      call abort
   endif
+  
   close(90)
 !$OMP END MASTER