Changeset 3266 for trunk/LMDZ.GENERIC/libf

Timestamp:

Mar 13, 2024, 5:13:20 PM (7 months ago)

Author:

mturbet

Message:

first cleanup in hydrol.F90

File:

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

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

@@ -55 +55 @@
       real snowlayer
       parameter (snowlayer=33.0)        ! 33 kg/m^2 of snow, equal to a layer of 3.3 cm 
-      real oceantime
+      
+      real oceantime ! this is a relaxation timescale for the oceanbulkavg parameterization
       parameter (oceantime=10*24*3600)
-
-      logical,save :: oceanbulkavg ! relax ocean temperatures to a GLOBAL mean value?
+      logical,save :: oceanbulkavg ! simple parameterization to relax ocean temperatures to the global mean value (crude, 0th order parameterization to mimick ocean heat transport)
+      
       logical,save :: activerunoff ! enable simple runoff scheme?
-      logical,save :: oceanalbvary ! ocean albedo varies with the diurnal cycle?
+      logical,save :: oceanalbvary ! simple parameterization to account for the effect of solar zenith angle on the ocean albedo (for the moment it is not used, but to be included in the future)
 !$OMP THREADPRIVATE(oceanbulkavg,activerunoff,oceanalbvary)
 
 !     Arguments
 !     ---------
-      real rnat(ngrid) ! I changed this to integer (RW)
+      real rnat(ngrid) ! rnat is terrain type: 0-ocean; 1-continent
       real,dimension(:),allocatable,save :: runoff
       real totalrunoff, tsea, oceanarea
@@ -175 +176 @@
       endif
 
-!       call writediagfi(ngrid,"hydrol_rnat"," "," ",2,rnat)
-!!       write (*,*) "oceanarea", oceanarea
+!      write (*,*) "oceanarea", oceanarea
 
 !     add physical tendencies already calculated
@@ -188 +188 @@
          enddo    
       enddo
-
-!      call writediagfi(ngrid,"hydrol_zqsurf_iliq_1"," "," ",2,zqsurf(1:ngrid,iliq))
  
       do ig=1,ngrid
@@ -350 +348 @@
       end do ! ig=1,ngrid
 
-!      call writediagfi(ngrid,"hydrol_zqsurf_iliq_2"," "," ",2,zqsurf(1:ngrid,iliq))
-
-!     perform crude bulk averaging of temperature in ocean
+
+!     simple parameterization to perform crude bulk averaging of temperature in ocean
 !     ----------------------------------------------------
-      if(oceanbulkavg)then
+      if(oceanbulkavg)then 
 
          oceanarea2=0.       
@@ -433 +430 @@
       endif !activerunoff         
 
-!      call writediagfi(ngrid,"hydrol_zqsurf_iliq_3"," "," ",2,zqsurf(1:ngrid,iliq))
-
 !     Re-add the albedo effects of CO2 ice if necessary
 !     -------------------------------------------------
@@ -455 +450 @@
       enddo
 
-!      call writediagfi(ngrid,"hydrol_dqs_hyd_iliq"," "," ",2,dqs_hyd(1:ngrid,iliq))  
-
       if (activerunoff) then
         call writediagfi(ngrid,'runoff','Runoff amount',' ',2,runoff)