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

Timestamp:

Feb 20, 2024, 4:46:35 PM (9 months ago)

Author:

tbertrand

Message:

Pluto GCM:
Cleaning code and adding extra options for newstart.F

File:

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

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

@@ -45 +45 @@
       use comcstfi_mod, only: pi, g, rcp, r, rad, mugaz, cpp
       use time_phylmdz_mod, only: daysec
-      use callkeys_mod, only: albedo_spectral_mode, calladj, calldifv, &
+      use callkeys_mod, only: fast,albedo_spectral_mode, calladj, calldifv, &
                               callrad, callsoil, nosurf, &
                               aerohaze, corrk, diagdtau,&
@@ -293 +293 @@
       ! For Surface Tracers : (kg/m2/s)
       real dqsurf(ngrid,nq)                 ! Cumulated tendencies.
-      real zdqsurfc(ngrid)                  ! Condense_n2 routine.
+      !real zdqsurfc(ngrid)                  ! Condense_n2 routine.
       REAL zdqsc(ngrid,nq)                  ! Condense_n2 routine.
       real zdqsdif(ngrid,nq)                ! Turbdiff/vdifc routines.
@@ -768 +768 @@
      ! w = F / (rho*area) and rho = P/(r*T)
      ! But first linearly interpolate mass flux to mid-layers
-      do l=1,nlayer-1
+      if (.not.fast) then 
+       do l=1,nlayer-1
          pw(1:ngrid,l)=0.5*(flxw(1:ngrid,l)+flxw(1:ngrid,l+1))
-      enddo
-      pw(1:ngrid,nlayer)=0.5*flxw(1:ngrid,nlayer) ! since flxw(nlayer+1)=0
-      do l=1,nlayer
+       enddo
+       pw(1:ngrid,nlayer)=0.5*flxw(1:ngrid,nlayer) ! since flxw(nlayer+1)=0
+       do l=1,nlayer
          pw(1:ngrid,l)=(pw(1:ngrid,l)*r*pt(1:ngrid,l)) /  &
                        (pplay(1:ngrid,l)*cell_area(1:ngrid))
-      enddo
-      ! omega in Pa/s
-      do l=1,nlayer-1
+       enddo
+       ! omega in Pa/s
+       do l=1,nlayer-1
          omega(1:ngrid,l)=0.5*(flxw(1:ngrid,l)+flxw(1:ngrid,l+1))
-      enddo
-      omega(1:ngrid,nlayer)=0.5*flxw(1:ngrid,nlayer) ! since flxw(nlayer+1)=0
-      do l=1,nlayer
+       enddo
+       omega(1:ngrid,nlayer)=0.5*flxw(1:ngrid,nlayer) ! since flxw(nlayer+1)=0
+       do l=1,nlayer
          omega(1:ngrid,l)=g*omega(1:ngrid,l)/cell_area(1:ngrid)
-      enddo
-
+       enddo
+      endif
 !---------------------------------
 ! II. Compute radiative tendencies
@@ -922 +923 @@
                albedo_equivalent(1:ngrid)=albedo(1:ngrid,1)
                fluxsurfabs_sw(1:ngrid) = fluxtop_dn(1:ngrid)*(1.-albedo(1:ngrid,1))
-               ! TB24:
                fluxabs_sw(1:ngrid)=fluxsurfabs_sw(1:ngrid)
                fluxrad_sky(1:ngrid)    = fluxsurfabs_sw(1:ngrid)
@@ -1138 +1138 @@
 
          pdt(1:ngrid,1:nlayer) = pdt(1:ngrid,1:nlayer)+zdtc(1:ngrid,1:nlayer)
+         pdv(1:ngrid,1:nlayer) = pdv(1:ngrid,1:nlayer)+zdvc(1:ngrid,1:nlayer)
+         pdu(1:ngrid,1:nlayer) = pdu(1:ngrid,1:nlayer)+zduc(1:ngrid,1:nlayer)
          zdtsurf(1:ngrid)      = zdtsurf(1:ngrid) + zdtsurfc(1:ngrid)
 
          pdq(1:ngrid,1:nlayer,1:nq)   = pdq(1:ngrid,1:nlayer,1:nq)+ zdqc(1:ngrid,1:nlayer,1:nq)
-         ! dqsurf(1:ngrid,igcm_n2_ice) = dqsurf(1:ngrid,igcm_n2_ice) + zdqsurfc(1:ngrid)
+         dqsurf(1:ngrid,igcm_n2) = dqsurf(1:ngrid,igcm_n2) + zdqsc(1:ngrid,igcm_n2)
 
 !!         call writediagfi(ngrid,"condense_n2_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
@@ -1324 +1326 @@
 
       tsurf(1:ngrid)=tsurf(1:ngrid)+ptimestep*zdtsurf(1:ngrid)
-
       ! Compute soil temperatures and subsurface heat flux.
       if (callsoil) then
@@ -1378 +1379 @@
       ! Surface pressure.
       ps(1:ngrid) = pplev(1:ngrid,1) + pdpsrf(1:ngrid)*ptimestep
-
 
 
@@ -1645 +1645 @@
       call writediagfi(ngrid,"tsurf","Surface temperature","K",2,tsurf)
       call writediagfi(ngrid,"ps","Surface pressure","Pa",2,ps)
-      call writediagfi(ngrid,"temp","temperature","K",3,zt)
-      call writediagfi(ngrid,"teta","potential temperature","K",3,zh)
-      call writediagfi(ngrid,"u","Zonal wind","m.s-1",3,zu)
-      call writediagfi(ngrid,"v","Meridional wind","m.s-1",3,zv)
-      call writediagfi(ngrid,"w","Vertical wind","m.s-1",3,pw)
-      call writediagfi(ngrid,"p","Pressure","Pa",3,pplay)
+
+      if (.not.fast) then
+       call writediagfi(ngrid,"temp","temperature","K",3,zt)
+       call writediagfi(ngrid,"teta","potential temperature","K",3,zh)
+       call writediagfi(ngrid,"u","Zonal wind","m.s-1",3,zu)
+       call writediagfi(ngrid,"v","Meridional wind","m.s-1",3,zv)
+       call writediagfi(ngrid,"w","Vertical wind","m.s-1",3,pw)
+       call writediagfi(ngrid,"p","Pressure","Pa",3,pplay)
+      endif
 
 !     Subsurface temperatures