Changeset 2422 for LMDZ5/trunk/libf/phylmd/time_phylmdz_mod.F90

Timestamp:

Jan 6, 2016, 12:37:41 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Small modification in the way time and calendar are handled: Now all the time keeping is done in the physics and only the timestep is transfered from the dynamics to the physics. Due to changes in computations and roundoffs this will change reference bench results.
The implementation of this change in phymar is left as future work.
EM

File:

• LMDZ5/trunk/libf/phylmd/time_phylmdz_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/time_phylmdz_mod.F90

@@ -77 +77 @@
   END SUBROUTINE init_iteration
 
+  SUBROUTINE update_time(pdtphys_)
+  ! This subroutine updates the module saved variables.
+  USE ioipsl, ONLY : ymds2ju
+  USE phys_cal_mod, ONLY: phys_cal_update
+  USE print_control_mod, ONLY: lunout
+  IMPLICIT NONE
+    REAL,INTENT(IN) :: pdtphys_
+    REAL            :: julian_date
+    
+    ! Check if the physics timestep has changed
+    IF ( ABS( (pdtphys-pdtphys_) / ((pdtphys+pdtphys_)/2))> 10.*EPSILON(pdtphys_)) THEN
+       WRITE(lunout,*) "WARNING ! Physics time step changes from a call to the next",pdtphys_,pdtphys
+       WRITE(lunout,*) "Not sure the physics parametrizations can handle this..."
+    ENDIF
+    pdtphys=pdtphys_
+    
+    ! Update elapsed time since begining of run:
+    current_time=current_time+pdtphys
+
+    ! Compute corresponding Julian date and update calendar
+    CALL ymds2ju(annee_ref,1,day_ini,start_time+current_time,julian_date)
+    CALL phys_cal_update(julian_date)
+    
+  END SUBROUTINE update_time
+
 END MODULE time_phylmdz_mod