Changeset 3013 for trunk

Timestamp:

Jul 24, 2023, 10:40:53 PM (16 months ago)

Author:

llange

Message:

MARS PCM
CO2 condens: Fixing bug introduced in r2977: as the cooling of the air before condensing is now done in the surface condensation loop, the boundary conditions for the VL scheme must be changed in order to avoid duplicating the cooling.

Also fixed a bug in surfini.F (2892): qsurf is now ngrid x nq x nslope: I added the dimension nslope.

LL

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• co2condens_mod.F (modified) (1 diff)
• physiq_mod.F (modified) (1 diff)
• surfini.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F

@@ -647 +647 @@
 c           Value transfert at the surface interface when condensation
 c           sublimation:
-            ztm(1) = ztcondsol(ig)
+            if (zmflux(1).lt.0) then 
+c               Surface condensation
+                ztm(1) = ztc(1)
+            else 
+c               Surface sublimation:
+               ztm(1) =  ztcondsol(ig)
+            endif
             zum(1) = 0  
             zvm(1) = 0  

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -694 +694 @@
 c        Initialize albedo and orbital calculation
 c        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-         CALL surfini(ngrid,qsurf)
+         CALL surfini(ngrid,nslope,qsurf)
          CALL iniorbit(aphelie,periheli,year_day,peri_day,obliquit)
 c        initialize soil 

trunk/LMDZ.MARS/libf/phymars/surfini.F

@@ -1 @@
-      SUBROUTINE surfini(ngrid,qsurf)
+      SUBROUTINE surfini(ngrid,nslope,qsurf)
 
       USE ioipsl_getin_p_mod, ONLY : getin_p
@@ -28 +28 @@
 
       integer,intent(in) :: ngrid ! number of atmospheric columns
-      real,intent(inout) :: qsurf(ngrid,nqmx) ! tracer on surface (kg/m2)
-
-      INTEGER ig,icap,iq,alternate
+      integer,intent(in) :: nslope ! number of sub-grid surfaces
+      real,intent(inout) :: qsurf(ngrid,nqmx,nslope) ! tracer on surface (kg/m2)
+
+      INTEGER ig,icap,iq,alternate,islope
       REAL icedryness ! ice dryness
       
@@ -78 +79 @@
          !write(*,*) "all qsurf to zero. dirty."
          do iq=1,nqmx
-         qsurf(ig,iq)=0.  !! on jette les inputs GCM
+            do islope = 1,nslope
+              qsurf(ig,iq,islope)=0.  !! on jette les inputs GCM
                           !! on regle juste watercaptag
                           !! il faudrait garder les inputs GCM 
                           !! si elles sont consequentes
+            enddo
          enddo
          if ( ( latitude(ig)*180./pi .gt. 70. ) .and.