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

Timestamp:

Mar 16, 2023, 2:56:08 PM (22 months ago)

Author:

emillour

Message:

Mars PCM:
Some corrections to the 1D model to take into account recent changes
(removing co2ice and added islope dimension to some arrays).
EM

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (5 diffs)

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

@@ -13 +13 @@
      &                     watercaptag, watercap, hmons, summit, base
       use slope_mod, only: theta_sl, psi_sl
+      use comslope_mod, only: def_slope,subslope_dist
       use phyredem, only: physdem0,physdem1
       use geometry_mod, only: init_geometry
@@ -486 +487 @@
       call getin("slope_orientation",psi_sl(1))
       
+      ! sub-slopes parameters (assuming no sub-slopes distribution for now). 
+      def_slope(1)=-90 ! minimum slope angle
+      def_slope(2)=90 ! maximum slope angle
+      subslope_dist(1,1)=1 ! fraction of subslopes in mesh 
 c
 c  for the gravity wave scheme
@@ -546 +551 @@
 c  CO2 ice on the surface
 c  -------------------
-
+      ! get the index of co2 tracer (not known at this stage)
+      igcm_co2=0
+      do iq=1,nq
+        if (trim(tname(iq))=="co2") then
+          igcm_co2=iq
+        endif
+      enddo
+      if (igcm_co2==0) then
+        write(*,*) "testphys1d error, missing co2 tracer!"
+        stop
+      endif
       qsurf(igcm_co2)=0.E+0 ! default value for co2ice
       PRINT *,'Initial CO2 ice on the surface (kg.m-2)'
@@ -656 +671 @@
 c Check if the surface is a water ice reservoir 
 c --------------------------------------------------
-      watercap(1)=0 ! Initialize watercap
+      watercap(1,:)=0 ! Initialize watercap
       watercaptag(1)=.false. ! Default: no water ice reservoir
       print *,'Water ice cap on ground ?'
@@ -700 +715 @@
       call physdem0("startfi.nc",longitude,latitude,nsoilmx,ngrid,llm,
      &              nq,dtphys,float(day0),0.,cell_area,
-     &              albedodat,inertiedat)
+     &              albedodat,inertiedat,def_slope,subslope_dist)
       call physdem1("startfi.nc",nsoilmx,ngrid,llm,nq,
      &              dtphys,time,