Changeset 1443 for LMDZ5/branches/LMDZ5V1.0-dev/libf/dyn3dpar/infotrac.F90
- Timestamp:
- Oct 21, 2010, 10:23:18 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/LMDZ5V1.0-dev/libf/dyn3dpar/infotrac.F90
r1403 r1443 61 61 CHARACTER(len=1), DIMENSION(3) :: txts 62 62 CHARACTER(len=2), DIMENSION(9) :: txtp 63 CHARACTER(len= 13) :: str1,str263 CHARACTER(len=23) :: str1,str2 64 64 65 65 INTEGER :: nqtrue ! number of tracers read from tracer.def, without higer order of moment 66 66 INTEGER :: iq, new_iq, iiq, jq, ierr 67 INTEGER, EXTERNAL :: lnblnk68 67 69 68 !----------------------------------------------------------------------- … … 217 216 new_iq=new_iq+10 ! 9 tracers added 218 217 ELSE 219 WRITE(lunout,*) 'This choice of advection schema is not available' 218 WRITE(lunout,*) 'This choice of advection schema is not available',iq,hadv(iq),vadv(iq) 220 219 CALL abort_gcm('infotrac_init','Bad choice of advection schema - 1',1) 221 220 END IF … … 258 257 iadv(new_iq)=11 259 258 ELSE 260 WRITE(lunout,*)'This choice of advection schema is not available' 259 WRITE(lunout,*)'This choice of advection schema is not available',iq,hadv(iq),vadv(iq) 260 261 261 CALL abort_gcm('infotrac_init','Bad choice of advection schema - 2',1) 262 262 END IF … … 265 265 tname(new_iq)= tnom_0(iq) 266 266 IF (iadv(new_iq)==0) THEN 267 ttext(new_iq)= str1(1:lnblnk(str1))267 ttext(new_iq)=trim(str1) 268 268 ELSE 269 ttext(new_iq)= str1(1:lnblnk(str1))//descrq(iadv(new_iq))269 ttext(new_iq)=trim(tnom_0(iq))//descrq(iadv(new_iq)) 270 270 END IF 271 271 … … 276 276 new_iq=new_iq+1 277 277 iadv(new_iq)=-20 278 ttext(new_iq)= str2(1:lnblnk(str2))//txts(jq)279 tname(new_iq)= str1(1:lnblnk(str1))//txts(jq)278 ttext(new_iq)=trim(str2)//txts(jq) 279 tname(new_iq)=trim(str1)//txts(jq) 280 280 END DO 281 281 ELSE IF (iadv(new_iq)==30) THEN … … 283 283 new_iq=new_iq+1 284 284 iadv(new_iq)=-30 285 ttext(new_iq)= str2(1:lnblnk(str2))//txtp(jq)286 tname(new_iq)= str1(1:lnblnk(str1))//txtp(jq)285 ttext(new_iq)=trim(str2)//txtp(jq) 286 tname(new_iq)=trim(str1)//txtp(jq) 287 287 END DO 288 288 END IF
Note: See TracChangeset
for help on using the changeset viewer.