Changeset 1790 for trunk/LMDZ.TITAN/libf/phytitan/dyn1d/rcm1d.F

Timestamp:

Sep 29, 2017, 1:57:44 PM (7 years ago)

Author:

jvatant

Message:

Update 1D model that hasn't been cleaned since migration from generic

+ remove co2,water, kcm ...
+ added moyzon

Also correct deftank according to r1788
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/dyn1d/rcm1d.F (modified) (13 diffs)

trunk/LMDZ.TITAN/libf/phytitan/dyn1d/rcm1d.F

@@ -6 +6 @@
       use mod_grid_phy_lmdz, only : regular_lonlat
       use infotrac, only: nqtot, tname
-      use surfdat_h, only: albedodat, phisfi, dryness, watercaptag,
+      use surfdat_h, only: albedodat, phisfi,
      &                     zmea, zstd, zsig, zgam, zthe,
      &                     emissiv, emisice, iceradius,
@@ -14 +14 @@
       use phyredem, only: physdem0,physdem1
       use geometry_mod, only: init_geometry
-      use slab_ice_h, only: noceanmx
       use planete_mod, only: apoastr,periastr,year_day,peri_day,
      &         obliquit,nres,z0,lmixmin,emin_turb,coefvis,coefir,
@@ -23 +22 @@
      &                            nday, iphysiq
       use callkeys_mod, only: tracer,check_cpp_match,rings_shadow,
-     &                        specOLR,water,pceil,ok_slab_ocean
+     &                        specOLR,pceil
       USE comvert_mod, ONLY: ap,bp,aps,bps,pa,preff, sig,
      &                       presnivs,pseudoalt,scaleheight
@@ -58 +57 @@
 !     A. Spiga
 !     J. Leconte (2012)
+!     J. Vatant d'Ollone (2017)
 !
 !==================================================================
@@ -88 +88 @@
       REAL,ALLOCATABLE :: qsurf(:)    ! tracer surface budget (e.g. kg.m-2)
       REAL,ALLOCATABLE :: tsoil(:)    ! subsurface soil temperature (K)
-!      REAL co2ice               ! co2ice layer (kg.m-2) !not used anymore
-      integer :: i_co2_ice=0     ! tracer index of co2 ice
-      integer :: i_h2o_ice=0     ! tracer index of h2o ice
-      integer :: i_h2o_vap=0     ! tracer index of h2o vapor
       REAL emis(1)               ! surface layer
       REAL q2(llm+1)             ! Turbulent Kinetic Energy
@@ -112 +108 @@
       integer :: nq !=1 ! number of tracers
  
-      character*2 str2
-      character (len=7) :: str7
       character(len=44) :: txt
 
-      logical oldcompare, earthhack,saveprofile
+      logical saveprofile
 
 !     added by RW for zlay computation
@@ -125 +119 @@
       real logplevs(llm)
 
-!     added by BC
-      REAL cloudfrac(1,llm)
-      REAL hice(1),totcloudfrac(1)
-
-!     added by BC for ocean
-      real rnat(1)
-      REAL tslab(1,noceanmx),tsea_ice(1),sea_ice(1)
-      real pctsrf_sic(1)
-
-
+!     added by JVO
+      REAL tankCH4(1)
 
 !     added by AS to avoid the use of adv trac common
@@ -157 +143 @@
       IF (.not. ALLOCATED(zgam)) ALLOCATE(zgam(1))
       IF (.not. ALLOCATED(zthe)) ALLOCATE(zthe(1))
-      IF (.not. ALLOCATED(dryness)) ALLOCATE(dryness(1))
-      IF (.not. ALLOCATED(watercaptag)) ALLOCATE(watercaptag(1))
       !! those are defined in comdiurn_h.F90
       IF (.not.ALLOCATED(sinlat)) ALLOCATE(sinlat(1))
@@ -272 +256 @@
             endif
           enddo !of do iq=1,nq
-! check for co2_ice / h2o_ice tracers:
-         i_co2_ice=0
-         i_h2o_ice=0
-         i_h2o_vap=0
-         do iq=1,nq
-           if (tname(iq)=="co2_ice") then
-             i_co2_ice=iq
-           elseif (tname(iq)=="h2o_ice") then
-             i_h2o_ice=iq
-           elseif (tname(iq)=="h2o_vap") then
-             i_h2o_vap=iq
-           endif
-         enddo
         else
           write(*,*) 'Cannot find required file "traceur.def"'
