Changeset 1984 for LMDZ5/trunk/libf/dyn3d_common

Timestamp:

Feb 18, 2014, 10:59:29 AM (10 years ago)

Author:

Ehouarn Millour

Message:

Cleanup and update of "grilles_gcm_netcdf_sub": removed old version, removed flag "grilles_gcm_netcdf" used to trigger output of grilles_gcm.nc file (file will now always be created) and fixed writing vertical levels presnivs in the file.
JBM+EM

Location:

LMDZ5/trunk/libf/dyn3d_common

Files:

• grilles_gcm_netcdf.F (deleted)
• grilles_gcm_netcdf_sub.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.F90

@@ -1 +1 @@
 !
-! $Header$
+! $Id: $
 !
 ! This subroutine creates the file grilles_gcm.nc containg longitudes and
 ! latitudes in degrees for grid u and v. This subroutine is called from
-! ce0l if grilles_gcm_netcdf=TRUE. This subroutine corresponds to the first 
+! ce0l. This subroutine corresponds to the first 
 ! part in the program create_fausse_var.
 !
@@ -48 +48 @@
   INTEGER :: area_id
   INTEGER :: mask_id
+  INTEGER :: presnivs_id
   
   rad = 6400000
@@ -80 +81 @@
   ! CREATION OUTPUT
   ! ouverture fichier netcdf de sortie out
-  status=NF_CREATE('grilles_gcm.nc',NF_NOCLOBBER,ncid_out)
+  status=NF_CREATE('grilles_gcm.nc',NF_CLOBBER,ncid_out)
+  CALL handle_err(status)
   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)
 
 
@@ -136 +142 @@
 
   status=NF_ENDDEF(ncid_out)
+  write(*,*) "COUCOU 6"
+  CALL handle_err(status)
   ! 5) ----- FERMETURE DES FICHIERS NETCDF------------------
   ! --------------------------------------------------------
@@ -143 +151 @@
 #ifdef NC_DOUBLE
   status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonuid,1,iim+1,rlonudeg)
+  CALL handle_err(status)
   status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonvid,1,iim+1,rlonvdeg)
+  CALL handle_err(status)
   status=NF_PUT_VARA_DOUBLE(ncid_out,out_latuid,1,jjm+1,rlatudeg)
+  CALL handle_err(status)
   status=NF_PUT_VARA_DOUBLE(ncid_out,out_latvid,1,jjm,rlatvdeg)
+  CALL handle_err(status)
 #else
   status=NF_PUT_VARA_REAL(ncid_out,out_lonuid,1,iim+1,rlonudeg)
+  CALL handle_err(status)
   status=NF_PUT_VARA_REAL(ncid_out,out_lonvid,1,iim+1,rlonvdeg)
+  CALL handle_err(status)
   status=NF_PUT_VARA_REAL(ncid_out,out_latuid,1,jjm+1,rlatudeg)
+  CALL handle_err(status)
   status=NF_PUT_VARA_REAL(ncid_out,out_latvid,1,jjm,rlatvdeg)
+  CALL handle_err(status)
 #endif
 
@@ -171 +187 @@
 #ifdef NC_DOUBLE
   status=NF_PUT_VARA_DOUBLE(ncid_out,out_varid,start, count,temp)
+  CALL handle_err(status)
 #else
   status=NF_PUT_VARA_REAL(ncid_out,out_varid,start, count,temp)
+  CALL handle_err(status)
 #endif
 
   ! On re-ouvre le fichier pour rajouter 4 nouvelles variables necessaire pour INCA
 ! lev - phis - aire - mask
-  rlevdeg(:) = presnivs
+! rlevdeg(:) = presnivs
+  rlevdeg(:) = presnivs(:)
   phis_loc(:,:) = phis(:,:)/g
 
 ! niveaux de pression verticaux
   status = NF_REDEF (ncid_out)
+  CALL handle_err(status)
   status=NF_DEF_DIM(ncid_out,'lev',llm,out_levdim)
+  CALL handle_err(status)
+  status=NF_DEF_VAR(ncid_out,'presnivs',NF_FLOAT,1,out_levdim,&
+                    presnivs_id)
+  CALL handle_err(status)
   
 ! fields
@@ -196 +220 @@
 
   status=NF_ENDDEF(ncid_out)
+  CALL handle_err(status)
 
   ! ecriture des variables
 #ifdef NC_DOUBLE
-  status=NF_PUT_VARA_DOUBLE(ncid_out,out_levid,1,llm,rlevdeg)
+  status=NF_PUT_VARA_DOUBLE(ncid_out,presnivs_id,1,llm,rlevdeg)
+  CALL handle_err(status)
 #else
   status=NF_PUT_VARA_REAL(ncid_out,out_levid,1,llm,rlevdeg)
+  CALL handle_err(status)
 #endif
 
@@ -214 +241 @@
 
   status = nf_put_vara_double(ncid_out, phis_id,start,count, phis_loc)
+  CALL handle_err(status)
   status = nf_put_vara_double(ncid_out, area_id,start,count, aire)
+  CALL handle_err(status)
   masque_int(:,:) = nINT(masque(:,:))
   status = nf_put_vara_int(ncid_out, mask_id,start,count,masque_int)