Changeset 1011 for IOIPSL/trunk/src/histcom.f90
- Timestamp:
- 05/07/10 15:05:34 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
IOIPSL/trunk/src/histcom.f90
r1005 r1011 1128 1128 IF (pzid > W_F(idf)%n_zax) THEN 1129 1129 WRITE(str70, & 1130 & '("The vertical coordinate chosen for variable ", a)') &1130 & '("The vertical coordinate chosen for variable ",A)') & 1131 1131 & TRIM(tmp_name) 1132 1132 str71 = " Does not exist." … … 1300 1300 IF (l_dbg) WRITE(*,*) "histdef : 6.0" 1301 1301 !- 1302 IF (W_F(idf)%W_V(iv)%freq_wrt > 0) THEN1303 WRITE(str10,'(I8.8)') INT(W_F(idf)%W_V(iv)%freq_wrt)1304 str40 = TRIM(W_F(idf)%W_V(iv)%topp)//"_"//TRIM(str10)1305 ELSE1306 WRITE(str10,'(I2.2,"month")') ABS(INT(W_F(idf)%W_V(iv)%freq_wrt))1307 str40 = TRIM(W_F(idf)%W_V(iv)%topp)//"_"//TRIM(str10)1308 ENDIF1309 CALL find_str (W_F(idf)%W_V(1:W_F(idf)%n_tax)%tax_name,str40,pos)1310 !-1311 1302 ! No time axis for once, l_max, l_min or never operation 1312 1303 !- … … 1315 1306 & .AND.(TRIM(W_F(idf)%W_V(iv)%topp) /= 'l_max') & 1316 1307 & .AND.(TRIM(W_F(idf)%W_V(iv)%topp) /= 'l_min') ) THEN 1308 IF (TRIM(W_F(idf)%W_V(iv)%topp) == 'inst') THEN 1309 str10 = 't_inst_' 1310 ELSE 1311 str10 = 't_op_' 1312 ENDIF 1313 IF (W_F(idf)%W_V(iv)%freq_wrt > 0) THEN 1314 WRITE (UNIT=str40,FMT='(A,I8.8)') & 1315 & TRIM(str10),INT(W_F(idf)%W_V(iv)%freq_wrt) 1316 ELSE 1317 WRITE (UNIT=str40,FMT='(A,I2.2,"month")') & 1318 & TRIM(str10),ABS(INT(W_F(idf)%W_V(iv)%freq_wrt)) 1319 ENDIF 1320 CALL find_str (W_F(idf)%W_V(1:W_F(idf)%n_tax)%tax_name,str40,pos) 1317 1321 IF (pos < 0) THEN 1318 1322 W_F(idf)%n_tax = W_F(idf)%n_tax+1 … … 1387 1391 dims(1) = W_F(idf)%tid 1388 1392 IF (W_F(idf)%n_tax > 1) THEN 1389 str30 = "t_"//W_F(idf)%W_V(itx)%tax_name1393 str30 = W_F(idf)%W_V(itx)%tax_name 1390 1394 ELSE 1391 1395 str30 = "time_counter" 1392 1396 ENDIF 1393 iret = NF90_DEF_VAR (nfid, str30,NF90_REAL8, &1394 & 1395 IF ( W_F(idf)%n_tax <= 1) THEN1397 iret = NF90_DEF_VAR (nfid,TRIM(str30),NF90_REAL8, & 1398 & dims(1),W_F(idf)%W_V(itx)%tdimid) 1399 IF (itx <= 1) THEN 1396 1400 iret = NF90_PUT_ATT (nfid,W_F(idf)%W_V(itx)%tdimid,'axis',"T") 1397 1401 ENDIF 1398 1402 iret = NF90_PUT_ATT (nfid,W_F(idf)%W_V(itx)%tdimid, & 1399 & 1403 & 'standard_name',"time") 1400 1404 !--- 1401 1405 ! To transform the current itau into a real date and take it
Note: See TracChangeset
for help on using the changeset viewer.