Changeset 127
- Timestamp:
- 08/09/07 10:17:53 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
IOIPSL/trunk/src/fliocom.f90
r10 r127 4950 4950 ENDIF 4951 4951 !--- 4952 !-- Rule 3 : we look for a correct "units" attribute 4953 !--- 4954 IF (i_v < 0) THEN 4955 SELECT CASE (c_ax) 4956 CASE('x') 4957 n_r = 2 4958 c_r(1)='degree_e'; c_r(2)='degrees_e'; 4959 CASE('y') 4960 n_r = 2 4961 c_r(1)='degree_n'; c_r(2)='degrees_n'; 4962 CASE('z') 4963 n_r = 3 4964 c_r(1)='m'; c_r(2)='km'; c_r(3)='hpa'; 4965 CASE('t') 4966 n_r = 6 4967 c_r(1)='week'; c_r(2)='day'; c_r(3)='hour'; 4968 c_r(4)='minute'; c_r(5)='second'; c_r(6)='timesteps'; 4969 END SELECT 4970 !----- 4971 L_R3: DO kv=1,nw_nv(f_i) 4972 i_rc = NF90_GET_ATT(nw_id(f_i),kv,'units',c_t1) 4973 IF (i_rc == NF90_NOERR) THEN 4974 CALL strlowercase (c_t1) 4975 i_rc = NF90_INQUIRE_VARIABLE(nw_id(f_i),kv,ndims=n_d) 4976 IF (n_d <= l_d) THEN 4977 DO k=1,n_r 4978 IF (INDEX(c_t1,TRIM(c_r(k))) == 1) THEN 4979 i_v = kv; nbd = n_d; 4980 EXIT L_R3 4981 ENDIF 4982 ENDDO 4983 ENDIF 4984 ENDIF 4985 ENDDO L_R3 4986 ENDIF 4987 !--- 4988 !-- Rule 4 : we look for a variable with one dimension 4952 !-- Rule 3 : we look for a variable with one dimension 4989 4953 !-- and which has the same name as its dimension 4990 4954 !--- … … 5009 4973 IF (i_rc == NF90_NOERR) THEN 5010 4974 CALL strlowercase (c_t1) 5011 L_R 4: DO kv=1,nw_nv(f_i)4975 L_R3: DO kv=1,nw_nv(f_i) 5012 4976 i_rc = NF90_INQUIRE_VARIABLE & 5013 4977 & (nw_id(f_i),kv,name=c_t2,ndims=n_d) … … 5016 4980 IF (TRIM(c_t1) == TRIM(c_t2)) THEN 5017 4981 i_v = kv; nbd = n_d; 5018 EXIT L_R 44982 EXIT L_R3 5019 4983 ENDIF 5020 4984 ENDIF 5021 ENDDO L_R 44985 ENDDO L_R3 5022 4986 ENDIF 5023 4987 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.