Changeset 3852 for LMDZ6/branches/LMDZ-tracers/libf/dyn3d/dynredem.F90
- Timestamp:
- Feb 22, 2021, 5:28:31 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ-tracers/libf/dyn3d/dynredem.F90
r3851 r3852 7 7 USE IOIPSL 8 8 #endif 9 USE infotrac 9 USE infotrac, ONLY: nqtot, tracers 10 10 USE netcdf, ONLY: NF90_CREATE, NF90_DEF_DIM, NF90_INQ_VARID, NF90_GLOBAL, & 11 11 NF90_CLOSE, NF90_PUT_ATT, NF90_UNLIMITED, NF90_CLOBBER, & … … 145 145 CALL cre_var(nid,"teta" ,"Temperature",[rlonvID,rlatuID,sID,timID]) 146 146 DO iq=1,nqtot 147 CALL cre_var(nid,t name(iq),ttext(iq),[rlonvID,rlatuID,sID,timID])147 CALL cre_var(nid,tracers(iq)%name,tracers(iq)%lnam,[rlonvID,rlatuID,sID,timID]) 148 148 END DO 149 149 CALL cre_var(nid,"masse","Masse d air" ,[rlonvID,rlatuID,sID,timID]) … … 166 166 ! Purpose: Write the NetCDF restart file (append). 167 167 !------------------------------------------------------------------------------- 168 USE infotrac 168 USE infotrac, ONLY: nqtot, tracers, type_trac 169 169 USE control_mod 170 170 USE netcdf, ONLY: NF90_OPEN, NF90_NOWRITE, NF90_GET_VAR, NF90_INQ_VARID, & … … 226 226 227 227 !--- Tracers in file "start_trac.nc" (added by Anne) 228 lread_inca=.FALSE.; fil="start_trac.nc" 229 IF(type_trac=='inca') INQUIRE(FILE=fil,EXIST=lread_inca) 228 fil="start_trac.nc" 229 INQUIRE(FILE=fil, EXIST=lread_inca) 230 lread_inca = lread_inca .AND. type_trac == 'inca' 230 231 IF(lread_inca) CALL err(NF90_OPEN(fil,NF90_NOWRITE,nid_trac),"open") 231 232 232 233 !--- Save tracers 233 DO iq=1,nqtot; var=t name(iq); ierr=-1234 DO iq=1,nqtot; var=tracers(iq)%name; ierr=-1 234 235 IF(lread_inca) THEN !--- Possibly read from "start_trac.nc" 235 236 fil="start_trac.nc" … … 237 238 dum='inq'; IF(ierr==NF90_NoErr) dum='fnd' 238 239 WRITE(lunout,*)msg(dum,var) 239 240 241 240 IF(ierr==NF90_NoErr) CALL dynredem_read_u(nid_trac,var,q(:,:,:,iq),llm) 242 241 END IF
Note: See TracChangeset
for help on using the changeset viewer.