Changeset 33 for LMDZ.3.3/trunk/libf
- Timestamp:
- Feb 7, 2000, 4:43:34 PM (25 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/phylmd/ajsec.F
r2 r33 1 SUBROUTINE ajsec(paprs, pplay, t,q, d_t,d_q , itr, tr, d_tr)1 SUBROUTINE ajsec(paprs, pplay, t,q, d_t,d_q) 2 2 IMPLICIT none 3 3 c====================================================================== … … 16 16 REAL t(klon,klev), q(klon,klev) 17 17 REAL d_t(klon,klev), d_q(klon,klev) 18 c19 INTEGER it, itr20 REAL tr(klon,klev,nbtr), d_tr(klon,klev,nbtr)21 REAL ztr(klon,klev,nbtr), trm(nbtr)22 18 c 23 19 INTEGER limbas, limhau ! les couches a ajuster … … 45 41 ENDDO 46 42 ENDDO 47 IF (itr.GE.1) THEN48 DO it = 1, itr49 DO k = 1, klev50 DO i = 1, klon51 d_tr(i,k,it) = 0.052 ENDDO53 ENDDO54 ENDDO55 ENDIF56 43 c------------------------------------- detection des profils a modifier 57 44 DO k = limbas, limhau … … 62 49 ENDDO 63 50 ENDDO 64 c65 IF (itr.GE.1) THEN66 DO it = 1, itr67 DO k = limbas, limhau68 DO i = 1, klon69 ztr(i,k,it) = tr(i,k,it)70 ENDDO71 ENDDO72 ENDDO73 ENDIF74 51 c 75 52 DO k = limbas, limhau … … 102 79 hm = zh(i,k2) 103 80 qm = zq(i,k2) 104 IF (itr.GE.1) THEN105 DO it = 1, itr106 trm(it) = ztr(i,k2,it)107 ENDDO108 ENDIF109 81 8020 CONTINUE 110 82 sm = sm +zpkdp(i,k) 111 83 hm = hm +zpkdp(i,k) * (zh(i,k)-hm) / sm 112 84 qm = qm +zpkdp(i,k) * (zq(i,k)-qm) / sm 113 IF (itr.GE.1) THEN114 DO it = 1, itr115 trm(it) = trm(it) +116 . zpkdp(i,k) * (ztr(i,k,it)-trm(it)) / sm117 ENDDO118 ENDIF119 85 down = .FALSE. 120 86 IF (k1 .ne. limbas) THEN … … 137 103 zq(i,k) = qm 138 104 ENDDO 139 IF (itr.GE.1) THEN140 DO it = 1, itr141 DO k = k1, k2142 ztr(i,k,it) = trm(it)143 ENDDO144 ENDDO145 ENDIF146 105 k2 = k2 + 1 147 106 ENDIF … … 173 132 d_q(i,k) = 0.0 174 133 ENDDO 175 ENDDO176 ENDIF177 c178 IF (itr.GE.1) THEN179 DO it = 1, itr180 c181 DO k = limbas, limhau182 DO i = 1, klon183 d_tr(i,k,it) = ztr(i,k,it) - tr(i,k,it)184 ENDDO185 ENDDO186 c187 IF (limbas.GT.1) THEN188 DO k = 1, limbas-1189 DO i = 1, klon190 d_tr(i,k,it) = 0.0191 ENDDO192 ENDDO193 ENDIF194 c195 IF (limhau.LT.klev) THEN196 DO k = limhau+1, klev197 DO i = 1, klon198 d_tr(i,k,it) = 0.0199 ENDDO200 ENDDO201 ENDIF202 c203 134 ENDDO 204 135 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.