Changeset 2878 for trunk/LMDZ.MARS/libf/phymars/dyn1d

Timestamp:

Jan 25, 2023, 11:39:35 AM (22 months ago)

Author:

llange

Message:

Mars PCM 1D

Following r-2827, fixing bugs

1. igcm_co2 is now initialized
2. remove the abort_gcm in read_profile: The qsurf readed at this line

is not useful at it is initialized later

LL

Location:

trunk/LMDZ.MARS/libf/phymars/dyn1d

Files:

• read_profile_mod.F90 (modified) (2 diffs)
• testphys1d.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/dyn1d/read_profile_mod.F90

@@ -39 +39 @@
 
   use infotrac, only: tname
-
+  use tracer_mod, only: igcm_co2
   implicit none
 !---------------------------------------------------------------------------------------------------------------------!
@@ -110 +110 @@
     if (ierr.eq.0) then
       read(91,*)qsurf(iq)
-      if (trim(name_tracer(iq)) == 'co2' .and. qsurf(iq).ne.0) then
-        call abort_physic("read_profile","The value of co2 profile on surface must be 0. It's value must be specified in the run.def as co2ice=XX",1)
+      if (trim(name_tracer(iq)) == 'co2') then
+         igcm_co2 = iq
       endif
       do ilayer = 1, nb_layer

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -546 +546 @@
 c  CO2 ice on the surface
 c  -------------------
-      print *, "RVVV: igcm_co2", igcm_co2
+
       qsurf(igcm_co2)=0.E+0 ! default value for co2ice
       PRINT *,'Initial CO2 ice on the surface (kg.m-2)'