Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/ce0l.F90 =================================================================== --- LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/ce0l.F90 (revision 1487) +++ LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/ce0l.F90 (revision 1488) @@ -91,4 +91,12 @@ END IF + IF (grilles_gcm_netcdf) THEN + WRITE(lunout,'(//)') + WRITE(lunout,*) ' *************************** ' + WRITE(lunout,*) ' *** grilles_gcm_netcdf *** ' + WRITE(lunout,*) ' *************************** ' + WRITE(lunout,'(//)') + CALL grilles_gcm_netcdf_sub() + END IF #endif ! of #ifndef CPP_EARTH #else Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/conf_gcm.F =================================================================== --- LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/conf_gcm.F (revision 1487) +++ LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/conf_gcm.F (revision 1488) @@ -841,4 +841,10 @@ ok_etat0 = .TRUE. CALL getin('ok_etat0',ok_etat0) + +!Config Key = grilles_gcm_netcdf +!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit +!Config Def = n + grilles_gcm_netcdf = .FALSE. + CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf) write(lunout,*)' #########################################' @@ -887,4 +893,5 @@ write(lunout,*)' ok_limit = ', ok_limit write(lunout,*)' ok_etat0 = ', ok_etat0 + write(lunout,*)' grilles_gcm_netcdf = ', grilles_gcm_netcdf c RETURN Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/grilles_gcm_netcdf_sub.F90 =================================================================== --- LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/grilles_gcm_netcdf_sub.F90 (revision 1488) +++ LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/grilles_gcm_netcdf_sub.F90 (revision 1488) @@ -0,0 +1,184 @@ +! +! $Header$ +! +! 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 +! part in the program create_fausse_var. +! +SUBROUTINE grilles_gcm_netcdf_sub + + IMPLICIT NONE + + INCLUDE "dimensions.h" + INCLUDE "paramet.h" + INCLUDE "comconst.h" + INCLUDE "comgeom.h" + INCLUDE "comvert.h" + INCLUDE "netcdf.inc" + INCLUDE "serre.h" + + 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_varid + INTEGER out_lonudim,out_lonvdim + INTEGER out_latudim,out_latvdim,out_dim(3) + + INTEGER, PARAMETER :: longcles = 20 + REAL clesphy0(longcles) + + INTEGER start(4),COUNT(4) + + INTEGER status,i,j + REAL rlatudeg(jjp1),rlatvdeg(jjm) + REAL rlonudeg(iip1),rlonvdeg(iip1) + + REAL dlon1(iip1),dlon2(iip1),dlat1(jjp1),dlat2(jjp1) + REAL acoslat,dxkm,dykm,resol(iip1,jjp1) + + rad = 6400000 + omeg = 7.272205e-05 + g = 9.8 + kappa = 0.285716 + daysec = 86400 + cpp = 1004.70885 + + preff = 101325. + pa= 50000. + + CALL conf_gcm( 99, .TRUE. , clesphy0 ) + CALL iniconst + CALL inigeom + + DO j=1,jjp1 + rlatudeg(j)=rlatu(j)*180./pi + ENDDO + DO j=1,jjm + rlatvdeg(j)=rlatv(j)*180./pi + ENDDO + + DO i=1,iip1 + rlonudeg(i)=rlonu(i)*180./pi + 360. + rlonvdeg(i)=rlonv(i)*180./pi + 360. + ENDDO + + + ! 2 ----- OUVERTURE DE LA SORTIE NETCDF + ! --------------------------------------------------- + ! CREATION OUTPUT + ! ouverture fichier netcdf de sortie out + status=NF_CREATE('grilles_gcm.nc',NF_NOCLOBBER,ncid_out) + status=NF_DEF_DIM(ncid_out,'lonu',iim+1,out_lonudim) + status=NF_DEF_DIM(ncid_out,'lonv',iim+1,out_lonvdim) + status=NF_DEF_DIM(ncid_out,'latu',jjm+1,out_latudim) + status=NF_DEF_DIM(ncid_out,'latv',jjm,out_latvdim) + + + ! Longitudes en u + status=NF_DEF_VAR(ncid_out,'lonu',NF_FLOAT,1,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',9,'Longitude en u') + + ! Longitudes en v + status=NF_DEF_VAR(ncid_out,'lonv',NF_FLOAT,1,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', 9,'Longitude en v') + + ! Latitude en u + status=NF_DEF_VAR(ncid_out,'latu',NF_FLOAT,1,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', 8,'Latitude en u') + + ! Latitude en v + status=NF_DEF_VAR(ncid_out,'latv',NF_FLOAT,1,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', 8,'Latitude en v') + + ! ecriture de la grille u + out_dim(1)=out_lonudim + out_dim(2)=out_latudim + status=NF_DEF_VAR(ncid_out,'grille_u',NF_FLOAT,2,out_dim, out_varid) + CALL handle_err(status) + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'units', 6,'Kelvin') + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'long_name', 16,'Grille aux point u') + + ! ecriture de la grille v + out_dim(1)=out_lonvdim + out_dim(2)=out_latvdim + status=NF_DEF_VAR(ncid_out,'grille_v',NF_FLOAT,2,out_dim, out_varid) + CALL handle_err(status) + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'units', 6,'Kelvin') + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'long_name', 16,'Grille aux point v') + + ! ecriture de la grille u + out_dim(1)=out_lonvdim + out_dim(2)=out_latudim + status=NF_DEF_VAR(ncid_out,'grille_s',NF_FLOAT,2,out_dim, out_varid) + CALL handle_err(status) + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'units', 6,'Kelvin') + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'long_name',16,'Grille aux point u') + + status=NF_ENDDEF(ncid_out) + ! 5) ----- FERMETURE DES FICHIERS NETCDF------------------ + ! -------------------------------------------------------- + ! 3-b- Ecriture de la grille pour la sortie + ! rajoute l'ecriture de la grille + +#ifdef NC_DOUBLE + status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonuid,1,iim+1,rlonudeg) + status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonvid,1,iim+1,rlonvdeg) + status=NF_PUT_VARA_DOUBLE(ncid_out,out_latuid,1,jjm+1,rlatudeg) + status=NF_PUT_VARA_DOUBLE(ncid_out,out_latvid,1,jjm,rlatvdeg) +#else + status=NF_PUT_VARA_REAL(ncid_out,out_lonuid,1,iim+1,rlonudeg) + status=NF_PUT_VARA_REAL(ncid_out,out_lonvid,1,iim+1,rlonvdeg) + status=NF_PUT_VARA_REAL(ncid_out,out_latuid,1,jjm+1,rlatudeg) + status=NF_PUT_VARA_REAL(ncid_out,out_latvid,1,jjm,rlatvdeg) +#endif + + start(1)=1 + start(2)=1 + start(3)=1 + start(4)=1 + + COUNT(1)=iim+1 + COUNT(2)=jjm+1 + COUNT(3)=1 + COUNT(4)=1 + + DO j=1,jjm+1 + DO i=1,iim+1 + temp(i,j)=MOD(i,2)+MOD(j,2) + ENDDO + ENDDO + +#ifdef NC_DOUBLE + status=NF_PUT_VARA_DOUBLE(ncid_out,out_varid,start, count,temp) +#else + status=NF_PUT_VARA_REAL(ncid_out,out_varid,start, count,temp) +#endif + + + ! fermeture du fichier netcdf + CALL ncclos(ncid_out,rcode_out) + +END SUBROUTINE grilles_gcm_netcdf_sub + + + +SUBROUTINE handle_err(status) + INCLUDE "netcdf.inc" + + INTEGER status + IF (status.NE.nf_noerr) THEN + PRINT *,NF_STRERROR(status) + CALL abort_gcm('grilles_gcm_netcdf','netcdf error',1) + ENDIF +END SUBROUTINE handle_err + Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/logic.h =================================================================== --- LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/logic.h (revision 1487) +++ LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/logic.h (revision 1488) @@ -10,10 +10,10 @@ & statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus & & ,read_start,ok_guide,ok_strato,ok_gradsfile & - & ,ok_limit,ok_etat0 + & ,ok_limit,ok_etat0,grilles_gcm_netcdf LOGICAL purmats,forward,leapf,apphys,statcl,conser, & & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus & & ,read_start,ok_guide,ok_strato,ok_gradsfile & - & ,ok_limit,ok_etat0 + & ,ok_limit,ok_etat0,grilles_gcm_netcdf INTEGER iflag_phys Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/ce0l.F90 =================================================================== --- LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/ce0l.F90 (revision 1487) +++ LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/ce0l.F90 (revision 1488) @@ -102,4 +102,13 @@ END IF + IF (grilles_gcm_netcdf) THEN + WRITE(lunout,'(//)') + WRITE(lunout,*) ' *************************** ' + WRITE(lunout,*) ' *** grilles_gcm_netcdf *** ' + WRITE(lunout,*) ' *************************** ' + WRITE(lunout,'(//)') + CALL grilles_gcm_netcdf_sub() + END IF + !$OMP MASTER CALL finalize_parallel Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/conf_gcm.F =================================================================== --- LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/conf_gcm.F (revision 1487) +++ LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/conf_gcm.F (revision 1488) @@ -888,4 +888,10 @@ ok_etat0 = .TRUE. CALL getin('ok_etat0',ok_etat0) + +!Config Key = grilles_gcm_netcdf +!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit +!Config Def = n + grilles_gcm_netcdf = .FALSE. + CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf) write(lunout,*)' #########################################' @@ -937,4 +943,5 @@ write(lunout,*)' ok_limit = ', ok_limit write(lunout,*)' ok_etat0 = ', ok_etat0 + write(lunout,*)' grilles_gcm_netcdf = ', grilles_gcm_netcdf c RETURN Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/grilles_gcm_netcdf_sub.F90 =================================================================== --- LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/grilles_gcm_netcdf_sub.F90 (revision 1488) +++ LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/grilles_gcm_netcdf_sub.F90 (revision 1488) @@ -0,0 +1,184 @@ +! +! $Header$ +! +! 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 +! part in the program create_fausse_var. +! +SUBROUTINE grilles_gcm_netcdf_sub + + IMPLICIT NONE + + INCLUDE "dimensions.h" + INCLUDE "paramet.h" + INCLUDE "comconst.h" + INCLUDE "comgeom.h" + INCLUDE "comvert.h" + INCLUDE "netcdf.inc" + INCLUDE "serre.h" + + 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_varid + INTEGER out_lonudim,out_lonvdim + INTEGER out_latudim,out_latvdim,out_dim(3) + + INTEGER, PARAMETER :: longcles = 20 + REAL clesphy0(longcles) + + INTEGER start(4),COUNT(4) + + INTEGER status,i,j + REAL rlatudeg(jjp1),rlatvdeg(jjm) + REAL rlonudeg(iip1),rlonvdeg(iip1) + + REAL dlon1(iip1),dlon2(iip1),dlat1(jjp1),dlat2(jjp1) + REAL acoslat,dxkm,dykm,resol(iip1,jjp1) + + rad = 6400000 + omeg = 7.272205e-05 + g = 9.8 + kappa = 0.285716 + daysec = 86400 + cpp = 1004.70885 + + preff = 101325. + pa= 50000. + + CALL conf_gcm( 99, .TRUE. , clesphy0 ) + CALL iniconst + CALL inigeom + + DO j=1,jjp1 + rlatudeg(j)=rlatu(j)*180./pi + ENDDO + DO j=1,jjm + rlatvdeg(j)=rlatv(j)*180./pi + ENDDO + + DO i=1,iip1 + rlonudeg(i)=rlonu(i)*180./pi + 360. + rlonvdeg(i)=rlonv(i)*180./pi + 360. + ENDDO + + + ! 2 ----- OUVERTURE DE LA SORTIE NETCDF + ! --------------------------------------------------- + ! CREATION OUTPUT + ! ouverture fichier netcdf de sortie out + status=NF_CREATE('grilles_gcm.nc',NF_NOCLOBBER,ncid_out) + status=NF_DEF_DIM(ncid_out,'lonu',iim+1,out_lonudim) + status=NF_DEF_DIM(ncid_out,'lonv',iim+1,out_lonvdim) + status=NF_DEF_DIM(ncid_out,'latu',jjm+1,out_latudim) + status=NF_DEF_DIM(ncid_out,'latv',jjm,out_latvdim) + + + ! Longitudes en u + status=NF_DEF_VAR(ncid_out,'lonu',NF_FLOAT,1,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',9,'Longitude en u') + + ! Longitudes en v + status=NF_DEF_VAR(ncid_out,'lonv',NF_FLOAT,1,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', 9,'Longitude en v') + + ! Latitude en u + status=NF_DEF_VAR(ncid_out,'latu',NF_FLOAT,1,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', 8,'Latitude en u') + + ! Latitude en v + status=NF_DEF_VAR(ncid_out,'latv',NF_FLOAT,1,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', 8,'Latitude en v') + + ! ecriture de la grille u + out_dim(1)=out_lonudim + out_dim(2)=out_latudim + status=NF_DEF_VAR(ncid_out,'grille_u',NF_FLOAT,2,out_dim, out_varid) + CALL handle_err(status) + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'units', 6,'Kelvin') + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'long_name', 16,'Grille aux point u') + + ! ecriture de la grille v + out_dim(1)=out_lonvdim + out_dim(2)=out_latvdim + status=NF_DEF_VAR(ncid_out,'grille_v',NF_FLOAT,2,out_dim, out_varid) + CALL handle_err(status) + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'units', 6,'Kelvin') + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'long_name', 16,'Grille aux point v') + + ! ecriture de la grille u + out_dim(1)=out_lonvdim + out_dim(2)=out_latudim + status=NF_DEF_VAR(ncid_out,'grille_s',NF_FLOAT,2,out_dim, out_varid) + CALL handle_err(status) + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'units', 6,'Kelvin') + status=NF_PUT_ATT_TEXT(ncid_out,out_varid,'long_name',16,'Grille aux point u') + + status=NF_ENDDEF(ncid_out) + ! 5) ----- FERMETURE DES FICHIERS NETCDF------------------ + ! -------------------------------------------------------- + ! 3-b- Ecriture de la grille pour la sortie + ! rajoute l'ecriture de la grille + +#ifdef NC_DOUBLE + status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonuid,1,iim+1,rlonudeg) + status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonvid,1,iim+1,rlonvdeg) + status=NF_PUT_VARA_DOUBLE(ncid_out,out_latuid,1,jjm+1,rlatudeg) + status=NF_PUT_VARA_DOUBLE(ncid_out,out_latvid,1,jjm,rlatvdeg) +#else + status=NF_PUT_VARA_REAL(ncid_out,out_lonuid,1,iim+1,rlonudeg) + status=NF_PUT_VARA_REAL(ncid_out,out_lonvid,1,iim+1,rlonvdeg) + status=NF_PUT_VARA_REAL(ncid_out,out_latuid,1,jjm+1,rlatudeg) + status=NF_PUT_VARA_REAL(ncid_out,out_latvid,1,jjm,rlatvdeg) +#endif + + start(1)=1 + start(2)=1 + start(3)=1 + start(4)=1 + + COUNT(1)=iim+1 + COUNT(2)=jjm+1 + COUNT(3)=1 + COUNT(4)=1 + + DO j=1,jjm+1 + DO i=1,iim+1 + temp(i,j)=MOD(i,2)+MOD(j,2) + ENDDO + ENDDO + +#ifdef NC_DOUBLE + status=NF_PUT_VARA_DOUBLE(ncid_out,out_varid,start, count,temp) +#else + status=NF_PUT_VARA_REAL(ncid_out,out_varid,start, count,temp) +#endif + + + ! fermeture du fichier netcdf + CALL ncclos(ncid_out,rcode_out) + +END SUBROUTINE grilles_gcm_netcdf_sub + + + +SUBROUTINE handle_err(status) + INCLUDE "netcdf.inc" + + INTEGER status + IF (status.NE.nf_noerr) THEN + PRINT *,NF_STRERROR(status) + CALL abort_gcm('grilles_gcm_netcdf','netcdf error',1) + ENDIF +END SUBROUTINE handle_err + Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/logic.h =================================================================== --- LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/logic.h (revision 1487) +++ LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/logic.h (revision 1488) @@ -10,10 +10,10 @@ & statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus & & ,read_start,ok_guide,ok_strato,ok_gradsfile & - & ,ok_limit,ok_etat0 + & ,ok_limit,ok_etat0,grilles_gcm_netcdf LOGICAL purmats,forward,leapf,apphys,statcl,conser, & & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus & & ,read_start,ok_guide,ok_strato,ok_gradsfile & - & ,ok_limit,ok_etat0 + & ,ok_limit,ok_etat0,grilles_gcm_netcdf INTEGER iflag_phys