Changeset 2696

Timestamp:

Jun 21, 2022, 11:05:33 AM (3 years ago)

Author:

aslmd

Message:

fixed dqevap=0 before iterative process

File:

• trunk/LMDZ.GENERIC/libf/phystd/largescale_generic_mod.F90 (modified) (3 diffs)

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

@@ -111 +111 @@
                      !      zt(i)=15.   ! check too low temperatures
                 endif
-                call Psat_generic(zt(i),local_p,metallicity,psat_tmp,zqs_temp)
-                zqs(i)=zqs_temp
+                call Psat_generic(zt(i),local_p,metallicity,psat_tmp,zqs_temp) !useless ? 
+                zqs(i)=zqs_temp ! useless ?
                 zt(i)=pt(i,k)+pdt(i,k)*ptimestep
                 zx_q(i) = pq(i,k,igcm_generic_vap)+pdq(i,k,igcm_generic_vap)*ptimestep
-                dqevap(i,k)=0.
+                ! dqevap(i,k)=0.
         !           iterative process to stabilize the scheme when large water amounts JL12
                 zcond(i) = 0.0d0
@@ -131 +131 @@
                 End do ! niter
                 zcond(i)=MAX(zcond(i),-(pq(i,k,igcm_generic_ice)+pdq(i,k,igcm_generic_ice)*ptimestep))
+                zcond(i) = zcond(i)/ptimestep
             ENDDO ! i=1,ngrid
 
@@ -139 +140 @@
 
         Enddo ! k= nlayer, 1, -1
-
+        print*, "pdqvaplsc = ",pdqvaplsc(:,:)
+        print*,"pdtlsc = ",pdtlsc(:,:)
     end subroutine largescale_generic
 end module largescale_generic_mod