Changeset 2903

Timestamp:

Feb 22, 2023, 7:08:12 PM (21 months ago)

Author:

llange

Message:

PCM
Fix a bug in CO2condens. Surface temperature used at the surface as a boundary condition for the VL scheme must be ztcondsol. I've directly put ztcondsol rather that ptsurf+pdtsurfc*timestep.
Remove stuff about the solar albedo dependant flux to r2901, sorry about the mistake
LL & FF

File:

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

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

@@ -170 +170 @@
       REAL masseq(nlayer),wq(nlayer+1)
       INTEGER ifils,iq2
-c LL: flux dependant albedo
-      real :: totflux(ngrid)
 c----------------------------------------------------------------------
 
@@ -177 +175 @@
 c   --------------
 c
-       totflux(:) = fluxsurf_sw(:,1)+fluxsurf_sw(:,2)
       ! AS: firstcall OK absolute
       IF (firstcall) THEN
@@ -502 +499 @@
 !     ----------------------------------------
       CALL albedocaps(zls,ngrid,piceco2,psolaralb,emisref)
-      DO ig = 1,ngrid
-        if(piceco2(ig).gt.0.) then
-          psolaralb(ig,1) = 0.267+0.00059*totflux(ig)
-          psolaralb(ig,2) = psolaralb(ig,1)
-          if(ig.eq.ngrid) then
-           psolaralb(ig,1) = 0.532+8.72e-4*totflux(ig)
-           psolaralb(ig,2) = psolaralb(ig,1)
-          endif
-        endif
-      ENDDO
-
-
 
 ! set pemisurf() to emissiv when there is bare surface (needed for co2snow)
@@ -577 +562 @@
 c           Value transfert at the surface interface when condensation
 c           sublimation:
-            ztm(1) = ztsrf(ig) + pdtsrfc(ig)*ptimestep
+            ztm(1) = ztcondsol(ig)
             zum(1) = 0  
             zvm(1) = 0  
@@ -751 +736 @@
            ! NB: Updated surface pressure, at grid point 'ngrid', is
            !     ps(ngrid)=pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep
-!             write(*,*) "co2condens: South pole: latitude(ngrid)=",
-!     &                                           latitude(ngrid)
-!             ztcondsol(ngrid)=
-!     &          1./(bcond-acond*log(.01*vmr_co2(ngrid,1)*
-!     &                    (pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep)))
-!             pdtsrfc(ngrid)=(ztcondsol(ngrid)-ztsrf(ngrid))/ptimestep
+             write(*,*) "co2condens: South pole: latitude(ngrid)=",
+     &                                           latitude(ngrid)
+             ztcondsol(ngrid)=
+     &          1./(bcond-acond*log(.01*vmr_co2(ngrid,1)*
+     &                    (pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep)))
+             pdtsrfc(ngrid)=(ztcondsol(ngrid)-ztsrf(ngrid))/ptimestep
            endif
          endif