Changeset 5270 for LMDZ6/trunk/libf/dyn3d_common

Timestamp:

Oct 24, 2024, 1:55:38 PM (5 weeks ago)

Author:

abarral

Message:

Replace F77 netcdf library by F90 netcdf library

File:

• LMDZ6/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.f90 (modified) (10 diffs)

LMDZ6/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.f90

@@ -14 +14 @@
   USE comconst_mod, ONLY: cpp, kappa, g, omeg, daysec, rad, pi
   USE comvert_mod, ONLY: presnivs, preff, pa
-  use netcdf, only: nf90_def_var, nf90_int, nf90_float, nf90_put_var
+  USE netcdf, ONLY: nf90_def_var, nf90_int, nf90_float, nf90_put_var, nf90_clobber, nf90_64bit_offset, nf90_def_dim, &
+          nf90_put_att, nf90_enddef, nf90_create
   
   IMPLICIT NONE
@@ -21 +22 @@
   INCLUDE "paramet.h"
   INCLUDE "comgeom.h"
-  INCLUDE "netcdf.inc"
 
 !======================== 
@@ -102 +102 @@
   ! OPEN output netcdf file
   !-------------------------
-  status=NF_CREATE('grilles_gcm.nc',IOR(NF_CLOBBER,NF_64BIT_OFFSET),ncid_out)
+  status=nf90_create('grilles_gcm.nc',IOR(nf90_clobber,nf90_64bit_offset),ncid_out)
   CALL handle_err(status)
   
   ! DEFINE output dimensions
   !-------------------------
-  status=NF_DEF_DIM(ncid_out,'lonu',iim+1,out_lonudim)
-  CALL handle_err(status)
-  status=NF_DEF_DIM(ncid_out,'lonv',iim+1,out_lonvdim)
-  CALL handle_err(status)
-  status=NF_DEF_DIM(ncid_out,'latu',jjm+1,out_latudim)
-  CALL handle_err(status)
-  status=NF_DEF_DIM(ncid_out,'latv',jjm,out_latvdim)
-  CALL handle_err(status)
-  !
-  status=NF_DEF_DIM(ncid_out,'lev',llm,out_levdim)
+  status=nf90_def_dim(ncid_out,'lonu',iim+1,out_lonudim)
+  CALL handle_err(status)
+  status=nf90_def_dim(ncid_out,'lonv',iim+1,out_lonvdim)
+  CALL handle_err(status)
+  status=nf90_def_dim(ncid_out,'latu',jjm+1,out_latudim)
+  CALL handle_err(status)
+  status=nf90_def_dim(ncid_out,'latv',jjm,out_latvdim)
+  CALL handle_err(status)
+  !
+  status=nf90_def_dim(ncid_out,'lev',llm,out_levdim)
   CALL handle_err(status)
   
@@ -124 +124 @@
   status=NF90_DEF_VAR(ncid_out,'lonu',NF90_FLOAT,out_lonudim, out_lonuid)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,out_lonuid,'units', 12,'degrees_east')
-  status=NF_PUT_ATT_TEXT(ncid_out,out_lonuid,'long_name',19,'Longitude on u grid')
+  status=nf90_put_att(ncid_out,out_lonuid,'units','degrees_east')
+  status=nf90_put_att(ncid_out,out_lonuid,'long_name','Longitude on u grid')
   !   Longitudes on "v" dynamical grid
   status=NF90_DEF_VAR(ncid_out,'lonv',NF90_FLOAT,out_lonvdim, out_lonvid)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,out_lonvid,'units', 12,'degrees_east')
-  status=NF_PUT_ATT_TEXT(ncid_out,out_lonvid,'long_name', 19,'Longitude on v grid')
+  status=nf90_put_att(ncid_out,out_lonvid,'units','degrees_east')
+  status=nf90_put_att(ncid_out,out_lonvid,'long_name','Longitude on v grid')
   !   Latitudes on "u" dynamical grid
   status=NF90_DEF_VAR(ncid_out,'latu',NF90_FLOAT,out_latudim, out_latuid)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,out_latuid,'units', 13,'degrees_north')
-  status=NF_PUT_ATT_TEXT(ncid_out,out_latuid,'long_name', 18,'Latitude on u grid')
+  status=nf90_put_att(ncid_out,out_latuid,'units','degrees_north')
+  status=nf90_put_att(ncid_out,out_latuid,'long_name','Latitude on u grid')
   !  Latitudes on "v" dynamical grid
   status=NF90_DEF_VAR(ncid_out,'latv',NF90_FLOAT,out_latvdim, out_latvid)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,out_latvid,'units', 13,'degrees_north')
