Changeset 3422 for trunk/LMDZ.PLUTO

Timestamp:

Aug 26, 2024, 6:11:23 PM (3 months ago)

Author:

afalco

Message:

Pluto PCM: no gcm works in parallel. example restart.def for old-to-new pluto start files.
AF

Location:

trunk/LMDZ.PLUTO

Files:

• deftank/README.md (added)
• deftank/restart.old_pluto_to_new.def (added)
• libf/dynphy_lonlat/phypluto/leapfrog_nogcm.F (modified) (2 diffs)
• libf/phypluto/callcorrk_pluto_mod.F90 (modified) (1 diff)
• libf/phypluto/inifis_mod.F90 (modified) (2 diffs)
• libf/phypluto/physiq_mod.F90 (modified) (5 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/leapfrog_nogcm.F

@@ -13 +13 @@
 #endif
       USE infotrac, ONLY: nqtot,ok_iso_verif,tname
-      USE guide_mod, ONLY : guide_main
       USE write_field, ONLY: writefield
       USE callkeys_mod, only: tau_n2, tau_ch4, tau_co
@@ -24 +23 @@
       use exner_milieu_m, only: exner_milieu
       use cpdet_mod, only: cpdet,tpot2t,t2tpot
-      use sponge_mod, only: callsponge,mode_sponge,sponge
        use comuforc_h
       USE comvert_mod, ONLY: ap,bp,pressure_exner,presnivs,

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

@@ -621 +621 @@
             print*,'transfer kmatrix bounds, exiting.'
             print*,'t(',k,')=',tlevrad(k),' < ',tgasmin
-            call abort
+            ! call abort
           elseif(tlevrad(k).gt.tgasmax)then
             print*,'Maximum temperature is outside the radiative'

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

@@ -1387 +1387 @@
      if (is_master) write(*,*)trim(rname)//": generic_condensation = ",generic_condensation
 
-! Test of incompatibility:
      if (is_master) write(*,*)trim(rname)//": Spectral Dependant albedo ?"
      albedo_spectral_mode=.false. ! default value
@@ -1436 +1435 @@
      call getin_p("cpp_mugaz_mode",cpp_mugaz_mode)
      if (is_master) write(*,*)trim(rname)//": cpp_mugaz_mode = ",cpp_mugaz_mode
+
+
+! Test of incompatibility:
+     if (condcosurf.and.no_n2frost) then
+       call abort_physic(rname, "CO surface condensation and no_n2frost are both active which may not be relevant", 1)
+     end if
 
      if ((cpp_mugaz_mode >= 1).and.(is_master).and.(ngrid>1)) then

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

@@ -2287 +2287 @@
 
 !       ! Oceanic layers !AF24: removed
-
       ! ! Thermal plume model  !AF24: removed
-
 !        GW non-oro outputs  !AF24: removed
 
@@ -2309 +2307 @@
 !           call writediagfi(ngrid,"fluxsurflwcs","lw back radiation (cs).","W m-2",2,fluxsurf_lw1)
 
-         ! if(ok_slab_ocean) then
-         !    call writediagfi(ngrid,"GND","heat flux from ground","W m-2",2,fluxgrdocean)
-         ! else
          call writediagfi(ngrid,"GND","heat flux from ground","W m-2",2,fluxgrd)
          ! endif
@@ -2366 +2361 @@
         call writediagfi(ngrid,"zdtlw","LW heating","T s-1",3,zdtlw)
         call writediagfi(ngrid,"dtrad","radiative heating","K s-1",3,dtrad)
-      ! call writediagfi(ngrid,"zdtdyn","Dyn. heating","T s-1",3,zdtdyn)
+        call writediagfi(ngrid,"zdtdyn","Dyn. heating","T s-1",3,zdtdyn)
 
         !Pluto specific
@@ -2396 +2391 @@
          do iq=1,nq
           call writediagfi(ngrid,noms(iq),noms(iq),'kg/kg',3,zq(1,1,iq))
-          !  call writediagfi(ngrid,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
-          !                   'kg m^-2',2,qsurf_hist(1,iq) )
+         !   call writediagfi(ngrid,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
+         !                    'kg m^-2',2,qsurf_hist(1,iq) )
            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) )
+         call writediagfi(ngrid,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
+                         'kg m^-2',2,qsurf(1,iq) )
 
             ! if(watercond.or.CLFvarying)then  !AF24: removed
@@ -2419 +2414 @@
 
          !Pluto specific
-         !call writediagfi(ngrid,'n2_iceflux','n2_iceflux',"kg m^-2 s^-1",2,flusurf(1,igcm_n2) )
-         ! call writediagfi(ngrid,'haze_reff','haze_reff','m',3,reffrad(1,1,1))
+         call writediagfi(ngrid,'n2_iceflux','n2_iceflux',"kg m^-2 s^-1",2,flusurf(1,igcm_n2) )
+         call writediagfi(ngrid,'haze_reff','haze_reff','m',3,reffrad(1,1,1))
          if (methane) then
             call writediagfi(ngrid,'ch4_iceflux','ch4_iceflux',&