Changeset 3102

Timestamp:

Oct 25, 2023, 1:46:14 PM (13 months ago)

Author:

emillour

Message:

comments on rain_generic

File:

• trunk/LMDZ.GENERIC/libf/phystd/rain_generic.F90 (modified) (5 diffs)

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

@@ -298 +298 @@
                      ! JL22. Accounts for gravitational energy of falling precipitations (probably not to be used in the GCM
                      !   where the counterpart is not included in the dynamics.)
+                     ! only for mesoscale simulations !!!
 #endif
                   endif ! of if (precip_rate(i) .GT.0.)
@@ -307 +308 @@
 
             if(precip_scheme_generic.eq.1)then
+            ! for now this is the only precipitation scheme working here and is therefore the default one (NC2023)
 
                do i = 1, ngrid
@@ -322 +324 @@
 
             elseif (precip_scheme_generic.ge.2) then
-            
+            ! all the precipitation schemes below have not been tested and still have water-dependant routines
+            ! these schemes should be tested, improved ...
+            write(*,*) "Be careful, you have chosen a precip_scheme_generic equal to or greater than 2. For now it has not been tested."
                do i = 1, ngrid
                   if (rneb(i,k).GT.0.0) then
@@ -412 +416 @@
                enddo
 
-            endif ! if precip_scheme_generic=1
+            endif ! if precip_scheme_generic
 
          enddo ! of do k = nlayer, 1, -1
@@ -422 +426 @@
                call abort
             endif
-            if (t(i,1) .LT. T_h2O_ice_liq) then
-               dqssnow_generic(i,i_ice_generic) = precip_rate(i)
-               dqsrain_generic(i,i_ice_generic) = 0.0
-            else
-               dqssnow_generic(i,i_ice_generic) = 0.0
-               dqsrain_generic(i,i_ice_generic) = precip_rate(i) ! liquid water = ice for now
-            endif
-
-            ! For now we force :
+            ! lines below should be uncommented and customized if you want to make a difference between rain and snow
+            !if (t(i,1) .LT. T_h2O_ice_liq) then
+            !   dqssnow_generic(i,i_ice_generic) = precip_rate(i)
+            !   dqsrain_generic(i,i_ice_generic) = 0.0
+            !else
+            !   dqssnow_generic(i,i_ice_generic) = 0.0
+            !   dqsrain_generic(i,i_ice_generic) = precip_rate(i) ! liquid water = ice for now
+            !endif
+
+            ! For now we make no distinction between rain and snow, and consider only rain:
             dqsrain_generic(i,i_ice_generic) = precip_rate(i)
             dqssnow_generic(i,i_ice_generic) = 0.0