Changeset 2635 for trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.F90

Timestamp:

Mar 4, 2022, 6:29:05 PM (3 years ago)

Author:

aslmd

Message:

Introduction of cpp_mugaz_mode to decide what source should be used for cpp and mugaz. Force_cpp and check_cpp_match are deprecated (see issue la-communaut-des-mod-les-atmosph-riques-plan-taires/git-trunk#27).

File:

• trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.F90 (modified) (4 diffs)

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

@@ -240 +240 @@
      if (is_master) write(*,*) trim(rname)//": enertest = ",enertest
 
-     if (is_master) write(*,*) trim(rname)//&
-       ": Check to see if cpp values used match gases.def ?"
-     check_cpp_match=.true. ! default value
-     call getin_p("check_cpp_match",check_cpp_match)
-     if (is_master) write(*,*) trim(rname)//&
-       ": check_cpp_match = ",check_cpp_match
-
      if (is_master) write(*,*) trim(rname)//": call radiative transfer ?"
      callrad=.true. ! default value
@@ -1008 +1001 @@
      if (is_master) write(*,*)trim(rname)//": kmixmin = ",kmixmin
 
-     if (is_master) write(*,*)trim(rname)//&
-       ": Does user want to force cpp and mugaz?"
+     if (is_master) write(*,*)'Predefined Cp from dynamics is ',cpp,'J kg^-1 K^-1'
+     if (is_master) write(*,*)'Predefined Mg from dynamics is ',mugaz,'amu'
+
      force_cpp=.false. ! default value
      call getin_p("force_cpp",force_cpp)
-     if (is_master) write(*,*)trim(rname)//": force_cpp = ",force_cpp
-
      if (force_cpp) then
+      if (is_master) write(*,*)trim(rname)//": force_cpp = ",force_cpp
+      if (is_master) write(*,*)trim(rname)//": force_cpp is deprecated.",&
+      "Set cpp_mugaz_mode=1 in callfis to emulate force_cpp=.true."
+      call abort_physic(rname,"Anyway, you need to set force_cpp=.false. to continue.",1)
+     endif
+
+     if (is_master) write(*,*)trim(rname)//&
+     ": where do you want your cpp/mugaz value to come from?",&
+     "=> 0: dynamics (3d), 1: forced in callfis (1d), 2: computed from gases.def (1d)?"
+     cpp_mugaz_mode = 0 ! default value
+     call getin_p("cpp_mugaz_mode",cpp_mugaz_mode)
+     if (is_master) write(*,*)trim(rname)//": cpp_mugaz_mode = ",cpp_mugaz_mode
+
+     if ((cpp_mugaz_mode >= 1).and.(is_master).and.(ngrid>1)) then
+        write(*,*)'    !!!!! Be aware that having different values for cpp and mugaz in the dynamics and physics'
+        write(*,*)'    in 3D can result in very pathological behavior. You have been warned !!!!!'
+     endif
+
+     if (cpp_mugaz_mode == 1) then
        mugaz = -99999.
        if (is_master) write(*,*)trim(rname)//&
@@ -1020 +1031 @@
        call getin_p("mugaz",mugaz)
        IF (mugaz.eq.-99999.) THEN
-         call abort_physic(rname,"mugaz must be set if force_cpp = T",1)
-       ELSE
-         if (is_master) write(*,*)trim(rname)//": mugaz=",mugaz
+         call abort_physic(rname,"mugaz must be set if cpp_mugaz_mode = 1",1)
        ENDIF
        cpp = -99999.
@@ -1029 +1038 @@
        call getin_p("cpp",cpp)
        IF (cpp.eq.-99999.) THEN
-           PRINT *, "cpp must be set if force_cpp = T"
+           PRINT *, "cpp must be set if cpp_mugaz_mode = 1"
            STOP
-       ELSE
-           if (is_master) write(*,*)trim(rname)//": cpp=",cpp
        ENDIF
-     endif ! of if (force_cpp)
+       if (is_master) write(*,*)'New Cp from callfis is ',cpp,'J kg^-1 K^-1'
+       if (is_master) write(*,*)'New Mg from callfis is ',mugaz,'amu'
+  
+     endif ! of if (cpp_mugaz_mode == 1)
      call su_gases
      call calc_cpp_mugaz