source: LMDZ5/trunk/libf/phydev/phyaqua_mod.F90 @ 3998

Last change on this file since 3998 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
RevLine 
[1671]1!
2! $Id: $
3!
[1994]4MODULE phyaqua_mod
[1615]5
[1994]6  IMPLICIT NONE
[1615]7
[1994]8CONTAINS
[1615]9
[2351]10  SUBROUTINE iniaqua(nlon, iflag_phys)
[1994]11
12  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
13  !  Create an initial state (startphy.nc) for the physics
14  !  Usefull for idealised cases (e.g. aquaplanets or testcases)
15  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
16
[2395]17  USE phys_state_var_mod, ONLY: phys_state_var_init
[1994]18  USE mod_phys_lmdz_para, ONLY: klon_omp
19  IMPLICIT NONE
[1671]20     
[1994]21  INTEGER,INTENT(IN) :: nlon,iflag_phys
[1615]22
[1994]23  CALL phys_state_var_init()
[1671]24
25
[1994]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
[1686]31      CALL phyredem ("startphy.nc")
[1671]32
[1994]33  END SUBROUTINE iniaqua
[1615]34
[1994]35END MODULE phyaqua_mod
Note: See TracBrowser for help on using the repository browser.