Changeset 1511 for LMDZ5/trunk/libf/dyn3d/grilles_gcm_netcdf_sub.F90

Timestamp:

Apr 28, 2011, 5:21:47 PM (14 years ago)

Author:

jghattas

Message:

Added 4 more variables in output file grilles_gcm.nc created by ce0l :
lev, phis, aire and mask. These variables are needed by INCA.

A Cozic

File:

• LMDZ5/trunk/libf/dyn3d/grilles_gcm_netcdf_sub.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/dyn3d/grilles_gcm_netcdf_sub.F90

@@ -7 +7 @@
 ! part in the program create_fausse_var.
 !
-SUBROUTINE grilles_gcm_netcdf_sub
+SUBROUTINE grilles_gcm_netcdf_sub(masque,phis)
 
   IMPLICIT NONE
@@ -19 +19 @@
   INCLUDE "serre.h"
 
+
+  REAL,DIMENSION(iip1,jjp1),INTENT(IN)  :: masque ! masque terre/mer
+  REAL,DIMENSION(iip1,jjp1),INTENT(IN)  :: phis   ! geopotentiel au sol
+
   REAL temp(iim+1,jjm+1)
   ! Attributs netcdf sortie
   INTEGER ncid_out,rcode_out
-  INTEGER out_lonuid,out_lonvid,out_latuid,out_latvid
+  INTEGER out_lonuid,out_lonvid,out_latuid,out_latvid,out_levid
   INTEGER out_varid
   INTEGER out_lonudim,out_lonvdim
   INTEGER out_latudim,out_latvdim,out_dim(3)
+  INTEGER out_levdim
 
   INTEGER, PARAMETER :: longcles = 20
@@ -33 +38 @@
 
   INTEGER status,i,j
-  REAL rlatudeg(jjp1),rlatvdeg(jjm)
+  REAL rlatudeg(jjp1),rlatvdeg(jjm),rlevdeg(llm)
   REAL rlonudeg(iip1),rlonvdeg(iip1)
 
   REAL dlon1(iip1),dlon2(iip1),dlat1(jjp1),dlat2(jjp1)
   REAL acoslat,dxkm,dykm,resol(iip1,jjp1)
-
+  REAL,DIMENSION(iip1,jjp1)  :: phis_loc
+  INTEGER masque_int(iip1,jjp1)
+  INTEGER :: phis_id
+  INTEGER :: area_id
+  INTEGER :: mask_id
+  
   rad = 6400000
   omeg = 7.272205e-05
@@ -165 +175 @@
 #endif
 
-
+  ! On re-ouvre le fichier pour rajouter 4 nouvelles variables necessaire pour INCA
+! lev - phis - aire - mask
+  rlevdeg(:) = presnivs
+  phis_loc(:,:) = phis(:,:)/g
+
+! niveaux de pression verticaux
+  status = NF_REDEF (ncid_out)
+  status=NF_DEF_DIM(ncid_out,'lev',llm,out_levdim)
+  
+! fields
+  out_dim(1)=out_lonvdim
+  out_dim(2)=out_latudim
+
+  status = nf_def_var(ncid_out,'phis',NF_FLOAT,2,out_dim,phis_id)
+  CALL handle_err(status)
+  status = nf_def_var(ncid_out,'aire',NF_FLOAT,2,out_dim,area_id)
+  CALL handle_err(status)
+  status = nf_def_var(ncid_out,'mask',NF_INT  ,2,out_dim,mask_id)
+  CALL handle_err(status)
+
+  status=NF_ENDDEF(ncid_out)
+
+  ! ecriture des variables
+#ifdef NC_DOUBLE
+  status=NF_PUT_VARA_DOUBLE(ncid_out,out_levid,1,llm,rlevdeg)
+#else
+  status=NF_PUT_VARA_REAL(ncid_out,out_levid,1,llm,rlevdeg)
+#endif
+
+  start(1)=1
+  start(2)=1
+  start(3)=1
+  start(4)=0
+  COUNT(1)=iip1
+  COUNT(2)=jjp1
+  COUNT(3)=1
+  COUNT(4)=0
+
+  status = nf_put_vara_double(ncid_out, phis_id,start,count, phis_loc)
+  status = nf_put_vara_double(ncid_out, area_id,start,count, aire)
+  masque_int(:,:) = nINT(masque(:,:))
+  status = nf_put_vara_int(ncid_out, mask_id,start,count,masque_int)
+  CALL handle_err(status)
+  
   ! fermeture du fichier netcdf
   CALL ncclos(ncid_out,rcode_out)