Changeset 3753 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Apr 30, 2025, 5:55:07 PM (2 months ago)

Author:

afalco

Message:

Generic PCM

Tracers in 1D can now be initialized with VMR set in traceur.def
(homogeneous value or profile) with or without chemistry. Inichim_1D
is still called after. Initialization of co2_ice, h2o_vap and h2o_ice
can still be initialized with profile_co2_ice etc. the old way.

MM

Location:

trunk/LMDZ.GENERIC/libf/phystd/dyn1d

Files:

• inichim_1D.F90 (modified) (3 diffs)
• initracer_1D.F90 (added)
• rcm1d.F (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/inichim_1D.F90

@@ +1 @@
+!module inichim_1D_mod
+
+!implicit none
+
       subroutine inichim_1D(nlayer, nq, pq, qsurf, play, &
                                   flagh2o,flagthermo)
@@ -130 +134 @@
          else if (qxf(iq) /= 'None') then
             ! Opening file
-            fil = trim(datadir)//'/chemical_profiles/'//qxf(iq)
-            print*, 'chemical pofile '//trim(noms(iq))//': ', fil
+            fil = qxf(iq)
+            print*, 'chemical profile '//trim(noms(iq))//': ', fil
             open(UNIT=90,FILE=fil,STATUS='old',iostat=ierr)
             if (ierr.eq.0) then
@@ -206 +210 @@
 ! Or use #ModernTrac-v1 and add another option section 2.
 
-      end
+   end subroutine inichim_1D
+
+!end module inichim_1D_mod

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

@@ -44 +44 @@
      &                  NF90_CLOSE
       use version_info_mod, only: print_version_info
+      !use inichim_1D_mod, only: inichim_1D
+      !use initracer_1D_mod, only: initracer_1D
 
       implicit none
@@ -1094 +1096 @@
 ! -----------------
 #ifndef MESOSCALE
-      if(tracer.and.photochem.and. .not.restart) then
-           print*, "Calling inichim_1D"
+      if(tracer.and..not.restart) then
+           print*, "Initializing tracers properties and values"
            call initracer(1,nq)
            allocate(nametmp(nq))
            nametmp(1:nq)=tname(1:nq)
-           call inichim_1D(nlayer, nq, q, qsurf, play, 0, 0)
+           call initracer_1D(nlayer, nq, q, qsurf, play)
+           if(photochem) then
+                  print*, "Calling inichim_1D"
+                  call inichim_1D(nlayer, nq, q, qsurf, play, 0, 0)
+           end if
            tname(1:nq)=nametmp(1:nq)
            noms(1:nq)=nametmp(1:nq)