Changeset 2618

Timestamp:

Jan 19, 2022, 2:52:13 PM (3 years ago)

Author:

bclmd

Message:

options to fix the temperature or to accelerate the convergence

File:

• trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -137 +137 @@
       real pctsrf_sic(1)
 
-
+      logical fixed_temperature, accelerate_temperature
 
 !     added by AS to avoid the use of adv trac common
@@ -590 +590 @@
       write(*,*)" specOLR = ",specOLR
 
+c option for temperature evolution?
+      write(*,*)"accelerate temperature?"
+      accelerate_temperature=.false.
+      call getin("accelerate_temperature",accelerate_temperature)
+      write(*,*)" accelerate_temperature = ",accelerate_temperature
+
+      write(*,*)"no thermal evolution?"
+      fixed_temperature=.false.
+      call getin("fixed_temperature",fixed_temperature)
+      write(*,*)" fixed_temperature = ",fixed_temperature
+
 c
 c  pour le schema d'ondes de gravite
@@ -1017 +1028 @@
            v(ilayer)=v(ilayer)+dtphys*dv(ilayer)
            temp(ilayer)=temp(ilayer)+dtphys*dtemp(ilayer)
+           if(fixed_temperature) then
+             temp(ilayer)=temp(ilayer)
+           elseif(accelerate_temperature) then
+             temp(ilayer)=temp(ilayer)+
+     &              dtphys*dtemp(ilayer)*max(1.0,play(ilayer)*1e-5)**1.8
+           else
+             temp(ilayer)=temp(ilayer)+dtphys*dtemp(ilayer)
+           endif
         ENDDO