-  status=NF_PUT_ATT_TEXT(ncid_out,out_latvid,'long_name', 18,'Latitude on v grid')
+  status=nf90_put_att(ncid_out,out_latvid,'units','degrees_north')
+  status=nf90_put_att(ncid_out,out_latvid,'long_name','Latitude on v grid')
   !  "u" lat/lon dynamical grid
   out_dim(1)=out_lonudim
@@ -146 +146 @@
   status=NF90_DEF_VAR(ncid_out,'grille_u',NF90_FLOAT,out_dim, out_uid)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,out_uid,'units', 3,'m/s')
-  status=NF_PUT_ATT_TEXT(ncid_out,out_uid,'long_name', 21,'u-wind dynamical grid')
+  status=nf90_put_att(ncid_out,out_uid,'units','m/s')
+  status=nf90_put_att(ncid_out,out_uid,'long_name','u-wind dynamical grid')
   !  "v" lat/lon dynamical grid
   out_dim(1)=out_lonvdim
@@ -153 +153 @@
   status=NF90_DEF_VAR(ncid_out,'grille_v',NF90_FLOAT,out_dim, out_vid)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,out_vid,'units', 3,'m/s')
-  status=NF_PUT_ATT_TEXT(ncid_out,out_vid,'long_name', 21,'v-wind dynamical grid')
+  status=nf90_put_att(ncid_out,out_vid,'units','m/s')
+  status=nf90_put_att(ncid_out,out_vid,'long_name','v-wind dynamical grid')
   !  "s" (scalar) lat/lon dynamical grid
   out_dim(1)=out_lonvdim
@@ -160 +160 @@
   status=NF90_DEF_VAR(ncid_out,'grille_s',NF90_FLOAT,out_dim, out_tempid)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,out_tempid,'units', 6,'Kelvin')
-  status=NF_PUT_ATT_TEXT(ncid_out,out_tempid,'long_name',21,'scalar dynamical grid')
+  status=nf90_put_att(ncid_out,out_tempid,'units','Kelvin')
+  status=nf90_put_att(ncid_out,out_tempid,'long_name','scalar dynamical grid')
   !
   ! for INCA :
@@ -167 +167 @@
   status=NF90_DEF_VAR(ncid_out,'presnivs',NF90_FLOAT,out_levdim, presnivs_id)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,presnivs_id,'units', 2,'Pa')
-  status=NF_PUT_ATT_TEXT(ncid_out,presnivs_id,'long_name',15,'Vertical levels')
+  status=nf90_put_att(ncid_out,presnivs_id,'units','Pa')
+  status=nf90_put_att(ncid_out,presnivs_id,'long_name','Vertical levels')
   ! surface geopotential height: named "phis" as the sfc geopotential, is actually phis/g
   out_dim(1)=out_lonvdim
@@ -174 +174 @@
   status = nf90_def_var(ncid_out,'phis',NF90_FLOAT,out_dim,phis_id)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,phis_id,'units', 1,'m')
-  status=NF_PUT_ATT_TEXT(ncid_out,phis_id,'long_name',27,'surface geopotential height')
+  status=nf90_put_att(ncid_out,phis_id,'units','m')
+  status=nf90_put_att(ncid_out,phis_id,'long_name','surface geopotential height')
   ! gridcell area
   status = nf90_def_var(ncid_out,'aire',NF90_FLOAT,out_dim,area_id)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,area_id,'units', 2,'m2')
-  status=NF_PUT_ATT_TEXT(ncid_out,area_id,'long_name',13,'gridcell area')
+  status=nf90_put_att(ncid_out,area_id,'units','m2')
+  status=nf90_put_att(ncid_out,area_id,'long_name','gridcell area')
   ! land-sea mask (nearest integer approx)
   status = nf90_def_var(ncid_out,'mask',NF90_INT,out_dim,mask_id)
   CALL handle_err(status)
-  status=NF_PUT_ATT_TEXT(ncid_out,mask_id,'long_name',27,'land-sea mask (nINT approx)')
+  status=nf90_put_att(ncid_out,mask_id,'long_name','land-sea mask (nINT approx)')
   
   ! END the 'define' mode in netCDF file
-  status=NF_ENDDEF(ncid_out)  
+  status=nf90_enddef(ncid_out)
   CALL handle_err(status)
   
@@ -232 +232 @@
 
 SUBROUTINE handle_err(status)
-  INCLUDE "netcdf.inc"
+  USE netcdf, ONLY: nf90_noerr, nf90_strerror
+  IMPLICIT NONE
 
   INTEGER status
-  IF (status.NE.nf_noerr) THEN
-     PRINT *,NF_STRERROR(status)
+  IF (status.NE.nf90_noerr) THEN
+     PRINT *,nf90_strerror(status)
      CALL abort_gcm('grilles_gcm_netcdf','netcdf error',1)
   ENDIF