Changeset 3267 for trunk

Timestamp:

Mar 13, 2024, 5:50:44 PM (8 months ago)

Author:

mturbet

Message:

second cleanup in hydrol.F90

File:

• trunk/LMDZ.GENERIC/libf/phystd/hydrol.F90 (modified) (15 diffs)

trunk/LMDZ.GENERIC/libf/phystd/hydrol.F90

@@ -28 +28 @@
 !     -------
 !     Calculate the surface hydrology and albedo changes.
+!     Both for oceanic and continental regions
 !     
 !     Authors
@@ -109 +110 @@
 !$OMP THREADPRIVATE(firstcall)
 
-     real :: runoffamount(ngrid)
+      real :: runoffamount(ngrid)
 !#ifdef CPP_PARA
       real :: runoffamount_glo(klon_glo)
@@ -130 +131 @@
          call getin_p("activerunoff",activerunoff)
          write(*,*)" activerunoff = ",activerunoff
-         
-         
          
          if (activerunoff) then
@@ -197 +196 @@
       do ig = 1, ngrid
 
-!     Ocean
-!     -----
+!     Ocean regions (rnat = 0)
+!     -----------------------
          if(nint(rnat(ig)).eq.0)then
 
-!     re-calculate oceanic albedo
+!     Parameterization (not used for the moment) to compute the effect of solar zenith angle on the albedo
+!     --------------------------
+!
 !            if(diurnal.and.oceanalbvary)then
 !               fauxo      = ( 1.47 - ACOS( mu0(ig) ) )/0.15 ! where does this come from (Benjamin)?
@@ -207 +208 @@
 !               albedo(ig) = MAX(MIN(albedo(ig),0.60),0.04)
 !            else
-               do nw=1,L_NSPECTV
-                  albedo(ig,nw) = alb_ocean ! For now, alb_ocean is defined in inifis_mod.F90. Later we could introduce spectral dependency for alb_ocean.
-               enddo
+!
+!               do nw=1,L_NSPECTV
+!                  albedo(ig,nw) = alb_ocean ! For now, alb_ocean is defined in inifis_mod.F90. Later we could introduce spectral dependency for alb_ocean.
+!              enddo
 !            end if
 
-
-            if(ok_slab_ocean) then
+            ! we first start by fixing the albedo of oceanic grid to that of the ocean
+            do nw=1,L_NSPECTV
+               albedo(ig,nw) = alb_ocean ! For now, alb_ocean is defined in inifis_mod.F90. Later we could introduce spectral dependency for alb_ocean.
+            enddo
+
+            if(ok_slab_ocean) then ! if ocean heat transport param activated
           
                zfra = MAX(0.0,MIN(1.0,zqsurf(ig,iice)/45.0))     ! Snow Fraction (Critical height 45kg/m2~15cm)
@@ -227 +233 @@
                ! Oceanic ice height, just for diagnostics
                hice(ig)    = MIN(10.,sea_ice(ig)/rhowater)
-            else !ok_slab_ocean
-
-
-!     calculate oceanic ice height including the latent heat of ice formation
-!     hice is the height of oceanic ice with a maximum of maxicethick.
+            else !ok_slab_ocean ; here this is the case where we are dealing with a static ocean
+
+
+               ! calculate oceanic ice height including the latent heat of ice formation
+               ! hice is the height of oceanic ice with a maximum of maxicethick.
                hice(ig)    = zqsurf(ig,iice)/rhowater ! update hice to include recent snowfall
-
-!              twater(ig)  = tsurf(ig) + ptimestep*zdtsurf(ig) &
-               twater(ig)  = ztsurf(ig) - hice(ig)*RLFTT*rhowater/pcapcal(ig)
-               ! this is temperature water would have if we melted entire ocean ice layer
+               twater(ig)  = ztsurf(ig) - hice(ig)*RLFTT*rhowater/pcapcal(ig) ! this is the temperature water would have if we melted the entire ocean ice layer
                hicebis(ig) = hice(ig)
                hice(ig)    = 0.
@@ -272 +275 @@
 
 
-!     Continent
-!     ---------
+!     Continental regions (rnat = 1)
+!     -----------------------
          elseif (nint(rnat(ig)).eq.1) then
 
-!     melt the snow
+            ! melt the snow
             if(ztsurf(ig).gt.T_h2O_ice_liq)then
                if(zqsurf(ig,iice).gt.1.0e-8)then
@@ -287 +290 @@
                   zqsurf(ig,iliq) = zqsurf(ig,iliq) + fsnoi
 
-!                 thermal effects
+                  ! thermal effects
                   pdtsurf_hyd(ig) = -fsnoi*RLFTT/pcapcal(ig)/ptimestep  
 
@@ -293 +296 @@
             else
 
-!     freeze the water
+               ! freeze the water
                if(zqsurf(ig,iliq).gt.1.0e-8)then
 
@@ -304 +307 @@
                   zqsurf(ig,iliq) = zqsurf(ig,iliq) - fsnoi
 
-!                 thermal effects
+                  ! thermal effects
                   pdtsurf_hyd(ig) = +fsnoi*RLFTT/pcapcal(ig)/ptimestep  
 
@@ -310 +313 @@
             endif
             
-!     deal with runoff
+            ! add runoff (to simulate transport of water from continental regions to oceanic regions ; in practice, this prevents liquid water to build up too much on continental regions)
             if(activerunoff)then
                runoff(ig) = max(zqsurf(ig,iliq) - mx_eau_sol, 0.0)
@@ -316 +319 @@
                   if(runoff(ig).ne.0.0)then
                      zqsurf(ig,iliq) = mx_eau_sol
-!                    runoff is added to ocean at end
+                    ! note: runoff is added to ocean at end
                   endif
                end if
@@ -322 +325 @@
             endif
 
-!     re-calculate continental albedo
+            ! re-calculate continental albedo
             DO nw=1,L_NSPECTV
                albedo(ig,nw) = albedo_bareground(ig)
@@ -389 +392 @@
          enddo
 
-!    collect on the full grid
+        ! collect on the full grid
         call gather(runoffamount,runoffamount_glo)
         call gather(zqsurf(1:ngrid,iliq),zqsurf_iliq_glo)
@@ -416 +419 @@
         endif! is_master
 
-!    scatter the field back on all processes
+        ! scatter the field back on all processes
         call scatter(zqsurf_iliq_glo,zqsurf(1:ngrid,iliq))