@@ -577 +548 @@
             write(*,*)"  tracer:",trim(txt)
               
-            ! CO2
-            if (txt.eq."co2_ice") then
-               q(:,iq)=0.   ! kg/kg of atmosphere
-               qsurf(iq)=0. ! kg/m2 at the surface               
-               ! Look for a "profile_co2_ice" input file
-               open(91,file='profile_co2_ice',status='old',
-     &         form='formatted',iostat=ierr)
-               if (ierr.eq.0) then
-                  read(91,*) qsurf(iq)
-                  do ilayer=1,nlayer
-                     read(91,*) q(ilayer,iq)
-                  enddo
-               else
-                  write(*,*) "No profile_co2_ice file!"
-               endif
-               close(91)
-            endif ! of if (txt.eq."co2")
-          
-            ! WATER VAPOUR
-            if (txt.eq."h2o_vap") then
-               q(:,iq)=0.   ! kg/kg of atmosphere
-               qsurf(iq)=0. ! kg/m2 at the surface
-               ! Look for a "profile_h2o_vap" input file   
-               open(91,file='profile_h2o_vap',status='old',
-     &         form='formatted',iostat=ierr)
-               if (ierr.eq.0) then
-                  read(91,*) qsurf(iq)
-                  do ilayer=1,nlayer
-                     read(91,*) q(ilayer,iq)
-                  enddo
-               else
-                  write(*,*) "No profile_h2o_vap file!"
-               endif
-               close(91)
-            endif ! of if (txt.eq."h2o_vap")
-            
-            ! WATER ICE
-            if (txt.eq."h2o_ice") then
-               q(:,iq)=0.   ! kg/kg of atmosphere
-               qsurf(iq)=0. ! kg/m2 at the surface
-               ! Look for a "profile_h2o_ice" input file
-               open(91,file='profile_h2o_ice',status='old',
-     &         form='formatted',iostat=ierr)
-               if (ierr.eq.0) then
-                  read(91,*) qsurf(iq)
-                  do ilayer=1,nlayer
-                     read(91,*) q(ilayer,iq)
-                  enddo
-               else
-                  write(*,*) "No profile_h2o_ice file!"
-               endif
-               close(91)
-            endif ! of if (txt.eq."h2o_ice")
-
          enddo ! of do iq=1,nq
          
@@ -673 +590 @@
       emissiv=emis(1) ! we do this so that condense_co2 sets things to the right
                    ! value if there is no snow
-
-      if(i_co2_ice.gt.0)then
-         qsurf(i_co2_ice)=0 ! default value for co2ice
-         print*,'Initial CO2 ice on the surface (kg.m-2)'
-         call getin("co2ice",qsurf(i_co2_ice))
-         write(*,*) " co2ice = ",qsurf(i_co2_ice)
-         IF (qsurf(i_co2_ice).ge.1.E+0) THEN
-            ! if we have some CO2 ice on the surface, change emissivity
-            if (latitude(1).ge.0) then ! northern hemisphere
-              emis(1)=emisice(1)
-            else ! southern hemisphere
-              emis(1)=emisice(2)
-            endif
-         ENDIF
-      endif
 
 c  calcul des pressions et altitudes en utilisant les niveaux sigma
@@ -823 +725 @@
 
 
-! Initialize slab ocean
-! -----------------
-      rnat=1. ! default value for rnat
-      if(inertiedat(1,1).GE.10000.)then
-         rnat=0.
-      endif
-      if(ok_slab_ocean)then
-      rnat=0.
-      tslab(1,1)=tsurf(1)
-      tslab(1,2)=tsurf(1)
-      tsea_ice=tsurf
-      pctsrf_sic=0.
-      sea_ice=0.
-      endif
-
-
-
 c  Write a "startfi" file
 c  --------------------
@@ -850 +735 @@
       call physdem1("startfi.nc",nsoilmx,1,llm,nq,
      &                dtphys,time,
-     &                tsurf,tsoil,emis,q2,qsurf,
-     &                cloudfrac,totcloudfrac,hice, 
-     &                rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
+     &                tsurf,tsoil,emis,q2,qsurf,tankCH4)
 
 c=======================================================================