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

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/dynphy_lonlat/phytitan

Files:

• newstart.F (modified) (4 diffs)
• start2archive.F (modified) (6 diffs)

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F

@@ -162 +162 @@
       real fact2
 
+!     added by JVO for Titan
+      REAL tankCH4(ngridmx)
 
 c sortie visu pour les champs dynamiques
@@ -341 +343 @@
         CALL phyetat0(.true.,ngridmx,llm,fichnom,tab0,Lmodif,nsoilmx,
      .        nqtot,day_ini,time,
-     .        tsurf,tsoil,emis,q2,qsurf)
+     .        tsurf,tsoil,emis,q2,qsurf,tankCH4)
 
         ! copy albedo and soil thermal inertia on (local) physics grid
@@ -1232 +1234 @@
  999  continue
 
-      
+
+c================================================================
+c   Initialisation for methane surface tank (added by JVO 2017)
+c================================================================
+      DO ig=1, ngridmx
+         tankCH4(ig)=2.0
+      ENDDO
 
 c=======================================================================
@@ -1339 +1347 @@
       call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
      &                dtphys,real(day_ini),
-     &                tsurf,tsoil,emis,q2,qsurf)
+     &                tsurf,tsoil,emis,q2,qsurf,tankCH4)
 
 c=======================================================================

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/start2archive.F

@@ -80 +80 @@
       INTEGER*4 day_ini_fi
 
+!     added by JVO for methane surface tank
+      REAL tankCH4(ngridmx)
 
 c Variable naturelle / grille scalaire
@@ -91 +93 @@
       REAL emisS(ip1jmp1)
 
+!     added by JVO for methane surface tank
+      REAL tankCH4S(ip1jmp1)
 
 c Variables intermediaires : vent naturel, mais pas coord scalaire
@@ -219 +223 @@
       CALL phyetat0(.true.,ngridmx,llm,fichnom,0,Lmodif,nsoilmx,nqtot,
      .      day_ini_fi,timefi,
-     .      tsurf,tsoil,emis,q2,qsurf)
+     .      tsurf,tsoil,emis,q2,qsurf,tankCH4)
 
 
@@ -310 +314 @@
 c q2 --> q2S
 c qsurf --> qsurfS
+c tankCH4 --> tankCH4S
 c
 c-----------------------------------------------------------------------
@@ -320 +325 @@
       call gr_fi_dyn(llm+1,ngridmx,iip1,jjp1,q2,q2S)
       call gr_fi_dyn(nqtot,ngridmx,iip1,jjp1,qsurf,qsurfS)
+      call gr_fi_dyn(1,ngridmx,iip1,jjp1,tankCH4,tankCH4S)
 
 c=======================================================================
@@ -458 +464 @@
 ! Note: no need to write volcapa, it is stored in "controle" table
 
+c-----------------------------------------------------------------
+c Ecriture du champs  tankCH4
+c-----------------------------------------------------------------
+      call write_archive(nid,ntime,'tankCH4',
+     &         'Depth of surface methane tank','m',2,tankCH4S)
 
 c Fin