Changeset 2902 for trunk/LMDZ.MARS/libf

Timestamp:

Feb 22, 2023, 6:56:01 PM (21 months ago)

Author:

llange

Message:

PEM
First implementation of a dynamic ice table following Schorgofer, Icarus (2010) / MSIM
The dynamic of the ice is implemented here. Impact on the soil properties (thermal, inertia) is currently implemented and will be commit later
Small fix for the ice table at equilibrium
LL

File:

• trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F (modified) (4 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----------------------------------------------------------------------
 
@@ -176 +177 @@
 c   --------------
 c
+       totflux(:) = fluxsurf_sw(:,1)+fluxsurf_sw(:,2)
       ! AS: firstcall OK absolute
       IF (firstcall) THEN
@@ -500 +502 @@
 !     ----------------------------------------
       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)
@@ -739 +753 @@
 !             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
+!             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