Changeset 3116 for trunk/LMDZ.GENERIC

Timestamp:

Nov 8, 2023, 2:33:52 PM (13 months ago)

Author:

mmaurice

Message:

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

File:

• trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/newstart.F (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/newstart.F

@@ -46 +46 @@
       use phyetat0_mod, only: phyetat0
       use exner_hyb_m, only: exner_hyb
+      use geometry_mod, only: longitude,  ! longitudes (rad)
+     &                         latitude,  ! latitudes (rad)                       
+     &                         cell_area ! physics grid area (m2)                        
       implicit none
 
@@ -112 +115 @@
       real,ALLOCATABLE :: ith(:,:,:),ithfi(:,:) ! thermal inertia (3D)
       real surfith(iip1,jjp1),surfithfi(ngridmx) ! surface thermal inertia (2D)
-      REAL latfi(ngridmx),lonfi(ngridmx),airefi(ngridmx)
+!      REAL latfi(ngridmx),lonfi(ngridmx),airefi(ngridmx)
 
       INTEGER i,j,l,isoil,ig,idum
@@ -1375 +1378 @@
 
           do ig=1, ngridmx
-             teque= 335.0-60.0*sin(latfi(ig))*sin(latfi(ig))-
-     &            10.0*cos(latfi(ig))*cos(latfi(ig))
+             teque= 335.0-60.0*sin(latitude(ig))*sin(latitude(ig))-
+     &            10.0*cos(latitude(ig))*cos(latitude(ig))
              tsurf(ig) = MAX(220.0,teque)
           end do
@@ -1653 +1656 @@
 C
 
-      call physdem0("restartfi.nc",lonfi,latfi,nsoilmx,ngridmx,llm,
+      call physdem0("restartfi.nc",longitude,latitude,nsoilmx,ngridmx,
+     &              llm,
      &              nqtot,dtphys,real(day_ini),0.0,
-     &              airefi,albfi,ithfi,zmea,zstd,zsig,zgam,zthe)
+     &              cell_area,albfi,ithfi,zmea,zstd,zsig,zgam,zthe)
       call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
      &                dtphys,real(day_ini),