Changeset 3116 for trunk/LMDZ.GENERIC


Ignore:
Timestamp:
Nov 8, 2023, 2:33:52 PM (14 months ago)
Author:
mmaurice
Message:

GENERIC PCM:
Fix newstart (longitude and latitude were not correctly written in restartfi.nc)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/newstart.F

    r3100 r3116  
    4646      use phyetat0_mod, only: phyetat0
    4747      use exner_hyb_m, only: exner_hyb
     48      use geometry_mod, only: longitude,  ! longitudes (rad)
     49     &                         latitude,  ! latitudes (rad)                       
     50     &                         cell_area ! physics grid area (m2)                       
    4851      implicit none
    4952
     
    112115      real,ALLOCATABLE :: ith(:,:,:),ithfi(:,:) ! thermal inertia (3D)
    113116      real surfith(iip1,jjp1),surfithfi(ngridmx) ! surface thermal inertia (2D)
    114       REAL latfi(ngridmx),lonfi(ngridmx),airefi(ngridmx)
     117!      REAL latfi(ngridmx),lonfi(ngridmx),airefi(ngridmx)
    115118
    116119      INTEGER i,j,l,isoil,ig,idum
     
    13751378
    13761379          do ig=1, ngridmx
    1377              teque= 335.0-60.0*sin(latfi(ig))*sin(latfi(ig))-
    1378      &            10.0*cos(latfi(ig))*cos(latfi(ig))
     1380             teque= 335.0-60.0*sin(latitude(ig))*sin(latitude(ig))-
     1381     &            10.0*cos(latitude(ig))*cos(latitude(ig))
    13791382             tsurf(ig) = MAX(220.0,teque)
    13801383          end do
     
    16531656C
    16541657
    1655       call physdem0("restartfi.nc",lonfi,latfi,nsoilmx,ngridmx,llm,
     1658      call physdem0("restartfi.nc",longitude,latitude,nsoilmx,ngridmx,
     1659     &              llm,
    16561660     &              nqtot,dtphys,real(day_ini),0.0,
    1657      &              airefi,albfi,ithfi,zmea,zstd,zsig,zgam,zthe)
     1661     &              cell_area,albfi,ithfi,zmea,zstd,zsig,zgam,zthe)
    16581662      call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
    16591663     &                dtphys,real(day_ini),
Note: See TracChangeset for help on using the changeset viewer.