Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/dynredem_mod.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dyn3d/dynredem_mod.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dyn3d/dynredem_mod.F90 (revision 5090) @@ -3,5 +3,5 @@ USE netcdf, ONLY: nf90_strerror,nf90_noerr,nf90_redef,nf90_put_var,nf90_enddef,nf90_put_att,& nf90_inq_varid,nf90_get_var,nf90_def_var - USE lmdz_netcdf_format, ONLY: nf90_format + USE lmdz_cppkeys_wrapper, ONLY: nf90_format IMPLICIT NONE; PRIVATE PUBLIC :: dynredem_write_u, dynredem_write_v, dynredem_read_u, err Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/gcm.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dyn3d/gcm.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dyn3d/gcm.F90 (revision 5090) @@ -26,11 +26,12 @@ USE logic_mod, ONLY: ecripar, iflag_phys, read_start + USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS + !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique ! A nettoyer. On ne veut qu'une ou deux routines d'interface ! dynamique -> physique pour l'initialisation -#ifdef CPP_PHYS + ! AB 2024/07/20: remplace CPP key by fortran logical, but ^ still relevant, see later use of iniphys later on USE iniphysiq_mod, ONLY: iniphysiq -#endif !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! @@ -406,12 +407,12 @@ IF ((iflag_phys==1).or.(iflag_phys>=100)) THEN - ! Physics: -#ifdef CPP_PHYS - CALL iniphysiq(iim,jjm,llm, & + ! Physics: + IF (CPPKEY_PHYS) THEN + CALL iniphysiq(iim,jjm,llm, & (jjm-1)*iim+2,comm_lmdz, & daysec,day_ini,dtphys/nsplit_phys, & rlatu,rlatv,rlonu,rlonv,aire,cu,cv,rad,g,r,cpp, & iflag_phys) -#endif + END IF ENDIF ! of IF ((iflag_phys==1).or.(iflag_phys>=100)) Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/leapfrog.F =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dyn3d/leapfrog.F (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dyn3d/leapfrog.F (revision 5090) @@ -27,4 +27,5 @@ & start_time,dt USE strings_mod, ONLY: msg + USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS IMPLICIT NONE @@ -462,10 +463,10 @@ #endif ! #endif of #ifdef CPP_IOIPSL -#ifdef CPP_PHYS - CALL calfis( lafin , jD_cur, jH_cur, + IF (CPPKEY_PHYS) THEN + CALL calfis( lafin , jD_cur, jH_cur, $ ucov,vcov,teta,q,masse,ps,p,pk,phis,phi , $ du,dv,dteta,dq, $ flxw,dufi,dvfi,dtetafi,dqfi,dpfi ) -#endif + END IF c ajout des tendances physiques: c ------------------------------ Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/replay3d.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dyn3d/replay3d.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dyn3d/replay3d.F90 (revision 5090) @@ -22,20 +22,9 @@ USE mod_const_mpi, ONLY: comm_lmdz - - -!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! - ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique - ! A nettoyer. On ne veut qu'une ou deux routines d'interface - ! dynamique -> physique pour l'initialisation -#ifdef CPP_PHYS - USE iniphysiq_mod, ONLY: iniphysiq -#endif -!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! - IMPLICIT NONE ! ...... Version du 10/01/98 .......... - ! avec coordonnees verticales hybrides + ! avec coordonnees verticales hybrides ! avec nouveaux operat. dissipation * ( gradiv2,divgrad2,nxgraro2 ) Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bands.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bands.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bands.F90 (revision 5090) @@ -381,8 +381,9 @@ subroutine AdjustBands_physic use times -#ifdef CPP_PHYS -! Ehouarn: what follows is only related to // physics + + ! Ehouarn: what follows is only related to // physics + USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS USE mod_phys_lmdz_para, only : klon_mpi_para_nb -#endif + USE parallel_lmdz implicit none @@ -408,5 +409,5 @@ medium=medium/mpi_size NbTot=0 -#ifdef CPP_PHYS + IF (CPPKEY_PHYS) THEN do i=0,mpi_size-1 Inc(i)=nint(klon_mpi_para_nb(i)*(medium-value(i))/value(i)) @@ -431,5 +432,5 @@ distrib_phys(i)=klon_mpi_para_nb(i)+inc(i) enddo -#endif + END IF end subroutine AdjustBands_physic Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynredem_mod.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynredem_mod.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynredem_mod.F90 (revision 5090) @@ -6,5 +6,5 @@ USE netcdf, ONLY: nf90_strerror,nf90_noerr,nf90_redef,nf90_put_var,nf90_inquire_dimension,& nf90_inq_varid,nf90_get_var,nf90_def_var,nf90_enddef,nf90_put_att - USE lmdz_netcdf_format, ONLY: nf90_format + USE lmdz_cppkeys_wrapper, ONLY: nf90_format PRIVATE PUBLIC :: dynredem_write_u, dynredem_write_v, dynredem_read_u, err Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/gcm.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/gcm.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/gcm.F90 (revision 5090) @@ -12,7 +12,4 @@ USE parallel_lmdz USE infotrac, ONLY: nqtot, init_infotrac -!#ifdef CPP_PHYS -! USE mod_interface_dyn_phys, ONLY: init_interface_dyn_phys -!#endif USE mod_hallo USE Bands @@ -20,7 +17,8 @@ USE control_mod -#ifdef CPP_PHYS + USE iniphysiq_mod, ONLY: iniphysiq -#endif + USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS + USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, g, r, rad USE logic_mod ! all of it, because of copyin clause when calling leapfrog @@ -415,11 +413,11 @@ IF ((iflag_phys==1).or.(iflag_phys>=100)) THEN ! Physics: -#ifdef CPP_PHYS - CALL iniphysiq(iim,jjm,llm, & - distrib_phys(mpi_rank),comm_lmdz, & - daysec,day_ini,dtphys/nsplit_phys, & - rlatu,rlatv,rlonu,rlonv,aire,cu,cv,rad,g,r,cpp, & - iflag_phys) -#endif + IF (CPPKEY_PHYS) THEN + CALL iniphysiq(iim,jjm,llm, & + distrib_phys(mpi_rank),comm_lmdz, & + daysec,day_ini,dtphys/nsplit_phys, & + rlatu,rlatv,rlonu,rlonv,aire,cu,cv,rad,g,r,cpp, & + iflag_phys) + END IF ENDIF ! of IF ((iflag_phys==1).or.(iflag_phys>=100)) Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/lmdz_call_calfis.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/lmdz_call_calfis.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/lmdz_call_calfis.F90 (revision 5090) @@ -88,7 +88,6 @@ USE comvert_mod, ONLY: ap, bp, pressure_exner USE temps_mod, ONLY: day_ini, day_ref, jd_ref, jh_ref, start_time -#ifdef CPP_PHYS USE lmdz_calfis_loc -#endif + USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS IMPLICIT NONE @@ -233,10 +232,10 @@ !$OMP BARRIER -#ifdef CPP_PHYS - CALL calfis_loc(lafin ,jD_cur, jH_cur, & + IF (CPPKEY_PHYS) THEN + CALL calfis_loc(lafin ,jD_cur, jH_cur, & ucov,vcov,teta,q,masse,ps,p,pk,phis,phi , & du,dv,dteta,dq, & flxw, dufi,dvfi,dtetafi,dqfi,dpfi ) -#endif + END IF ijb=ij_begin ije=ij_end Index: LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.F =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.F (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.F (revision 5090) @@ -27,14 +27,13 @@ $ pdpsfi) c -c Auteur : P. Le Van, F. Hourdin +c Auteur : P. Le Van, F. Hourdin c ......... USE infotrac, ONLY: nqtot, tracers USE control_mod, ONLY: planet_type, nsplit_phys -#ifdef CPP_PHYS USE callphysiq_mod, ONLY: call_physiq -#endif + USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, kappa, pi USE comvert_mod, ONLY: preff, presnivs - + IMPLICIT NONE c======================================================================= @@ -48,14 +47,14 @@ c ---------- c -c - les vents sont donnes dans la physique par leurs composantes +c - les vents sont donnes dans la physique par leurs composantes c naturelles. c - la variable thermodynamique de la physique est une variable -c intensive : T +c intensive : T c pour la dynamique on prend T * ( preff / p(l) ) **kappa c - les deux seules variables dependant de la geometrie necessaires -c pour la physique sont la latitude pour le rayonnement et -c l'aire de la maille quand on veut integrer une grandeur +c pour la physique sont la latitude pour le rayonnement et +c l'aire de la maille quand on veut integrer une grandeur c horizontalement. -c - les points de la physique sont les points scalaires de la +c - les points de la physique sont les points scalaires de la c la dynamique; numerotation: c 1 pour le pole nord @@ -67,5 +66,5 @@ c ------- c pucov covariant zonal velocity -c pvcov covariant meridional velocity +c pvcov covariant meridional velocity c pteta potential temperature c pps surface pressure @@ -77,5 +76,5 @@ c -------- c pdufi tendency for the natural zonal velocity (ms-1) -c pdvfi tendency for the natural meridional velocity +c pdvfi tendency for the natural meridional velocity c pdhfi tendency for the potential temperature c pdtsfi tendency for the surface temperature @@ -165,5 +164,5 @@ REAL flxwfi(ngridmx,llm) ! Flux de masse verticale sur la grille physiq c - + REAL SSUM @@ -201,5 +200,5 @@ c ---------------------------------- - + zpsrf(1) = pps(1,1) @@ -217,5 +216,5 @@ c ----------------------------------------------------------------- c .... zplev definis aux (llm +1) interfaces des couches .... -c .... zplay definis aux ( llm ) milieux des couches .... +c .... zplay definis aux ( llm ) milieux des couches .... c ----------------------------------------------------------------- @@ -329,6 +328,6 @@ c .... Calcul de la vitesse verticale ( en Pa*m*s ou Kg/s ) .... -c JG : ancien calcule de omega utilise dans physiq.F. Maintenant le flux -c de masse est calclue dans advtrac.F +c JG : ancien calcule de omega utilise dans physiq.F. Maintenant le flux +c de masse est calclue dans advtrac.F c DO l=1,llm c pvervel(1,l)=pw(1,1,l) * g /apoln @@ -351,7 +350,7 @@ DO j=2,jjm ig0 = 1+(j-2)*iim - zufi(ig0+1,l)= 0.5 * + zufi(ig0+1,l)= 0.5 * $ ( pucov(iim,j,l)/cu(iim,j) + pucov(1,j,l)/cu(1,j) ) - pcvgu(ig0+1,l)= 0.5 * + pcvgu(ig0+1,l)= 0.5 * $ ( pducov(iim,j,l)/cu(iim,j) + pducov(1,j,l)/cu(1,j) ) DO i=2,iim @@ -373,6 +372,6 @@ do j=1,jjm zrot(i,j,l) = (pvcov(i+1,j,l) - pvcov(i,j,l) - $ + pucov(i,j+1,l) - pucov(i,j,l)) - $ / (cu(i,j)+cu(i,j+1)) + $ + pucov(i,j+1,l) - pucov(i,j,l)) + $ / (cu(i,j)+cu(i,j+1)) $ / (cv(i+1,j)+cv(i,j)) *4 enddo @@ -402,5 +401,5 @@ -c 47. champs de vents aux pole nord +c 47. champs de vents aux pole nord c ------------------------------ c U = 1 / pi * integrale [ v * cos(long) * d long ] @@ -475,5 +474,5 @@ zdqfic(:,:,:)=0. -#ifdef CPP_PHYS + IF (CPPKEY_PHYS) THEN do isplit=1,nsplit_phys @@ -500,5 +499,5 @@ ! . nqtot, !! nq ! . tracers(:)%name,!! tracer names from dynamical core (given in infotrac) -! . debut_split, !! firstcall +! . debut_split, !! firstcall ! . lafin_split, !! lastcall ! . jD_cur, !! pday. see leapfrog @@ -534,6 +533,5 @@ enddo ! of do isplit=1,nsplit_phys -#endif -! of #ifdef CPP_PHYS + END IF zdufi(:,:)=zdufic(:,:)/nsplit_phys @@ -541,7 +539,4 @@ zdtfi(:,:)=zdtfic(:,:)/nsplit_phys zdqfi(:,:,:)=zdqfic(:,:,:)/nsplit_phys - - -500 CONTINUE c----------------------------------------------------------------------- @@ -681,7 +676,4 @@ c----------------------------------------------------------------------- - -700 CONTINUE - firstcal = .FALSE. Index: LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90 (revision 5090) @@ -28,9 +28,5 @@ pdqfi, & pdpsfi) -#ifdef CPP_PHYS - ! If using physics - ! - ! Auteur : P. Le Van, F. Hourdin - ! ......... + USE dimphy USE mod_phys_lmdz_mpi_data, mpi_root_xx=>mpi_master @@ -39,5 +35,5 @@ USE mod_interface_dyn_phys USE IOPHY -#endif + USE lmdz_mpi @@ -49,7 +45,7 @@ USE infotrac, ONLY : nqtot, tracers USE control_mod, ONLY : planet_type, nsplit_phys -#ifdef CPP_PHYS USE callphysiq_mod, ONLY: call_physiq -#endif + USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS + USE comvert_mod, ONLY : preff, presnivs USE comconst_mod, ONLY : cpp, daysec, dtphys, dtvr, kappa, pi @@ -146,5 +142,5 @@ REAL, INTENT(OUT) :: pdpsfi(iip1, jjb_u:jje_u) ! tendency on surface pressure (Pa/s) -#ifdef CPP_PHYS + ! Ehouarn: for now calfis_p needs some informations from physics to compile ! Local variables : @@ -237,4 +233,6 @@ INTEGER :: jjb,jje + IF (CPPKEY_PHYS) THEN + ! !----------------------------------------------------------------------- @@ -249,5 +247,5 @@ IF ( firstcal ) THEN debut = .TRUE. - IF (ngridmx.NE.2+(jjm-1)*iim) THEN + IF (ngridmx/=2+(jjm-1)*iim) THEN write(lunout,*) 'STOP dans calfis' write(lunout,*) & @@ -719,5 +717,4 @@ zdqfic_omp(:,:,:)=0. -#ifdef CPP_PHYS do isplit=1,nsplit_phys @@ -750,7 +747,4 @@ enddo -#endif - ! of #ifdef CPP_PHYS - zdufi_omp(:,:)=zdufic_omp(:,:)/nsplit_phys @@ -847,5 +841,4 @@ klon=klon_mpi -500 CONTINUE !$OMP BARRIER @@ -1204,9 +1197,6 @@ !----------------------------------------------------------------------- -700 CONTINUE - firstcal = .FALSE. -#endif - ! of #ifdef CPP_PHYS +END IF END SUBROUTINE calfis_loc Index: LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/limit_netcdf.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/limit_netcdf.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/limit_netcdf.F90 (revision 5090) @@ -76,5 +76,5 @@ NF90_CLOBBER, NF90_ENDDEF, NF90_UNLIMITED, NF90_FLOAT, & NF90_64BIT_OFFSET - USE lmdz_netcdf_format, ONLY: nf90_format + USE lmdz_cppkeys_wrapper, ONLY: nf90_format USE inter_barxy_m, ONLY: inter_barxy USE netcdf95, ONLY: nf95_def_var, nf95_put_att, nf95_put_var Index: LMDZ6/branches/Amaury_dev/libf/misc/lmdz_cppkeys_wrapper.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/misc/lmdz_cppkeys_wrapper.F90 (revision 5090) +++ LMDZ6/branches/Amaury_dev/libf/misc/lmdz_cppkeys_wrapper.F90 (revision 5090) @@ -0,0 +1,29 @@ +! --------------------------------------------- +! This file is part of an effort to replace most uses of preprocessor CPP keys by fortran variables, +! to improve readability, compilation coverage, and linting. +! CPP keys used here should ONLY be used here, +! and imported through USE ..., ONLY: ... elsewhere +! CPP keys supported (key -> fortran variables associated): +! NC_DOUBLE -> nf90_format +! CPP_PHYS -> CPPKEY_PHYS +! --------------------------------------------- + +MODULE lmdz_cppkeys_wrapper + USE, INTRINSIC :: ISO_FORTRAN_ENV, ONLY : REAL64, REAL32 + USE netcdf, ONLY: nf90_float, nf90_double + IMPLICIT NONE; PRIVATE + PUBLIC nf90_format, CPPKEY_PHYS + +#ifdef NC_DOUBLE + INTEGER, PARAMETER :: nf90_format = nf90_double +#else + INTEGER, PARAMETER :: nf90_format = nf90_float +#endif + +#ifdef CPP_PHYS + LOGICAL, PARAMETER :: CPPKEY_PHYS = .TRUE. +#else + LOGICAL, PARAMETER :: CPPKEY_PHYS = .FALSE. +#endif + +END MODULE lmdz_cppkeys_wrapper Index: LMDZ6/branches/Amaury_dev/libf/misc/lmdz_netcdf_format.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/misc/lmdz_netcdf_format.F90 (revision 5089) +++ (revision ) @@ -1,20 +1,0 @@ -! --------------------------------------------- -! This module serves as a wrapper around netcdf. -! It serves two primary functions: -! 1) Turn netcdf into a "real" fortran module, without the INCLUDE call -! 2) Handle the NC_DOUBLE CPP key. This key should ONLY be used here. -! The "real" netcdf module/headers should ONLY be called here. -! --------------------------------------------- - -MODULE lmdz_netcdf_format - USE, INTRINSIC :: ISO_FORTRAN_ENV, ONLY : REAL64, REAL32 - USE netcdf, ONLY: nf90_float, nf90_double - IMPLICIT NONE; PRIVATE - PUBLIC nf90_format - -#ifdef NC_DOUBLE - INTEGER, PARAMETER :: nf90_format = nf90_double -#else - INTEGER, PARAMETER :: nf90_format = nf90_float -#endif -END MODULE lmdz_netcdf_format Index: LMDZ6/branches/Amaury_dev/libf/misc/write_field.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/misc/write_field.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/misc/write_field.F90 (revision 5090) @@ -2,5 +2,5 @@ USE netcdf, ONLY: nf90_sync, nf90_put_var, nf90_enddef, nf90_def_dim, nf90_unlimited, & nf90_clobber, nf90_create, nf90_def_var - USE lmdz_netcdf_format, ONLY: nf90_format + USE lmdz_cppkeys_wrapper, ONLY: nf90_format implicit none Index: LMDZ6/branches/Amaury_dev/libf/phylmd/iostart.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/phylmd/iostart.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/phylmd/iostart.F90 (revision 5090) @@ -387,5 +387,5 @@ SUBROUTINE put_field_rgen(pass, field_name,title,field,field_size) USE netcdf, ONLY: nf90_def_var,nf90_put_att,nf90_inq_varid,nf90_put_var - USE lmdz_netcdf_format, ONLY: nf90_format + USE lmdz_cppkeys_wrapper, ONLY: nf90_format USE dimphy USE geometry_mod @@ -511,5 +511,5 @@ SUBROUTINE put_var_rgen(pass, var_name,title,var,var_size) USE netcdf, ONLY: nf90_def_var,nf90_put_var,nf90_inq_varid,nf90_put_att - USE lmdz_netcdf_format, ONLY: nf90_format + USE lmdz_cppkeys_wrapper, ONLY: nf90_format USE dimphy USE mod_phys_lmdz_para Index: LMDZ6/branches/Amaury_dev/libf/phylmd/phyaqua_mod.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/phylmd/phyaqua_mod.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/phylmd/phyaqua_mod.F90 (revision 5090) @@ -573,5 +573,5 @@ nf90_enddef, nf90_put_att, nf90_unlimited, nf90_noerr, nf90_global, nf90_clobber, & nf90_64bit_offset, nf90_def_dim, nf90_create - USE lmdz_netcdf_format, ONLY: nf90_format + USE lmdz_cppkeys_wrapper, ONLY: nf90_format IMPLICIT NONE Index: LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyaqua_mod.F90 =================================================================== --- LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyaqua_mod.F90 (revision 5089) +++ LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyaqua_mod.F90 (revision 5090) @@ -594,5 +594,5 @@ nf90_def_dim,nf90_create,nf90_put_var,nf90_unlimited,nf90_global,nf90_64bit_offset,& nf90_enddef - USE lmdz_netcdf_format, ONLY: nf90_format + USE lmdz_cppkeys_wrapper, ONLY: nf90_format IMPLICIT NONE