source: LMDZ5/trunk/libf/phydev/phyredem.F90 @ 3719

Last change on this file since 3719 was 2395, checked in by Ehouarn Millour, 9 years ago

Bug fix for aquaplanets: "rlat" and "rlon" were uninitialized when written to startphy.nc; phyredem should write "longitude_deg" and "latitude_deg" to the file. Overall "rlat" and "rlon" should not be used and "latitude_deg" and "longitude_deg" should be used instead in the physics. To be further cleaned up.
Some cleanup on that mater also made in phydev.
Note that this change will make bench test results different for longitudes and latitudes, because of roundoff effects.
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 735 bytes
RevLine 
[1671]1!
2! $Id: $
3!
[1994]4SUBROUTINE phyredem (fichnom)
[1671]5
[2395]6  USE geometry_mod, ONLY : longitude_deg, latitude_deg
[1994]7  USE iostart, ONLY: open_restartphy, close_restartphy, put_var, put_field
[1671]8
[1994]9  IMPLICIT NONE
[1671]10
[1994]11  CHARACTER(LEN=*),INTENT(IN) :: fichnom
[1671]12
[1994]13  INTEGER,PARAMETER :: tab_cntrl_len=100
14  REAL :: tab_cntrl(tab_cntrl_len)
[1671]15
[1994]16  ! open file
[1671]17
[1994]18  CALL open_restartphy(fichnom)
[1671]19
[1994]20  ! tab_cntrl() contains run parameters
[1671]21
[1994]22  tab_cntrl(:)=0.0
[1671]23 
24
[1994]25  CALL put_var("controle", "Control parameters", tab_cntrl)
[1671]26
[1994]27  ! coordinates
[1671]28
[2395]29  CALL put_field("longitude", "Longitudes on physics grid", longitude_deg)
[1671]30     
[2395]31  CALL put_field("latitude", "Latitudes on physics grid", latitude_deg)
[1671]32
[1994]33  ! close file
[1671]34
[1994]35  CALL close_restartphy
36  !$OMP BARRIER
[1671]37
[1994]38END SUBROUTINE phyredem
Note: See TracBrowser for help on using the repository browser.