Changeset 2802 for trunk/LMDZ.GENERIC/libf

Timestamp:

Oct 19, 2022, 11:24:43 AM (2 years ago)

Author:

aslmd

Message:

Fixed initialisation of RH_generic &cloudfrac

File:

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

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

@@ -746 +746 @@
             call getin_p("corrkdir",corrkdir)
             print*,"corrkdir = ", corrkdir
-            write (tmp1, '(i3)') L_NSPECTI
-            write (tmp2, '(i3)') L_NSPECTV
+            write (tmp1, '(i4)') L_NSPECTI
+            write (tmp2, '(i4)') L_NSPECTV
             banddir=trim(trim(adjustl(tmp1))//'x'//trim(adjustl(tmp2)))
             banddir=trim(trim(adjustl(corrkdir))//'/'//trim(adjustl(banddir)))
@@ -1630 +1630 @@
 
                ! Let's loop on tracers
+               cloudfrac(:,:)=0.0
                do iq=1,nq
-
                   call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
-                  
                   if (call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
                      do l = 1, nlayer
                         do ig=1,ngrid
-                           cloudfrac(ig,l)=MAX(rneb_generic(ig,l,iq),rneb_man(ig,l)) 
+                           cloudfrac(ig,l)=rneb_generic(ig,l,iq)
                         enddo
                      enddo
@@ -1718 +1717 @@
          ! Sedimentation.
          if (sedimentation) then
-        
+
             zdqsed(1:ngrid,1:nlayer,1:nq) = 0.0
             zdqssed(1:ngrid,1:nq)  = 0.0
@@ -2156 +2155 @@
       ! Calculate RH_generic (Generic Relative Humidity) for diagnostic.
       if(generic_condensation)then
+         RH_generic(:,:,:)=0.0
          do iq=1,nq
 
@@ -2499 +2499 @@
             call writediagfi(ngrid,trim(noms(iq))//'_col',trim(noms(iq))//'_col',    & 
                             'kg m^-2',2,qcol(1,iq) )
-                          
 !          call writediagfi(ngrid,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  & 
 !                          'kg m^-2',2,qsurf(1,iq) )