Changeset 3353 for trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

Timestamp:

May 31, 2024, 10:03:40 PM (6 months ago)

Author:

afalco

Message:

Pluto PCM:
Added zrecast & old sedim ;
Choose haze file ;
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (8 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -55 +55 @@
                               n2cond,nearn2cond,noseason_day,conservn2, &
                               convergeps,kbo,triton,paleo,paleoyears, &
-                              carbox, methane, oldplutovdifc, oldplutocorrk, &
-                              aerohaze,haze_proffix,source_haze, &
+                              carbox, methane,&
+                              oldplutovdifc,oldplutocorrk,oldplutosedim, &
+                              aerohaze,haze_proffix,source_haze,&
                               season, sedimentation,generic_condensation, &
                               specOLR, &
@@ -979 +980 @@
                                cloudfrac,totcloudfrac,.false.,                   &
                                firstcall,lastcall)
+                  albedo_equivalent(1:ngrid)=albedo(1:ngrid,1)
                else
                 call callcorrk(ngrid,nlayer,pq,nq,qsurf,  &
@@ -1480 +1482 @@
             end if
 
-            ! if (.not. water) then
-               ! Compute GCS (Generic Condensable Specie) cloud fraction. For now we can not have both water cloud fraction and GCS cloud fraction
-               ! Water is the priority
-               ! If you have set water and generic_condensation, then cloudfrac will be water cloudfrac
-               !
-               ! If you have set generic_condensation (and not water) and you have set several GCS
-               ! then cloudfrac will be only the cloudfrac of the last generic tracer
-               ! (Because it is rewritten every tracer in the loop)
-               !
-               ! Maybe one should create a cloudfrac_generic(ngrid,nlayer,nq) with 3 dimensions, the last one for tracers
-
                ! Let's loop on tracers
                cloudfrac(:,:)=0.0
@@ -1507 +1498 @@
          endif !generic_condensation
 
-         !Generic Rain !AF24: removed
-
          ! Sedimentation.
          if (sedimentation) then
@@ -1515 +1504 @@
             zdqssed(1:ngrid,1:nq)  = 0.0
 
-            ! if(watertest)then !AF24: removed
-
-            call callsedim(ngrid,nlayer,ptimestep,           &
+            if (oldplutosedim)then
+               call callsedim_pluto(ngrid,nlayer,ptimestep,    &
+                          pplev,zzlev,pt,pdt,rice_ch4,rice_co, &
+                          pq,pdq,zdqsed,zdqssed,nq,pphi)
+
+            else
+               call callsedim(ngrid,nlayer,ptimestep,       &
                           pplev,zzlev,pt,pdt,pq,pdq,        &
                           zdqsed,zdqssed,nq)
-
-            ! if(watertest)then !AF24: removed
+            endif
 
             ! Whether it falls as rain or snow depends only on the surface temperature
@@ -1527 +1519 @@
             dqsurf(1:ngrid,1:nq) = dqsurf(1:ngrid,1:nq) + zdqssed(1:ngrid,1:nq)
 
-!!            call writediagfi(ngrid,"callsedim_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
-
-            ! ! Test water conservation !AF24: removed
-
          end if ! end of 'sedimentation'
-
-!!         call writediagfi(ngrid,"mass_redist_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
-!!         call writediagfi(ngrid,"mass_redist_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_gas))
 
   ! ---------------
@@ -2096 +2081 @@
          !call wstats(ngrid,"ALB_1st","First Band Surface albedo"," ",2,albedo(:,1))
          call wstats(ngrid,"p","Pressure","Pa",3,pplay)
+         call wstats(ngrid,"emis","Emissivity","",2,emis)
          call wstats(ngrid,"temp","Atmospheric temperature","K",3,zt)
          call wstats(ngrid,"u","Zonal (East-West) wind","m.s-1",3,zu)
@@ -2157 +2143 @@
          endif
 
+       call writediagfi(ngrid,"emis","Emissivity","",2,emis)
        call writediagfi(ngrid,"temp","temperature","K",3,zt)
        call writediagfi(ngrid,"teta","potential temperature","K",3,zh)