Changeset 1789 for trunk/LMDZ.TITAN/libf/phytitan/physiq_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/physiq_mod.F90 (modified) (4 diffs)

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