Changeset 3073 for trunk/LMDZ.GENERIC

Timestamp:

Oct 5, 2023, 2:56:44 PM (14 months ago)

Author:

yjaziri

Message:

Generic PCM:

Chemistry: correction missing initialisation
Add escape values in outputs

YJ

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/aeronostd/calchim_asis.F90 (modified) (3 diffs)
• libf/dynphy_lonlat/phystd/inichim_newstart.F90 (modified) (1 diff)
• libf/phystd/dyn1d/inichim_1D.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/README

@@ -1820 +1820 @@
 == 23/07/2023 == YJ
 Chemistry: correction in photolysis online with level instead of layer
+
+== 04/10/2023 == YJ
+Chemistry: correction missing initialisation
+Add escape values in outputs

trunk/LMDZ.GENERIC/libf/aeronostd/calchim_asis.F90

@@ -480 +480 @@
                   write(*,*) ' Parameter value (default)     : SF_mode=', SF_mode(iesp)
                end if
-               if (index(tracline,'SF_value='    ) /= 0) then
+               if (index(tracline,'SF_value='   ) /= 0) then
                   read(tracline(index(tracline,'SF_value=')+len('SF_value='):),*) SF_value(iesp)
                   write(*,*) ' Parameter value (traceur.def) : SF_value=', SF_value(iesp)
@@ -486 +486 @@
                   write(*,*) ' Parameter value (default)     : SF_value=', SF_value(iesp)
                end if
-               if (index(tracline,'prod_rate='    ) /= 0) then
+               if (index(tracline,'prod_rate='  ) /= 0) then
                   read(tracline(index(tracline,'prod_rate=')+len('prod_rate='):),*) prod_rate(iesp)
                   write(*,*) ' Parameter value (traceur.def) : prod_rate=', prod_rate(iesp)
@@ -492 +492 @@
                   write(*,*) ' Parameter value (default)     : prod_rate=', prod_rate(iesp)
                end if
+               if (index(tracline,'escape='     ) /= 0) then
+                  read(tracline(index(tracline,'escape=')+len('escape='):),*) escape(:,iesp)
+                  write(*,*) ' Parameter value (traceur.def) : escape=', escape(:,iesp)
+               else
+                  write(*,*) ' Parameter value (default)     : escape=', escape(:,iesp)
+               end if
             end if
          end do

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/inichim_newstart.F90

@@ -76 +76 @@
       qx(:)        = 0.
       qxf(:)       = 'None'
+      nlines       = 0
 
 ! 2. load in traceur.def chemistry data for initialization:

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

@@ -74 +74 @@
       qx(:)    = 0.
       qxf(:)   = 'None'
+      nlines   = 0
 
 ! 2. load in traceur.def chemistry data for initialization: