Changeset 623 for trunk/NEMO/OPA_SRC/DIA/diafwb.F90
- Timestamp:
- 2007-02-26T18:26:08+01:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DIA/diafwb.F90
r474 r623 67 67 68 68 !! * Local declarations 69 INTEGER :: inum ! temporary logical unit 69 70 INTEGER :: ji, jj, jk, jt ! dummy loop indices 70 71 INTEGER :: ii0, ii1, ij0, ij1 … … 387 388 IF ( kt == nitend ) THEN 388 389 389 OPEN(111,FILE='STRAIT.dat') 390 WRITE(111,*) 391 WRITE(111,*) 'Net freshwater budget ' 392 WRITE(111,9010) ' emp = ',a_emp, ' m3 =', a_emp /(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 393 WRITE(111,9010) ' precip = ',a_precip,' m3 =', a_precip/(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 394 WRITE(111,9010) ' a_rnf = ',a_rnf, ' m3 =', a_rnf /(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 395 WRITE(111,*) 396 WRITE(111,9010) ' zarea =',zarea 397 WRITE(111,9010) ' zvol =',zvol 398 WRITE(111,*) 399 WRITE(111,*) 'Mean sea level : ' 400 WRITE(111,9010) ' at nit000 = ',a_sshb ,' m3 ' 401 WRITE(111,9010) ' at nitend = ',a_sshn ,' m3 ' 402 WRITE(111,9010) ' diff = ',(a_sshn-a_sshb),' m3 =', (a_sshn-a_sshb)/(FLOAT(nitend-nit000+1)*rdt) * 1.e-6,' Sv' 403 WRITE(111,9020) ' mean sea level elevation =', a_sshn/zarea,' m' 404 WRITE(111,*) 405 WRITE(111,*) 'Anomaly of salinity content : ' 406 WRITE(111,9010) ' at nit000 = ',a_salb ,' psu.m3 ' 407 WRITE(111,9010) ' at nitend = ',a_saln ,' psu.m3 ' 408 WRITE(111,9010) ' diff = ',(a_saln-a_salb),' psu.m3' 409 WRITE(111,*) 410 WRITE(111,*) 'Mean salinity : ' 411 WRITE(111,9020) ' at nit000 =',a_salb/zvol+zsm0 ,' psu ' 412 WRITE(111,9020) ' at nitend =',a_saln/zvol+zsm0 ,' psu ' 413 WRITE(111,9020) ' diff =',(a_saln-a_salb)/zvol,' psu' 414 WRITE(111,9020) ' S-SLevitus=',a_saln/zvol,' psu' 415 WRITE(111,*) 416 WRITE(111,*) 'Coeff : ' 417 WRITE(111,9030) ' Alpha+ = ', a_aplus 418 WRITE(111,9030) ' Alpha- = ', a_aminus 419 WRITE(111,*) 420 WRITE(111,*) 421 WRITE(111,*) 'Gibraltar : ' 422 WRITE(111,9030) ' Flux entrant (Sv) :', a_flxi(1) 423 WRITE(111,9030) ' Flux sortant (Sv) :', a_flxo(1) 424 WRITE(111,9030) ' T entrant (deg) :', a_temi(1) 425 WRITE(111,9030) ' T sortant (deg) :', a_temo(1) 426 WRITE(111,9030) ' S entrant (psu) :', a_sali(1) 427 WRITE(111,9030) ' S sortant (psu) :', a_salo(1) 428 WRITE(111,*) 429 WRITE(111,*) 'Cadiz : ' 430 WRITE(111,9030) ' Flux entrant (Sv) :', a_flxi(2) 431 WRITE(111,9030) ' Flux sortant (Sv) :', a_flxo(2) 432 WRITE(111,9030) ' T entrant (deg) :', a_temi(2) 433 WRITE(111,9030) ' T sortant (deg) :', a_temo(2) 434 WRITE(111,9030) ' S entrant (psu) :', a_sali(2) 435 WRITE(111,9030) ' S sortant (psu) :', a_salo(2) 436 WRITE(111,*) 437 WRITE(111,*) 'Bab el Mandeb : ' 438 WRITE(111,9030) ' Flux entrant (Sv) :', a_flxi(3) 439 WRITE(111,9030) ' Flux sortant (Sv) :', a_flxo(3) 440 WRITE(111,9030) ' T entrant (deg) :', a_temi(3) 441 WRITE(111,9030) ' T sortant (deg) :', a_temo(3) 442 WRITE(111,9030) ' S entrant (psu) :', a_sali(3) 443 WRITE(111,9030) ' S sortant (psu) :', a_salo(3) 444 WRITE(111,*) 445 WRITE(111,*) 'Baltic : ' 446 WRITE(111,9030) ' Flux entrant (Sv) :', a_flxi(4) 447 WRITE(111,9030) ' Flux sortant (Sv) :', a_flxo(4) 448 WRITE(111,9030) ' T entrant (deg) :', a_temi(4) 449 WRITE(111,9030) ' T sortant (deg) :', a_temo(4) 450 WRITE(111,9030) ' S entrant (psu) :', a_sali(4) 451 WRITE(111,9030) ' S sortant (psu) :', a_salo(4) 452 CLOSE(111) 390 CALL ctlopn( inum, 'STAIT.dat', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 391 & 1, numout, .TRUE., 1 ) 392 WRITE(inum,*) 393 WRITE(inum,*) 'Net freshwater budget ' 394 WRITE(inum,9010) ' emp = ',a_emp, ' m3 =', a_emp /(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 395 WRITE(inum,9010) ' precip = ',a_precip,' m3 =', a_precip/(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 396 WRITE(inum,9010) ' a_rnf = ',a_rnf, ' m3 =', a_rnf /(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 397 WRITE(inum,*) 398 WRITE(inum,9010) ' zarea =',zarea 399 WRITE(inum,9010) ' zvol =',zvol 400 WRITE(inum,*) 401 WRITE(inum,*) 'Mean sea level : ' 402 WRITE(inum,9010) ' at nit000 = ',a_sshb ,' m3 ' 403 WRITE(inum,9010) ' at nitend = ',a_sshn ,' m3 ' 404 WRITE(inum,9010) ' diff = ',(a_sshn-a_sshb),' m3 =', (a_sshn-a_sshb)/(FLOAT(nitend-nit000+1)*rdt) * 1.e-6,' Sv' 405 WRITE(inum,9020) ' mean sea level elevation =', a_sshn/zarea,' m' 406 WRITE(inum,*) 407 WRITE(inum,*) 'Anomaly of salinity content : ' 408 WRITE(inum,9010) ' at nit000 = ',a_salb ,' psu.m3 ' 409 WRITE(inum,9010) ' at nitend = ',a_saln ,' psu.m3 ' 410 WRITE(inum,9010) ' diff = ',(a_saln-a_salb),' psu.m3' 411 WRITE(inum,*) 412 WRITE(inum,*) 'Mean salinity : ' 413 WRITE(inum,9020) ' at nit000 =',a_salb/zvol+zsm0 ,' psu ' 414 WRITE(inum,9020) ' at nitend =',a_saln/zvol+zsm0 ,' psu ' 415 WRITE(inum,9020) ' diff =',(a_saln-a_salb)/zvol,' psu' 416 WRITE(inum,9020) ' S-SLevitus=',a_saln/zvol,' psu' 417 WRITE(inum,*) 418 WRITE(inum,*) 'Coeff : ' 419 WRITE(inum,9030) ' Alpha+ = ', a_aplus 420 WRITE(inum,9030) ' Alpha- = ', a_aminus 421 WRITE(inum,*) 422 WRITE(inum,*) 423 WRITE(inum,*) 'Gibraltar : ' 424 WRITE(inum,9030) ' Flux entrant (Sv) :', a_flxi(1) 425 WRITE(inum,9030) ' Flux sortant (Sv) :', a_flxo(1) 426 WRITE(inum,9030) ' T entrant (deg) :', a_temi(1) 427 WRITE(inum,9030) ' T sortant (deg) :', a_temo(1) 428 WRITE(inum,9030) ' S entrant (psu) :', a_sali(1) 429 WRITE(inum,9030) ' S sortant (psu) :', a_salo(1) 430 WRITE(inum,*) 431 WRITE(inum,*) 'Cadiz : ' 432 WRITE(inum,9030) ' Flux entrant (Sv) :', a_flxi(2) 433 WRITE(inum,9030) ' Flux sortant (Sv) :', a_flxo(2) 434 WRITE(inum,9030) ' T entrant (deg) :', a_temi(2) 435 WRITE(inum,9030) ' T sortant (deg) :', a_temo(2) 436 WRITE(inum,9030) ' S entrant (psu) :', a_sali(2) 437 WRITE(inum,9030) ' S sortant (psu) :', a_salo(2) 438 WRITE(inum,*) 439 WRITE(inum,*) 'Bab el Mandeb : ' 440 WRITE(inum,9030) ' Flux entrant (Sv) :', a_flxi(3) 441 WRITE(inum,9030) ' Flux sortant (Sv) :', a_flxo(3) 442 WRITE(inum,9030) ' T entrant (deg) :', a_temi(3) 443 WRITE(inum,9030) ' T sortant (deg) :', a_temo(3) 444 WRITE(inum,9030) ' S entrant (psu) :', a_sali(3) 445 WRITE(inum,9030) ' S sortant (psu) :', a_salo(3) 446 WRITE(inum,*) 447 WRITE(inum,*) 'Baltic : ' 448 WRITE(inum,9030) ' Flux entrant (Sv) :', a_flxi(4) 449 WRITE(inum,9030) ' Flux sortant (Sv) :', a_flxo(4) 450 WRITE(inum,9030) ' T entrant (deg) :', a_temi(4) 451 WRITE(inum,9030) ' T sortant (deg) :', a_temo(4) 452 WRITE(inum,9030) ' S entrant (psu) :', a_sali(4) 453 WRITE(inum,9030) ' S sortant (psu) :', a_salo(4) 454 CLOSE(inum) 453 455 ENDIF 454 456
Note: See TracChangeset
for help on using the changeset viewer.