Changeset 3132

Timestamp:

Nov 21, 2023, 8:25:02 AM (20 months ago)

Author:

emillour

Message:

Mars PCM:
Bug fix in phyetat0, igcm_co2_tmp was not initialed when paleoclimate==.false.
but used nonetheless. Added an if(paleoclimate) around computation of
perennial_co2ice() as it should only be computed in that context; also added
an if(paleoclimate) around writting perennial_co2ice() to restartfi.nc in
physdem1 for the same reason.
EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/phyetat0_mod.F90 (modified) (3 diffs)
• libf/phymars/phyredem.F90 (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -4349 +4349 @@
 Correction of a bug related to r3126: 'choice_ads' did not have a default value which made the model crash. 'choice_ads = 0' is the default value (no adsorption).
 Cleaning of the 1D initialization: any reference of the PEM has been removed from "init_testphys1D_mod.F90". This way is much cleaner even though it needs more code.
+
+== 21/11/2023 == EM
+Bug fix in phyetat0, igcm_co2_tmp was not initialed when paleoclimate==.false.
+but used nonetheless. Added an if(paleoclimate) around computation of 
+perennial_co2ice() as it should only be computed in that context; also added
+an if(paleoclimate) around writting perennial_co2ice() to restartfi.nc in
+physdem1 for the same reason.

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -781 +781 @@
 ! Perennial CO2 ice
 perennial_co2ice(:,:) = 0.
-if (startphy_file) then
+if (paleoclimate) then
+  if (startphy_file) then
     call get_field("perennial_co2ice",perennial_co2ice,found,indextime)
     if (.not. found) then
@@ -793 +794 @@
         endif
     endif ! not found
-else
+  else
     if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
         do islope = 1,nslope
@@ -800 +801 @@
         enddo
     endif
-endif !startphy_file
+  endif !startphy_file
+endif ! of if (paleoclimate)
 
 ! close file:

trunk/LMDZ.MARS/libf/phymars/phyredem.F90

@@ -227 +227 @@
  
   ! Perennial CO2 ice layer
-  call put_field("perennial_co2ice","CO2 ice cover",perennial_co2ice,time)
+  if(paleoclimate) then
+    call put_field("perennial_co2ice","CO2 ice cover",perennial_co2ice,time)
+  endif
 
   ! Surface temperature