Changeset 5100 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem
- Timestamp:
- Jul 23, 2024, 7:00:20 AM (2 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/dyn3dmem
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynetat0_loc.F90
r5099 r5100 9 9 USE infotrac, ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName 10 10 USE strings_mod, ONLY: maxlen, msg, strStack, real2str, int2str, strIdx 11 USE netcdf, ONLY: NF90_OPEN, NF90_NOWRITE, nf90_inquire_dimension, NF90_INQ_VARID, &12 NF90_CLOSE, nf90_get_var, NF90_INQUIRE_VARIABLE, nf90_noerr11 USE netcdf, ONLY: nf90_open, nf90_nowrite, nf90_inquire_dimension, NF90_INQ_VARID, & 12 nf90_close, nf90_get_var, NF90_INQUIRE_VARIABLE, nf90_noerr 13 13 USE readTracFiles_mod, ONLY: new2oldH2O, newHNO3, oldHNO3, getKey 14 14 USE control_mod, ONLY: planet_type … … 55 55 !--- Initial state file opening 56 56 var=fichnom 57 CALL err( NF90_OPEN(var,NF90_NOWRITE,fID),"open",var)57 CALL err(nf90_open(var,nf90_nowrite,fID),"open",var) 58 58 CALL get_var1("controle",tab_cntrl) 59 59 … … 215 215 END DO 216 216 DEALLOCATE(q_glo) 217 CALL err( NF90_CLOSE(fID),"close",fichnom)217 CALL err(nf90_close(fID),"close",fichnom) 218 218 day_ini=day_ini+INT(time) 219 219 time=time-INT(time) -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynredem_loc.F90
r5099 r5100 11 11 USE strings_mod, ONLY: maxlen 12 12 USE infotrac, ONLY: nqtot, tracers 13 USE netcdf, ONLY: NF90_CREATE, NF90_DEF_DIM, NF90_INQ_VARID, NF90_GLOBAL, &14 NF90_CLOSE, NF90_PUT_ATT, NF90_UNLIMITED, NF90_CLOBBER, &15 NF90_64BIT_OFFSET13 USE netcdf, ONLY: nf90_create, nf90_def_dim, NF90_INQ_VARID, nf90_global, & 14 nf90_close, nf90_put_att, nf90_unlimited, nf90_clobber, & 15 nf90_64bit_offset 16 16 USE dynredem_mod, ONLY: cre_var, put_var, err, modname, fil 17 17 USE comvert_mod, ONLY: ap, bp, presnivs, pa, preff, nivsig, nivsigs,& … … 112 112 113 113 !--- File creation 114 CALL err( NF90_CREATE(fichnom,IOR(NF90_CLOBBER,NF90_64BIT_OFFSET),nid))114 CALL err(nf90_create(fichnom,IOR(nf90_clobber,nf90_64bit_offset),nid)) 115 115 116 116 !--- Some global attributes 117 CALL err( NF90_PUT_ATT(nid,NF90_GLOBAL,"title","Fichier demarrage dynamique"))117 CALL err(nf90_put_att(nid,nf90_global,"title","Fichier demarrage dynamique")) 118 118 119 119 !--- Dimensions 120 CALL err( NF90_DEF_DIM(nid,"index", length, indexID))121 CALL err( NF90_DEF_DIM(nid,"rlonu", iip1, rlonuID))122 CALL err( NF90_DEF_DIM(nid,"rlatu", jjp1, rlatuID))123 CALL err( NF90_DEF_DIM(nid,"rlonv", iip1, rlonvID))124 CALL err( NF90_DEF_DIM(nid,"rlatv", jjm, rlatvID))125 CALL err( NF90_DEF_DIM(nid,"sigs", llm, sID))126 CALL err( NF90_DEF_DIM(nid,"sig", llmp1, sigID))127 CALL err( NF90_DEF_DIM(nid,"temps", NF90_UNLIMITED, timID))120 CALL err(nf90_def_dim(nid,"index", length, indexID)) 121 CALL err(nf90_def_dim(nid,"rlonu", iip1, rlonuID)) 122 CALL err(nf90_def_dim(nid,"rlatu", jjp1, rlatuID)) 123 CALL err(nf90_def_dim(nid,"rlonv", iip1, rlonvID)) 124 CALL err(nf90_def_dim(nid,"rlatv", jjm, rlatvID)) 125 CALL err(nf90_def_dim(nid,"sigs", llm, sID)) 126 CALL err(nf90_def_dim(nid,"sig", llmp1, sigID)) 127 CALL err(nf90_def_dim(nid,"temps", nf90_unlimited, timID)) 128 128 129 129 !--- Define and save invariant fields … … 156 156 CALL cre_var(nid,"masse","Masse d air" ,[rlonvID,rlatuID,sID,timID]) 157 157 CALL cre_var(nid,"ps" ,"Pression au sol",[rlonvID,rlatuID ,timID]) 158 CALL err( NF90_CLOSE(nid))158 CALL err(nf90_close (nid)) 159 159 160 160 WRITE(lunout,*)TRIM(modname)//': iim,jjm,llm,iday_end',iim,jjm,llm,iday_end … … 178 178 USE infotrac, ONLY: nqtot, tracers, type_trac 179 179 USE control_mod 180 USE netcdf, ONLY: NF90_OPEN, NF90_NOWRITE, nf90_get_var, NF90_INQ_VARID, &181 NF90_CLOSE, NF90_WRITE, NF90_PUT_VAR, nf90_noerr180 USE netcdf, ONLY: nf90_open, nf90_nowrite, nf90_get_var, NF90_INQ_VARID, & 181 nf90_close, NF90_WRITE, nf90_put_var, nf90_noerr 182 182 USE dynredem_mod, ONLY: dynredem_write_u, dynredem_write_v, dynredem_read_u, & 183 183 err, modname, fil, msg … … 213 213 IF(mpi_rank==0) THEN !++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 214 214 modname='dynredem1_loc'; fil=fichnom 215 CALL err( NF90_OPEN(fil,NF90_WRITE,nid),"open",fil)215 CALL err(nf90_open(fil,NF90_WRITE,nid),"open",fil) 216 216 217 217 !--- Write/extend time coordinate … … 219 219 var="temps" 220 220 CALL err(NF90_INQ_VARID(nid,var,vID),"inq",var) 221 CALL err( NF90_PUT_VAR(nid,vID,[time]),"put",var)221 CALL err(nf90_put_var(nid,vID,[time]),"put",var) 222 222 WRITE(lunout,*)TRIM(modname)//": Saving for ", nb, time 223 223 … … 228 228 tab_cntrl(31)=DBLE(itau_dyn + itaufin) 229 229 CALL err(NF90_INQ_VARID(nid,var,vID),"inq",var) 230 CALL err( NF90_PUT_VAR(nid,vID,tab_cntrl),"put",var)230 CALL err(nf90_put_var(nid,vID,tab_cntrl),"put",var) 231 231 END IF !++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 232 232 !$OMP END MASTER … … 244 244 fil="start_trac.nc" 245 245 IF(ANY(type_trac == ['inca','inco'])) INQUIRE(FILE=fil,EXIST=lread_inca) 246 IF(lread_inca) CALL err( NF90_OPEN(fil,NF90_NOWRITE,nid_trac),"open")246 IF(lread_inca) CALL err(nf90_open(fil,nf90_nowrite,nid_trac),"open") 247 247 !$OMP END MASTER 248 248 !$OMP BARRIER … … 266 266 !$OMP MASTER 267 267 IF(mpi_rank==0) THEN !++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 268 CALL err( NF90_CLOSE(nid),"close")268 CALL err(nf90_close(nid),"close") 269 269 fil="start_trac.nc" 270 IF(lread_inca) CALL err( NF90_CLOSE(nid_trac),"close")270 IF(lread_inca) CALL err(nf90_close(nid_trac),"close") 271 271 END IF !++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 272 272 !$OMP END MASTER -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynredem_mod.F90
r5099 r5100 52 52 !$OMP MASTER 53 53 start(3)=l 54 CALL err( NF90_PUT_VAR(ncid,nvarid,var_glo,start,count),"put",id)54 CALL err(nf90_put_var(ncid,nvarid,var_glo,start,count),"put",id) 55 55 !$OMP END MASTER 56 56 END IF … … 102 102 !$OMP MASTER 103 103 start(3)=l 104 CALL err( NF90_PUT_VAR(ncid,nvarid,var_glo,start,count),"put",id)104 CALL err(nf90_put_var(ncid,nvarid,var_glo,start,count),"put",id) 105 105 !$OMP END MASTER 106 106 END IF … … 182 182 CHARACTER(LEN=*), OPTIONAL, INTENT(IN) :: units 183 183 !=============================================================================== 184 CALL err( NF90_DEF_VAR(ncid,var,nf90_format,did,nvarid),"inq",var)185 IF(title/="") CALL err( NF90_PUT_ATT(ncid,nvarid,"title",title),var)186 IF(PRESENT(units)) CALL err( NF90_PUT_ATT(ncid,nvarid,"units",units),var)184 CALL err(nf90_def_var(ncid,var,nf90_format,did,nvarid),"inq",var) 185 IF(title/="") CALL err(nf90_put_att(ncid,nvarid,"title",title),var) 186 IF(PRESENT(units)) CALL err(nf90_put_att(ncid,nvarid,"units",units),var) 187 187 188 188 END SUBROUTINE cre_var … … 208 208 IF( PRESENT(units)) CALL cre_var(ncid,var,title,did,units) 209 209 IF(.NOT.PRESENT(units)) CALL cre_var(ncid,var,title,did) 210 CALL err( NF90_ENDDEF(ncid))210 CALL err(nf90_enddef(ncid)) 211 211 nd=SIZE(did) 212 212 DO k=1,nd; CALL err(nf90_inquire_dimension(ncid,did(k),len=nn(k))); END DO 213 IF(nd==1) CALL err( NF90_PUT_VAR(ncid,nvarid,RESHAPE(v,nn(1:1))),var)214 IF(nd==2) CALL err( NF90_PUT_VAR(ncid,nvarid,RESHAPE(v,nn(1:2))),var)213 IF(nd==1) CALL err(nf90_put_var(ncid,nvarid,RESHAPE(v,nn(1:1))),var) 214 IF(nd==2) CALL err(nf90_put_var(ncid,nvarid,RESHAPE(v,nn(1:2))),var) 215 215 CALL err(NF90_REDEF(ncid)) 216 216 END SUBROUTINE put_var -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/guide_loc_mod.F90
r5099 r5100 175 175 if (guide_plevs==1) then 176 176 if (ncidpl==-99) then 177 rcod=nf90_open('apbp.nc', Nf90_NOWRITe, ncidpl)177 rcod=nf90_open('apbp.nc',nf90_nowrite, ncidpl) 178 178 if (rcod/=nf90_noerr) THEN 179 179 abort_message=' Nudging error -> no file apbp.nc' … … 183 183 elseif (guide_plevs==2) then 184 184 if (ncidpl==-99) then 185 rcod=nf90_open('P.nc', Nf90_NOWRITe,ncidpl)185 rcod=nf90_open('P.nc',nf90_nowrite,ncidpl) 186 186 if (rcod/=nf90_noerr) THEN 187 187 abort_message=' Nudging error -> no file P.nc' … … 192 192 elseif (guide_u) then 193 193 if (ncidpl==-99) then 194 rcod=nf90_open('u.nc', Nf90_NOWRITe,ncidpl)194 rcod=nf90_open('u.nc',nf90_nowrite,ncidpl) 195 195 if (rcod/=nf90_noerr) THEN 196 196 abort_message=' Nudging error -> no file u.nc' -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/iniacademic_loc.F90
r5099 r5100 23 23 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 24 24 USE readTracFiles_mod, ONLY: addPhase 25 use netcdf, only : NF90_NOWRITE,NF90_OPEN,nf90_noerr,NF90_INQ_VARID,NF90_CLOSE, nf90_get_var25 use netcdf, only : nf90_nowrite,nf90_open,nf90_noerr,NF90_INQ_VARID,nf90_close, nf90_get_var 26 26 27 27 ! Author: Frederic Hourdin original: 15/01/93 … … 152 152 153 153 relief=0. 154 ierr = NF90_OPEN ('relief_in.nc', NF90_NOWRITE,nid_relief)154 ierr = nf90_open ('relief_in.nc', nf90_nowrite,nid_relief) 155 155 if (ierr==nf90_noerr) THEN 156 156 ierr=NF90_INQ_VARID(nid_relief,'RELIEF',varid) … … 162 162 endif 163 163 endif 164 ierr = NF90_CLOSE(nid_relief)164 ierr = nf90_close (nid_relief) 165 165 166 166
Note: See TracChangeset
for help on using the changeset viewer.