Changeset 1789

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

Files:

• README (modified) (1 diff)
• libf/dynphy_lonlat/phytitan/newstart.F (modified) (4 diffs)
• libf/dynphy_lonlat/phytitan/start2archive.F (modified) (6 diffs)
• libf/phytitan/phyetat0_mod.F90 (modified) (3 diffs)
• libf/phytitan/phyredem.F90 (modified) (3 diffs)
• libf/phytitan/physiq_mod.F90 (modified) (4 diffs)

trunk/LMDZ.TITAN/README

@@ -1333 +1333 @@
 Removed : aerave.F aeroptproperties.F90 radii_mod.F90 aerosol_mod.F90 aerave_new.F stokes.F90 
 vlz_fi.F callsedim.F iniaerosol.F suaer_corrk.F90 newsedim.F
+
+== 28/09/2017 == JVO
+Added the surface methane tank and put it in start files 

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 

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