source: LMDZ6/branches/Optimisation_LMDZ/libf/phydev/phyaqua_mod.F90 @ 3698

Last change on this file since 3698 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: 794 bytes
Line 
1!
2! $Id: $
3!
4MODULE phyaqua_mod
5
6  IMPLICIT NONE
7
8CONTAINS
9
10  SUBROUTINE iniaqua(nlon, iflag_phys)
11
12  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
13  !  Create an initial state (startphy.nc) for the physics
14  !  Usefull for idealised cases (e.g. aquaplanets or testcases)
15  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
16
17  USE phys_state_var_mod, ONLY: phys_state_var_init
18  USE mod_phys_lmdz_para, ONLY: klon_omp
19  IMPLICIT NONE
20     
21  INTEGER,INTENT(IN) :: nlon,iflag_phys
22
23  CALL phys_state_var_init()
24
25
26  ! Here you could create an initial condition for the physics
27  ! ...
28  ! ... fill in the fields...
29  ! ...
30  ! ... and create a "startphy.nc" file
31      CALL phyredem ("startphy.nc")
32
33  END SUBROUTINE iniaqua
34
35END MODULE phyaqua_mod
Note: See TracBrowser for help on using the repository browser.