Changeset 2950

Timestamp:

Apr 27, 2023, 11:14:41 AM (19 months ago)

Author:

evos

Message:

updating evol_h2o_ice_routine to handle 1D cases

File:

• trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod_slope.F90 (modified) (6 diffs)

trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod_slope.F90

@@ -42 +42 @@
   pos_tend=0.
   neg_tend=0.
-
+if (ngrid.NE.1) then ! to make sure we are not in 1D 
   do i=1,ngrid
      do islope=1,nslope
@@ -54 +54 @@
      enddo
   enddo
-
-  if(neg_tend.GT.pos_tend .and. pos_tend.GT.0) then
+   if(neg_tend.GT.pos_tend .and. pos_tend.GT.0) then
      do i=1,ngrid
        do islope=1,nslope
@@ -65 +64 @@
        enddo
      enddo
-  elseif(neg_tend.LT.pos_tend .and. neg_tend.GT.0) then
+   elseif(neg_tend.LT.pos_tend .and. neg_tend.GT.0) then
      do i=1,ngrid
        do islope=1,nslope
@@ -75 +74 @@
        enddo
      enddo
-  elseif(pos_tend.EQ.0 .OR. neg_tend.EQ.0) then 
+   elseif(pos_tend.EQ.0 .OR. neg_tend.EQ.0) then 
     print *, "Reason of stopping : There is either no water ice sublimating or no water ice increasing !!"
     print *, "Tendencies on ice sublimating=", neg_tend
@@ -86 +85 @@
        enddo
     enddo
+   endif
   endif
-
   negative_part = 0.
 
@@ -115 +114 @@
     enddo
   enddo
-
-
+else ! ngrid==1;
+        qsurf(i,islope)=qsurf(i,islope)+new_tendencies(i,islope)*dt_pem
+endif
 
 END SUBROUTINE evol_h2o_ice_s_slope