Changeset 1648 for trunk/LMDZ.TITAN/libf/phytitan/sugas_corrk.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/sugas_corrk.F90 (modified) (6 diffs)

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

@@ -24 +24 @@
       use radcommon_h, only : tgasref,tgasmin,tgasmax
       use radcommon_h, only : gasv,gasi,FZEROI,FZEROV,gweight
-      use radcommon_h, only : wrefvar,WNOI,WNOV
+      use radcommon_h, only : WNOI,WNOV
       use datafile_mod, only: datadir
       use comcstfi_mod, only: mugaz
@@ -44 +44 @@
       ! ALLOCATABLE ARRAYS -- AS 12/2011
       REAL*8, DIMENSION(:,:,:,:,:), ALLOCATABLE,SAVE :: gasi8, gasv8    !read by master
-      character*20,allocatable,DIMENSION(:),SAVE :: gastype ! for check with gnom, read by master
 
       real*8 x, xi(4), yi(4), ans, p
@@ -79 +78 @@
       read(111,*) ngas
 
-      if(ngas.ne.ngasmx)then
-         print*,'Number of gases in radiative transfer data (',ngas,') does not', &
-                'match that in gases.def (',ngasmx,'), exiting.'
-         call abort
-      endif 
-
       if(ngas.gt.5 .or. ngas.lt.1)then
          print*,ngas,' species in database [',               &
@@ -90 +83 @@
                 '], radiative code cannot handle this.'
          call abort
-      endif 
-
-      ! dynamically allocate gastype and read from Q.dat
-      IF ( .NOT. ALLOCATED( gastype ) ) ALLOCATE( gastype( ngas ) )
-
-      do igas=1,ngas
-         read(111,*) gastype(igas)
-         print*,'Gas ',igas,' is ',gastype(igas)
-      enddo
-
-      ! get array size, load the coefficients
-      open(111,file=TRIM(file_path),form='formatted')
-      read(111,*) L_REFVAR
-      IF( .NOT. ALLOCATED( wrefvar ) ) ALLOCATE( WREFVAR(L_REFVAR) )
-      read(111,*) wrefvar
-      close(111)
-
-      ! Check that gastype and gnom match
-      do igas=1,ngas
-         print*,'Gas ',igas,' is ',trim(gnom(igas))
-         if (trim(gnom(igas)).ne.trim(gastype(igas))) then
-            print*,'Name of a gas in radiative transfer data (',trim(gastype(igas)),') does not ', &
-                 'match that in gases.def (',trim(gnom(igas)),'), exiting. You should compare ', &
-                 'gases.def with Q.dat in your radiative transfer directory.' 
-            call abort
-         endif
-      enddo
-      print*,'Confirmed gas match in radiative transfer and gases.def!'
-
-      ! display the values
-      print*,'Variable gas volume mixing ratios:'
-      do n=1,L_REFVAR
-         !print*,n,'.',wrefvar(n),' kg/kg' ! pay attention!
-         print*,n,'.',wrefvar(n),' mol/mol'
-      end do
-      print*,''
+      endif
+      
+      L_REFVAR = 1 ! JVO 2017 : set to 1 to keep the code running until the new variable species treatment
 
 !=======================================================================
@@ -628 +588 @@
             enddo
 
-         endif
-
+         elseif (igas .eq. igas_CH4) then
+
+            ! first do self-induced absorption
+            dummy = -9999
+            call interpolateCH4CH4(200.D+0,200.D+0,7500.D+0,testcont,.true.,dummy)
+            ! then cross-interactions with other gases
+            do jgas=1,ngasmx
+               if (jgas .eq. igas_N2) then
+                  dummy = -9999
+                  call interpolateN2CH4(200.D+0,250.0D+0,100000.D+0,5000.D+0,testcont,.true.,dummy)
+               endif
+            enddo
+
+         endif  
+         
       enddo
       endif
@@ -645 +618 @@
       IF( ALLOCATED( gasv8 ) ) DEALLOCATE( gasv8 )
       IF( ALLOCATED( pgasref ) ) DEALLOCATE( pgasref )
-      IF( ALLOCATED( gastype ) ) DEALLOCATE( gastype )
 !$OMP END MASTER
 !$OMP BARRIER