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

Timestamp:

Nov 23, 2022, 11:13:18 AM (2 years ago)

Author:

romain.vande

Message:

Mars 1D PCM:
Same as previous commit in 1D.
The variable co2ice is deleted. All the co2 ice on surface is now in qsurf(:,igcm_co2).
CO2 tracer is now mandatory. diagfi output is unchanged.
RV

Location:

trunk/LMDZ.MARS/libf/phymars/dyn1d

Files:

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

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

@@ -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)
+      endif
       do ilayer = 1, nb_layer
         read(91,*)q(ilayer,iq)

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

@@ -16 +16 @@
       use geometry_mod, only: init_geometry
       use watersat_mod, only: watersat
-      use tracer_mod, only: igcm_h2o_vap,igcm_h2o_ice
+      use tracer_mod, only: igcm_h2o_vap,igcm_h2o_ice,igcm_co2
       use planete_h, only: year_day, periheli, aphelie, peri_day,
      &                     obliquit, emin_turb, lmixmin
@@ -95 +95 @@
       REAL,ALLOCATABLE :: qsurf(:) ! tracer surface budget (e.g. kg.m-2)
       REAL tsoil(nsoilmx)   ! subsurface soik temperature (K)
-      REAL co2ice(1)        ! co2ice layer (kg.m-2)
       REAL emis(1)          ! surface layer
       REAL albedo(1,1)      ! surface albedo
@@ -547 +546 @@
 c  CO2 ice on the surface
 c  -------------------
-      co2ice(1)=0.E+0 ! default value for co2ice
+      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)'
-      call getin("co2ice",co2ice)
-      write(*,*) " co2ice = ",co2ice
+      call getin("co2ice",qsurf(igcm_co2))
+      write(*,*) " co2ice = ",qsurf(igcm_co2)
 
 c
@@ -556 +556 @@
 c  ----------
       emis=emissiv
-      IF (co2ice(1).eq.1.E+0) THEN
+      IF (qsurf(igcm_co2).eq.1.E+0) THEN
          emis=emisice(1) ! northern hemisphere
          IF(latitude(1).LT.0) emis=emisice(2) ! southern hemisphere
@@ -704 +704 @@
       call physdem1("startfi.nc",nsoilmx,ngrid,llm,nq,
      &              dtphys,time,
-     &              tsurf,tsoil,co2ice,albedo,emis,q2,qsurf,tauscaling,
+     &              tsurf,tsoil,albedo,emis,q2,qsurf,tauscaling,
      &              totcloudfrac,wstar,watercap)