Changeset 1789 for trunk/LMDZ.TITAN/libf/phytitan/phyetat0_mod.F90

Timestamp:

Sep 28, 2017, 3:43:23 PM (7 years ago)

Author:

jvatant

Message:

Added the surface methane tank and put it in start files
--JVO

File:

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

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

@@ -8 +8 @@
                      ngrid,nlayer,fichnom,tab0,Lmodif,nsoil,nq, &
                      day_ini,time,tsurf,tsoil, &
-                     emis,q2,qsurf)
+                     emis,q2,qsurf,tankCH4)
 
 ! to use  'getin_p'
@@ -44 +44 @@
   real,intent(out) :: q2(ngrid,nlayer+1) ! 
   real,intent(out) :: qsurf(ngrid,nq) ! tracers on surface
+  real,intent(out) :: tankCH4(ngrid)  ! depth of CH4 tank
 
 !======================================================================
@@ -217 +218 @@
 write(*,*) "phyetat0: Surface emissivity <emis> range:", &
              minval(emis), maxval(emis)
+
+! Depth of methane tank (added by JVO 2017)
+ if (startphy_file) then
+   call get_field("tankCH4",tankCH4,found,indextime)
+   if (.not.found) then
+     write(*,*) "phyetat0: Failed loading <tankCH4>"
+     !  call abort
+     tankCH4(:)=2.
+   endif
+ else
+   tankCH4(:)=2.
+ endif ! of if (startphy_file)
+ write(*,*) "phyetat0: Depth of methane tank <tankCH4> range:", &
+              minval(tankCH4), maxval(tankCH4)
 
 ! pbl wind variance