Changeset 2621 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Jan 19, 2022, 6:12:03 PM (3 years ago)

Author:

bclmd

Message:

add aceleration factor and remove the option for fixed temperature (not fully correct)

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 BC to accelerate convergence
+      logical accelerate_temperature
+      real factor_acceleration
 
 !     added by AS to avoid the use of adv trac common
@@ -596 +598 @@
       write(*,*)" accelerate_temperature = ",accelerate_temperature
 
-      write(*,*)"no thermal evolution?"
-      fixed_temperature=.false.
-      call getin("fixed_temperature",fixed_temperature)
-      write(*,*)" fixed_temperature = ",fixed_temperature
+      write(*,*)"factor for temperature acceleration"
+      factor_acceleration=1.0
+      call getin("factor_acceleration",factor_acceleration)
+      write(*,*)" factor_acceleration = ",factor_acceleration
+
+
 
 c
@@ -1027 +1031 @@
            u(ilayer)=u(ilayer)+dtphys*du(ilayer)
            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
+           if(accelerate_temperature) then
+             temp(ilayer)=temp(ilayer)+dtphys*dtemp(ilayer) *
+     &         max(1.0,play(ilayer)*1e-5)**factor_acceleration
            else
              temp(ilayer)=temp(ilayer)+dtphys*dtemp(ilayer)