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

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

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• phyetat0_mod.F90 (modified) (3 diffs)
• phyredem.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (4 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

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

@@ -135 +135 @@
 
 subroutine physdem1(filename,nsoil,ngrid,nlay,nq, &
-                    phystep,time,tsurf,tsoil,emis,q2,qsurf)
+                    phystep,time,tsurf,tsoil,emis,q2,qsurf,tankCH4)
   ! write time-dependent variable to restart file
   use iostart, only : open_restartphy, close_restartphy, & 
@@ -155 +155 @@
   real,intent(in) :: q2(ngrid,nlay+1)
   real,intent(in) :: qsurf(ngrid,nq)
+  real,intent(in) :: tankCH4(ngrid)
 
   integer :: iq
@@ -176 +177 @@
   ! Planetary Boundary Layer
   call put_field("q2","pbl wind variance",q2)
+
+  ! Methane tank depth
+  call put_field("tankCH4","Depth of methane tank",tankCH4)
   
 ! tracers

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

@@ -375 +375 @@
 !$OMP THREADPRIVATE(dycchi,qysat,nomqy)
 
+      real,dimension(:),allocatable,save :: tankCH4 ! Depth of surface methane tank
+!$OMP THREADPRIVATE(tankCH4)
+
 !-----------------------------------------------------------------------------
       
@@ -402 +405 @@
         ALLOCATE(qsurf(ngrid,nq))  
         ALLOCATE(q2(ngrid,nlayer+1)) 
+        ALLOCATE(tankCH4(ngrid)) 
         ALLOCATE(ztprevious(ngrid,nlayer))
         ALLOCATE(zuprevious(ngrid,nlayer))
@@ -448 +452 @@
 !        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
          call phyetat0(startphy_file,ngrid,nlayer,"startfi.nc",0,0,nsoilmx,nq,      &
-                       day_ini,time_phys,tsurf,tsoil,emis,q2,qsurf)                       
+                       day_ini,time_phys,tsurf,tsoil,emis,q2,qsurf,tankCH4)                       
          if (.not.startphy_file) then
            ! additionnal "academic" initialization of physics
@@ -1241 +1245 @@
             call physdem1("restartfi.nc",nsoilmx,ngrid,nlayer,nq, &
                           ptimestep,ztime_fin,                    &
-                          tsurf,tsoil,emis,q2,qsurf_hist)
+                          tsurf,tsoil,emis,q2,qsurf_hist,tankCH4)
          endif