[297] | 1 | module output_eurasie_mod |
---|
| 2 | |
---|
| 3 | USE module3D_phy |
---|
| 4 | use bilan_eau_mod |
---|
| 5 | |
---|
| 6 | implicit none |
---|
| 7 | |
---|
| 8 | real :: bmean ! |
---|
| 9 | real :: accmean ! accumulation moyenne |
---|
| 10 | real :: ablmean ! ablation moyenne |
---|
| 11 | real :: calvmean ! moyenne calving |
---|
| 12 | real :: ablbordmean ! |
---|
| 13 | real :: ablatotmean |
---|
| 14 | real :: bmeltmean ! moyenne fusion basale |
---|
| 15 | real :: tbmean ! temperature basale moyenne |
---|
| 16 | real :: tbdotmean ! moyenne variation / temps temperature basale |
---|
| 17 | real :: vsmean ! vitesse de surface moyenne |
---|
| 18 | !real :: vsdotmean ! moyenne variation / temps vitesse de surface |
---|
| 19 | real :: uzsmean !!!! utilise ? ! vitesse verticale de surface moyenne |
---|
| 20 | real :: uzsdotmean ! moyenne variation / temps vitesse verticale de surface |
---|
| 21 | real :: uzkmean ! moyenne vitesse verticale de surface |
---|
| 22 | real :: hdotmean ! moyenne derivee / temps de H |
---|
| 23 | real :: bdotmean ! moyenne bedrock derive / temps |
---|
| 24 | !real :: pf0mean ! moyenne de PF0 |
---|
| 25 | !real :: pf1mean ! moyenne de PF1 |
---|
| 26 | !real :: evapmean ! moyenne de EVAP |
---|
| 27 | !real :: pwmean ! moyenne PW |
---|
| 28 | !real :: pdfmean ! moyenne PDF |
---|
[336] | 29 | integer, parameter :: nmask=7 ! nbr de masques |
---|
| 30 | logical,dimension(nx,ny,nmask) :: mask_cal ! masque regions calotte |
---|
[297] | 31 | integer, dimension(nx,ny) :: write_mask |
---|
| 32 | CONTAINS |
---|
| 33 | !_________________________________________________________________________ |
---|
| 34 | subroutine init_outshort |
---|
| 35 | |
---|
| 36 | !ndisp sorite courte tous les ndisp |
---|
| 37 | NDISP=100 |
---|
| 38 | mask_cal(:,:,:)=.false. |
---|
| 39 | mask_cal(:,:,1)=.true. |
---|
| 40 | ! calcul d'un masque pour les regions des calottes |
---|
| 41 | do j=1,ny |
---|
| 42 | do i=1,nx |
---|
| 43 | ! -- Fennoscandie |
---|
| 44 | IF( (((xlong(i,j).ge.345).AND.(xlong(i,j).lt.360)).AND. & !British |
---|
| 45 | ((ylat(i,j).ge.50).AND.(ylat(i,j).le.60))).OR. & |
---|
| 46 | (((xlong(i,j).ge.0).AND.(xlong(i,j).le.40)).AND. & !Scand |
---|
| 47 | ((ylat(i,j).ge.50).AND.(ylat(i,j).le.70))).OR. & |
---|
| 48 | ! (((xlong(i,j).ge.50).AND.(xlong(i,j).le.120)).AND. & |
---|
| 49 | (((xlong(i,j).ge.10).AND.(xlong(i,j).le.122)).AND. & !BK |
---|
| 50 | ((ylat(i,j).ge.65).AND.(ylat(i,j).le.85))) ) mask_cal(i,j,2)=.true. |
---|
| 51 | !! -- Scandinavian Ice Sheet Pauline |
---|
| 52 | IF( ((xlong(i,j).gt.0.).AND.(xlong(i,j).le.40)).AND. & |
---|
| 53 | ((ylat(i,j).ge.50).AND.(ylat(i,j).le.70))) mask_cal(i,j,3)=.true. |
---|
| 54 | ! -- Barents Ice Sheet Pauline |
---|
| 55 | IF ( (((xlong(i,j).gt.10.).AND.(xlong(i,j).le.60.)).AND. & |
---|
| 56 | ((ylat(i,j).gt.70).AND.(ylat(i,j).le.85.))).OR. & |
---|
| 57 | (((xlong(i,j).gt.40.).AND.(xlong(i,j).le.60.)).AND. & |
---|
| 58 | ((ylat(i,j).ge.65.).AND.(ylat(i,j).lt.70.)))) mask_cal(i,j,4)=.true. |
---|
| 59 | ! -- Kara Ice Sheet Pauline |
---|
| 60 | IF( (((xlong(i,j).gt.60.).AND.(xlong(i,j).le.122)).AND. & |
---|
| 61 | ((ylat(i,j).ge.70).AND.(ylat(i,j).le.85))).OR. & |
---|
| 62 | (((xlong(i,j).gt.60.).AND.(xlong(i,j).le.80.)).AND. & |
---|
| 63 | ((ylat(i,j).ge.65.).AND.(ylat(i,j).lt.70.)))) mask_cal(i,j,5)=.true. |
---|
| 64 | ! -- British Ice Sheet Pauline |
---|
| 65 | IF( ((xlong(i,j).gt.350.).AND.(xlong(i,j).le.360)).AND. & |
---|
| 66 | ((ylat(i,j).ge.50).AND.(ylat(i,j).le.60))) mask_cal(i,j,6)=.true. |
---|
| 67 | ! -- Siberian Ice Sheet |
---|
| 68 | IF( ((xlong(i,j).gt.122.).AND.(xlong(i,j).lt.190)).AND. & |
---|
| 69 | ((ylat(i,j).ge.60).AND.(ylat(i,j).le.85))) mask_cal(i,j,7)=.true. |
---|
| 70 | enddo |
---|
| 71 | enddo |
---|
| 72 | ! open (UNIT=7777,file='mask_hemin40-02.dat') |
---|
| 73 | ! do j=1,ny |
---|
| 74 | ! do i=1,nx |
---|
| 75 | ! if (mask_cal(i,j,2)) then |
---|
| 76 | ! write_mask(i,j)=2 |
---|
| 77 | ! else |
---|
| 78 | ! write_mask(i,j)=0 |
---|
| 79 | ! endif |
---|
| 80 | ! write(7777,*) i,j,write_mask(i,j) |
---|
| 81 | ! enddo |
---|
| 82 | ! enddo |
---|
| 83 | ! close(UNIT=7777) |
---|
| 84 | ! stop |
---|
| 85 | |
---|
| 86 | ! ecriture entete avec les colonnes du fichier short : |
---|
| 87 | WRITE(num_ritz, '(a91)')"! hemin40 : time,isvol,deltavol,inp,isvolf,inf,isacc,isabl,ISABLBORD,ABLATOT,ISCALV,ISBM" |
---|
| 88 | WRITE(num_ritz, '(a90)')"! colonne : 1 2 3 4 5 6 7 8 9 10 11 12" |
---|
| 89 | WRITE(num_ritz,'(a157)')"! calotte : Fennoscandie ScandinavianIS BarentsIS KaraIS SiberianIS BritishIS" |
---|
| 90 | WRITE(num_ritz,'(a153)')"! isvol : 13 26 39 52 65 78" |
---|
| 91 | WRITE(num_ritz,'(a153)')"! inp : 14 27 40 53 66 79" |
---|
| 92 | WRITE(num_ritz,'(a153)')"! isvolf : 15 28 41 54 67 80" |
---|
| 93 | WRITE(num_ritz,'(a153)')"! inf : 16 29 42 55 68 81" |
---|
| 94 | WRITE(num_ritz,'(a153)')"! isacc : 17 30 43 56 69 82" |
---|
| 95 | WRITE(num_ritz,'(a153)')"! isabl : 18 31 44 57 70 83" |
---|
| 96 | WRITE(num_ritz,'(a153)')"! isablbord : 19 32 45 58 71 84" |
---|
| 97 | WRITE(num_ritz,'(a153)')"! ablatot : 20 33 46 59 72 85" |
---|
| 98 | WRITE(num_ritz,'(a153)')"! iscalv : 21 34 47 60 73 86" |
---|
| 99 | WRITE(num_ritz,'(a153)')"! isbm : 22 35 48 61 74 87" |
---|
| 100 | WRITE(num_ritz,'(a153)')"! itjjamean_ : 23 36 49 62 75 88" |
---|
| 101 | WRITE(num_ritz,'(a153)')"! hmean_ : 24 37 50 63 76 89" |
---|
| 102 | WRITE(num_ritz,'(a153)')"! Hmax_ : 25 38 51 64 77 90" |
---|
| 103 | WRITE(num_ritz,*) |
---|
| 104 | end subroutine init_outshort |
---|
| 105 | |
---|
| 106 | !_________________________________________________________________________ |
---|
| 107 | subroutine shortoutput |
---|
| 108 | |
---|
| 109 | ! 1_initialization |
---|
| 110 | !------------------ |
---|
| 111 | integer KK |
---|
| 112 | ! integer inpl, INPG, INPF |
---|
| 113 | ! integer inplab,inpkeew,inpinn,inpcord |
---|
| 114 | ! integer inpscaN,inpbar,inpkara,inpsib,inpbrit |
---|
[336] | 115 | integer inp(nmask) ! Surface posee (nb de noeuds) |
---|
| 116 | integer inf(nmask) ! surface flottante |
---|
[297] | 117 | |
---|
| 118 | !Variables pour sommer |
---|
[336] | 119 | real isvol(nmask),isvolf(nmask) ! volume posé et flottants |
---|
| 120 | REAL ISCALV(nmask),ISACC(nmask),ISBM(nmask),ISABL(nmask) |
---|
| 121 | REAL ISABLBORD(nmask),ABLATOT(nmask),TACC(nmask),TBM(nmask) |
---|
| 122 | REAL ITJJA(nmask) |
---|
[297] | 123 | ! REAL ABLAMEAN |
---|
| 124 | !moyennes utilisées en output |
---|
[336] | 125 | REAL HMAX_(nmask) , HMEAN_(nmask) |
---|
| 126 | REAL BMEAN_(nmask) , ACCMEAN_(nmask) , ABLMEAN_(nmask) , ABLBORDMEAN_(nmask) |
---|
| 127 | REAL ABLATOTMEAN_(nmask) , CALVMEAN_(nmask) , ITJJAMEAN_(nmask) |
---|
[297] | 128 | |
---|
| 129 | |
---|
| 130 | |
---|
| 131 | ! REAL ITJJAMEAN_L, ITJJAMEAN_G, ITJJAMEAN_F ... |
---|
| 132 | ! REAL ITJJAMEAN_LAB,ITJJAMEAN_KEEW,ITJJAMEAN_INN,ITJJAMEAN_CORD |
---|
| 133 | ! REAL ITJJAMEAN_SCAN,ITJJAMEAN_BAR,ITJJAMEAN_KARA,ITJJAMEAN_SIB,ITJJAMEAN_BRIT |
---|
| 134 | |
---|
| 135 | REAL HDOTMEAN_G |
---|
| 136 | ! REAL ABLA(NX,NY) |
---|
| 137 | REAL DELTAVOL |
---|
| 138 | |
---|
| 139 | |
---|
| 140 | ! open(unit=4145,file='reg_output_nord.dat') |
---|
| 141 | |
---|
| 142 | BMELTMEAN=0. |
---|
| 143 | TBMEAN=0. |
---|
| 144 | TBDOTMEAN=0. |
---|
| 145 | VSMEAN=0. |
---|
| 146 | UZSMEAN=0. |
---|
| 147 | UZSDOTMEAN=0. |
---|
| 148 | UZKMEAN=0. |
---|
| 149 | HDOTMEAN=0. |
---|
| 150 | HDOTMEAN_G = 0. |
---|
| 151 | BDOTMEAN=0. |
---|
| 152 | |
---|
[336] | 153 | DO kk = 1,nmask |
---|
[297] | 154 | INP(kk) = 0 |
---|
| 155 | INF(kk) = 0 |
---|
| 156 | ISVOL(kk) = 0. |
---|
| 157 | ISVOLF(kk) = 0. |
---|
| 158 | ISBM(kk) = 0. |
---|
| 159 | ISACC(kk) = 0. |
---|
| 160 | ISCALV(kk) = 0. |
---|
| 161 | ISABL(kk) = 0. |
---|
| 162 | ISABLBORD(kk) = 0. |
---|
| 163 | ABLATOT(kk)= 0. |
---|
| 164 | TACC(kk) = 0. |
---|
| 165 | TBM(kk) = 0. |
---|
| 166 | ITJJA(kk) = 0. |
---|
| 167 | |
---|
| 168 | ! nouveau tof mai 2009 |
---|
| 169 | ! where (mask_cal(:,:,kk).and.H(:,:).gt.2..and.flot(:,:)) ISVOLF(kk) = ISVOLF(kk) + H(I,J) |
---|
| 170 | isvolf(kk) = sum(H(:,:),mask=(mask_cal(:,:,kk).and.(H(:,:).gt.2.).and.flot(:,:))) |
---|
| 171 | INF(kk) = count(mask_cal(:,:,kk).and.(H(:,:).gt.2.).and.flot(:,:)) |
---|
| 172 | |
---|
| 173 | INP(kk) = count(mask_cal(:,:,kk).and.(H(:,:).gt.2.).and..not.flot(:,:)) |
---|
| 174 | isvol(kk) = sum(H(:,:),mask=(mask_cal(:,:,kk).and.(H(:,:).gt.2.).and..not.flot(:,:))) |
---|
| 175 | isacc(kk) = sum(Acc(:,:),mask=(mask_cal(:,:,kk).and.(H(:,:).gt.2.).and..not.flot(:,:))) |
---|
| 176 | isabl(kk) = sum(Abl(:,:),mask=(mask_cal(:,:,kk).and.(H(:,:).gt.2.).and..not.flot(:,:))) |
---|
| 177 | itjja(kk) = sum(Tjuly(:,:),mask=(mask_cal(:,:,kk).and.(H(:,:).gt.2.).and..not.flot(:,:))) |
---|
| 178 | Hmax_(kk) = maxval(H(:,:),mask=(mask_cal(:,:,kk).and.(H(:,:).gt.2.).and..not.flot(:,:))) |
---|
| 179 | Tacc(kk) = sum(Acc(:,:),mask=(mask_cal(:,:,kk).and.(H(:,:).gt.2.))) |
---|
| 180 | isablbord(kk) = sum(ablbord(:,:),mask=(mask_cal(:,:,kk))) |
---|
| 181 | iscalv(kk) = sum(calv(:,:),mask=(mask_cal(:,:,kk).and.(H(:,:).gt.2.)))/dt |
---|
| 182 | |
---|
| 183 | ablatot(kk) = isabl(kk) + isablbord(kk) |
---|
| 184 | isbm(kk) = isacc(kk)+isabl(kk)+iscalv(kk)+isablbord(kk) |
---|
| 185 | Tbm(kk) = Tacc(kk)+isabl(kk)+iscalv(kk)+isablbord(kk) |
---|
| 186 | |
---|
| 187 | |
---|
| 188 | enddo |
---|
| 189 | |
---|
| 190 | |
---|
| 191 | !!$ |
---|
| 192 | !!$ |
---|
| 193 | !!$ IF(H(i,j).gt.2.) THEN |
---|
| 194 | !!$ if (mk(i,j).eq.1) then |
---|
| 195 | !!$ INF(2) = INF(2) + 1 |
---|
| 196 | !!$ ISVOLF(2) = ISVOLF(2) + H(I,J) |
---|
| 197 | !!$ else |
---|
| 198 | !!$ INP(2) = INP(2) + 1 |
---|
| 199 | !!$ ISVOL(2) = ISVOL(2) + H(I,J) |
---|
| 200 | !!$ ISACC(2) = ISACC(2) + ACC(I,J) |
---|
| 201 | !!$ ISABL(2) = ISABL(2) + ABL(I,J) |
---|
| 202 | !!$ ITJJA(2) = ITJJA(2) + TJULY(I,J) |
---|
| 203 | !!$ if (H(I,J).gt.HMAX_(2)) HMAX_(2)=H(I,J) |
---|
| 204 | !!$ endif |
---|
| 205 | !!$ ENDIF |
---|
| 206 | !!$ TACC(2) = TACC(2) + ACC(I,J) |
---|
| 207 | !!$ ISABLBORD(2) = ISABLBORD(2) + ABLBORD(I,J) |
---|
| 208 | !!$ ABLATOT(2) = ISABL(2) + ISABLBORD(2) |
---|
| 209 | !!$ ISCALV(2) = ISCALV(2) + CALV(I,J) |
---|
| 210 | !!$ ISBM(2) = ISACC(2)+ISABL(2)+ISCALV(2)+ISABLBORD(2) |
---|
| 211 | !!$ TBM(2) = TACC(2)+ISABL(2)+ISCALV(2)+ISABLBORD(2) |
---|
| 212 | !!$ ENDIF |
---|
| 213 | !!$ |
---|
| 214 | !!$ |
---|
| 215 | !!$ END DO |
---|
| 216 | |
---|
| 217 | |
---|
| 218 | |
---|
[336] | 219 | do K=1,nmask |
---|
[297] | 220 | |
---|
| 221 | ! == Les moyennes |
---|
| 222 | IF(INP(K).ne.0)THEN |
---|
| 223 | HMEAN_(K) = ISVOL(K) /INP(K) ! /(DX*DY*INP(K)) |
---|
| 224 | BMEAN_(K) = ISBM(K) /INP(K) |
---|
| 225 | ACCMEAN_(K) = ISACC(K) /INP(K) |
---|
| 226 | ABLMEAN_(K) = ISABL(K) /INP(K) |
---|
| 227 | ABLBORDMEAN_(K) = ISABLBORD(K) /INP(K) |
---|
| 228 | CALVMEAN_(K) = ISCALV(K) /INP(K) |
---|
| 229 | ABLATOTMEAN_(K) = ABLATOT(K) /INP(K) |
---|
| 230 | ITJJAMEAN_(K) = ITJJA(K) /INP(K) |
---|
| 231 | ENDIF |
---|
| 232 | ! == Les volmes intergrées (3D) |
---|
| 233 | ISVOL(K) = ISVOL(K)*DX*DY |
---|
| 234 | isacc(k)=isacc(k)*dx*dy |
---|
| 235 | isabl(k)=isabl(k)*dx*dy |
---|
| 236 | isablbord(k)=isablbord(k)*dx*dy |
---|
| 237 | ablatot(k)=ablatot(k)*dx*dy |
---|
| 238 | iscalv(k)=iscalv(k)*dx*dy |
---|
| 239 | isbm(k)=isbm(k)*dx*dy |
---|
| 240 | tacc(k)=tacc(k)*dx*dy |
---|
| 241 | tbm(k)=tbm(k)*dx*dy |
---|
| 242 | enddo |
---|
| 243 | !nom_table='abl' |
---|
| 244 | !call printtable_r(abl,nom_table) |
---|
| 245 | |
---|
| 246 | !==================================================== |
---|
| 247 | ! -- Hemin40 : : 1 |
---|
| 248 | |
---|
| 249 | if (NP.ne.0) then |
---|
| 250 | ! HMEAN=VOL/NP |
---|
| 251 | ! VOL=VOL*DX*DY |
---|
| 252 | ! ISVOL(1) = VOL |
---|
| 253 | ! BMEAN=BMEAN/NP |
---|
| 254 | ! ACCMEAN=ACCMEAN/NP |
---|
| 255 | ! ITJJAMEAN = ITJJA(1)/NP |
---|
| 256 | ! ABLMEAN=BMEAN-ACCMEAN |
---|
| 257 | !!! ABLMEAN = ABLMEAN/NP |
---|
| 258 | ! CALVMEAN=CALVMEAN/NP |
---|
| 259 | !!! BMELTMEAN=BMELTMEAN/NP |
---|
| 260 | !!! ABLBORDMEAN=ABLBORDMEAN/NP |
---|
| 261 | !!! ABLAMEAN = ABLAMEAN/NP |
---|
| 262 | TBMEAN=TBMEAN/NP |
---|
| 263 | TBDOTMEAN=TBDOTMEAN/NP |
---|
| 264 | VSMEAN=VSMEAN/NP |
---|
| 265 | ! VSDOTMEAN=VSDOTMEAN/NP |
---|
| 266 | UZSMEAN=UZSMEAN/NP |
---|
| 267 | UZSDOTMEAN=UZSDOTMEAN/NP |
---|
| 268 | UZKMEAN=UZKMEAN/NP |
---|
| 269 | DELTAVOL = HDOTMEAN*DX*DY |
---|
| 270 | HDOTMEAN=HDOTMEAN/NP |
---|
| 271 | ELSE |
---|
| 272 | TBMEAN=0. |
---|
| 273 | TBDOTMEAN=0. |
---|
| 274 | VSMEAN=0. |
---|
| 275 | UZSMEAN=0. |
---|
| 276 | UZSDOTMEAN=0. |
---|
| 277 | UZKMEAN=0. |
---|
| 278 | DELTAVOL = 0. |
---|
| 279 | HDOTMEAN=0. |
---|
| 280 | endif |
---|
| 281 | |
---|
| 282 | BDOTMEAN=BDOTMEAN/NX/NY |
---|
| 283 | |
---|
| 284 | ! 2_writing outputs |
---|
| 285 | !------------------ |
---|
| 286 | |
---|
| 287 | ! WRITE(35,904)TIME,VOL,DELTAVOL,NP,ISACC(1),ISABL(1), & |
---|
| 288 | ! ISABLBORD(1),ABLATOT(1),ISCALV(1),ISBM(1), & |
---|
| 289 | ! ISVOL(2),INPL,ISACC(2), & |
---|
| 290 | ! ISABL(2),ISABLBORD(2),ABLATOT(2),ISCALV(2),ISBM(2), & |
---|
| 291 | ! ISVOL(3),INPG,ISACC(3),ISABL(3),ISABLBORD(3), & |
---|
| 292 | ! ABLATOT(3),ISCALV(3),ISBM(3),ISVOL(4),INPF,ISACC(4), & |
---|
| 293 | ! ISABL(4),ISABLBORD(4),ABLATOT(4),ISCALV(4),ISBM(4), & |
---|
| 294 | ! nint(HMEAN),nint(HMAX),BMEAN,TBMEAN,nint(VSMEAN), & |
---|
| 295 | ! TBDOTMEAN,HDOTMEAN,BDOTMEAN,BMELTMEAN,NPAB, & |
---|
| 296 | !! TBDOTMEAN,VSDOTMEAN,HDOTMEAN,BDOTMEAN,BMELTMEAN,NPAB, & |
---|
| 297 | ! NPCAL,HDOTMEAN_G,TACC(1),TBM(1),TACC(2),TBM(2), & |
---|
| 298 | ! TACC(3),TBM(3),TACC(4),TBM(4), & |
---|
| 299 | ! ITJJAMEAN,ITJJAMEAN_L,ITJJAMEAN_G,ITJJAMEAN_F, & |
---|
| 300 | ! ISVOL(5),INPLAB,ISACC(5),ABLATOT(5),ITJJAMEAN_LAB, & |
---|
| 301 | ! ISVOL(6),INPKEEW,ISACC(6),ABLATOT(6),ITJJAMEAN_KEEW, & |
---|
| 302 | ! ISVOL(7),INPINN,ISACC(7),ABLATOT(7),ITJJAMEAN_INN, & |
---|
| 303 | ! ISVOL(8),INPCORD,ISACC(8),ABLATOT(8),ITJJAMEAN_CORD, & |
---|
| 304 | ! ISVOL(9),INPSCAN,ISACC(9),ABLATOT(9),ITJJAMEAN_SCAN, & |
---|
| 305 | ! ISVOL(10),INPBAR,ISACC(10),ABLATOT(10),ITJJAMEAN_BAR, & |
---|
| 306 | ! ISVOL(11),INPKARA,ISACC(11),ABLATOT(11),ITJJAMEAN_KARA, & |
---|
| 307 | ! ISVOL(12),INPSIB,ISACC(12),ABLATOT(12),ITJJAMEAN_SIB, & |
---|
| 308 | ! ISVOL(13),INPBRIT,ISACC(13),ABLATOT(13),ITJJAMEAN_BRIT, & |
---|
| 309 | ! HMEAN,HMEAN_L,HMEAN_G,HMEAN_F,HMEAN_LAB,HMEAN_KEEW, & |
---|
| 310 | ! HMEAN_INN,HMEAN_CORD,HMEAN_SCAN,HMEAN_BAR,HMEAN_KARA, & |
---|
| 311 | ! HMEAN_SIB,HMEAN_BRIT, & |
---|
| 312 | ! HMAX,HMAX_L,HMAX_G,HMAX_F,HMAX_LAB,HMAX_KEEW,HMAX_INN, & |
---|
| 313 | ! HMAX_CORD, HMAX_SCAN,HMAX_BAR,HMAX_KARA,HMAX_SIB,HMAX_BRIT |
---|
| 314 | |
---|
| 315 | |
---|
| 316 | |
---|
| 317 | |
---|
| 318 | !print*,'=========================' |
---|
| 319 | !print*,'ecrit short, time = ',time,dt |
---|
| 320 | !print*,'hdot ISABLBORD',hdot(30:40,100),ISABLBORD(2) |
---|
| 321 | |
---|
| 322 | write(num_ritz,905) time,isvol(1),deltavol,inp(1),isvolf(1),inf(1),isacc(1),isabl(1), & ! 8 |
---|
| 323 | ISABLBORD(1),ABLATOT(1),ISCALV(1),ISBM(1), & ! 4 |
---|
| 324 | isvol(2),inp(2),isvolf(2),inf(2),isacc(2),isabl(2),isablbord(2),ablatot(2),iscalv(2),isbm(2),itjjamean_(2),hmean_(2),Hmax_(2), & ! 13 |
---|
| 325 | isvol(3),inp(3),isvolf(3),inf(3),isacc(3),isabl(3),isablbord(3),ablatot(3),iscalv(3),isbm(3),itjjamean_(3),hmean_(3),hmax_(3), & ! 13 |
---|
| 326 | isvol(4),inp(4),isvolf(4),inf(4),isacc(4),isabl(4),isablbord(4),ablatot(4),iscalv(4),isbm(4),itjjamean_(4),hmean_(4),hmax_(4), & ! 13 |
---|
| 327 | isvol(5),inp(5),isvolf(5),inf(5),isacc(5),isabl(5),isablbord(5),ablatot(5),iscalv(5),isbm(5),itjjamean_(5),hmean_(5),hmax_(5), & ! 13 |
---|
| 328 | isvol(6),inp(6),isvolf(6),inf(6),isacc(6),isabl(6),isablbord(6),ablatot(6),iscalv(6),isbm(6),itjjamean_(6),hmean_(6),hmax_(6), & ! 13 |
---|
[336] | 329 | isvol(7),inp(7),isvolf(7),inf(7),isacc(7),isabl(7),isablbord(7),ablatot(7),iscalv(7),isbm(7),itjjamean_(7),hmean_(7),hmax_(7), & |
---|
[297] | 330 | diff_H, water_bilan, sum(calv_dtt(:,:))/dtt, sum(ablbord_dtt(:,:))/dtt, sum(Bm_dtt(:,:))/dtt, sum(bmelt_dtt(:,:))/dtt !6 |
---|
| 331 | |
---|
| 332 | |
---|
| 333 | ! !! format 900 faux, a reprendre |
---|
| 334 | |
---|
| 335 | 905 format(f9.1,1x, e11.4,1x,e11.5, 1x,i5, 1x,e10.4,1x,i5 , 6(1x,e12.5), & |
---|
| 336 | 12(2(1x,e10.4,1x,i5), 9(1x,e11.4) ) , 6(1x,e15.8)) |
---|
| 337 | |
---|
| 338 | |
---|
| 339 | 940 format('%%%% ',a,' time=',f8.0,' %%%%') |
---|
| 340 | end subroutine shortoutput |
---|
| 341 | |
---|
| 342 | |
---|
| 343 | end module output_eurasie_mod |
---|