Changeset 1802 for trunk/LMDZ.GENERIC/libf/phystd/dyn1d

Timestamp:

Oct 19, 2017, 11:27:15 PM (7 years ago)

Author:

bclmd

Message:

Adding photochemistry to LMDZ Generic

File:

• trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F (modified) (8 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -6 +6 @@
       use mod_grid_phy_lmdz, only : regular_lonlat
       use infotrac, only: nqtot, tname
+      use tracer_h, only: noms
       use surfdat_h, only: albedodat, phisfi, dryness, watercaptag,
      &                     zmea, zstd, zsig, zgam, zthe,
@@ -23 +24 @@
      &                            nday, iphysiq
       use callkeys_mod, only: tracer,check_cpp_match,rings_shadow,
-     &                        specOLR,water,pceil,ok_slab_ocean
+     &                        specOLR,water,pceil,ok_slab_ocean,photochem
       USE comvert_mod, ONLY: ap,bp,aps,bps,pa,preff, sig,
      &                       presnivs,pseudoalt,scaleheight
@@ -138 +139 @@
 !     added by AS to avoid the use of adv trac common
       character*20,allocatable :: nametrac(:)   ! name of the tracer (no need for adv trac common)
+      character*20,allocatable :: nametmp(:)   !
 
       real :: latitude(1), longitude(1), cell_area(1)
@@ -254 +256 @@
           if (nq>0) then
             allocate(tname(nq))
+            allocate(noms(nq))
             allocate(q(llm,nq))
             allocate(qsurf(nq))
@@ -267 +270 @@
           ! minimal version, just read in the tracer names, 1 per line
             read(90,*,iostat=ierr) tname(iq)
+            noms(iq)=tname(iq)
             if (ierr.ne.0) then
               write(*,*) 'rcm1d: error reading tracer names...'
@@ -292 +296 @@
         endif
         close(90)
+
 
       else ! of if (tracer)
@@ -816 +821 @@
 ! -----------------
       do isoil=0,nsoil-1
-        mlayer(isoil)=2.e-4*(2.**(isoil-0.5)) ! mid-layer depth
+        mlayer(isoil)=3.e-2*(2.**(isoil-0.5)) ! mid-layer depth
       enddo
       do isoil=1,nsoil
-        layer(isoil)=2.e-4*(2.**(isoil-1)) ! layer depth
+        layer(isoil)=3.e-2*(2.**(isoil-1)) ! layer depth
       enddo
 
@@ -838 +843 @@
       endif
 
+
+! Initialize chemical species
+! -----------------
+      if(tracer.and.photochem) then
+           call initracer(1,nq,tname)
+           allocate(nametmp(nq))
+           nametmp(1:nq)=tname(1:nq)
+           call inichim_1D(nq, q, qsurf, psurf, 0, 0)
+           tname(1:nq)=nametmp(1:nq)
+           noms(1:nq)=nametmp(1:nq)
+      endif ! tracer and photochem