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

Last change on this file since 3605 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
Line 
1!
2! $Id: $
3!
4SUBROUTINE phyredem (fichnom)
5
6  USE geometry_mod, ONLY : longitude_deg, latitude_deg
7  USE iostart, ONLY: open_restartphy, close_restartphy, put_var, put_field
8
9  IMPLICIT NONE
10
11  CHARACTER(LEN=*),INTENT(IN) :: fichnom
12
13  INTEGER,PARAMETER :: tab_cntrl_len=100
14  REAL :: tab_cntrl(tab_cntrl_len)
15
16  ! open file
17
18  CALL open_restartphy(fichnom)
19
20  ! tab_cntrl() contains run parameters
21
22  tab_cntrl(:)=0.0
23 
24
25  CALL put_var("controle", "Control parameters", tab_cntrl)
26
27  ! coordinates
28
29  CALL put_field("longitude", "Longitudes on physics grid", longitude_deg)
30     
31  CALL put_field("latitude", "Latitudes on physics grid", latitude_deg)
32
33  ! close file
34
35  CALL close_restartphy
36  !$OMP BARRIER
37
38END SUBROUTINE phyredem
Note: See TracBrowser for help on using the repository browser.