Changeset 167 for trunk/NEMO
- Timestamp:
- 2004-11-05T14:39:09+01:00 (20 years ago)
- Location:
- trunk/NEMO/OPA_SRC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DOM/domhgr.F90
r93 r167 309 309 ! ========== 310 310 311 ! angle 45 ° and ze1=106.e+3 / jp_cfg forced -> zlam1 = -85°, zphi1 = 29°N311 ! angle 45deg and ze1=106.e+3 / jp_cfg forced -> zlam1 = -85deg, zphi1 = 29degN 312 312 zlam1 = -85 313 313 zphi1 = 29 314 ze1 = 106000. / FLOAT(jp_cfg) ! resolution in meters 315 IF( nbench /= 0 ) ze1 = 106000.e0 ! benchmark: forced the resolution to be about 100 km 314 ! resolution in meters 315 ze1 = 106000. / FLOAT(jp_cfg) 316 ! benchmark: forced the resolution to be about 100 km 317 IF( nbench /= 0 ) ze1 = 106000.e0 316 318 zsin_alpha = - SQRT( 2. ) / 2. 317 319 zcos_alpha = SQRT( 2. ) / 2. 318 320 ze1deg = ze1 / (ra * rad) 319 321 IF( nbench /= 0 ) ze1deg = ze1deg / FLOAT(jp_cfg) ! benchmark: keep the lat/+lon 320 ! 322 ! ! at the right jp_cfg resolution 321 323 glam0 = zlam1 + zcos_alpha * ze1deg * FLOAT( jpjglo-2 ) 322 324 gphi0 = zphi1 + zsin_alpha * ze1deg * FLOAT( jpjglo-2 ) -
trunk/NEMO/OPA_SRC/SBC/ocesbc.F90
r160 r167 655 655 ztimemax = zdaymax * 24 656 656 ! 21th december day in hours 657 ztimemin = ztimemax + rjjhh * raajj / 2 ! rjjhh * raajj / 4 = 1 seasonal cycle in hours 657 ! rjjhh * raajj / 4 = 1 seasonal cycle in hours 658 ztimemin = ztimemax + rjjhh * raajj / 2 658 659 ! amplitudes 659 660 zemp_S = 0.7 ! intensity of COS in the South 660 661 zemp_N = 0.8 ! intensity of COS in the North 661 662 zemp_sais= 0.1 662 zTstar = 28.3 ! intemsity from 28.3 à -5° 663 zcos_sais = COS( (ztime - ztimemax) / (ztimemin - ztimemax) * rpi ) ! 1/2 period between 21th June and 21th December 663 zTstar = 28.3 ! intemsity from 28.3 a -5 deg 664 ! 1/2 period between 21th June and 21th December 665 zcos_sais = COS( (ztime - ztimemax) / (ztimemin - ztimemax) * rpi ) 664 666 ztrp= - 40. ! retroaction term (W/m2/K) 665 667 zconv = 3.16e-5 ! convert 1m/yr->3.16e-5mm/s 666 668 DO jj = 1, jpj 667 DO ji = 1, jpi 668 t_star (ji,jj) = zTstar * ( 1 + 1. / 50. * zcos_sais ) & ! domain from 15° to 50° 669 & * COS( rpi * (gphit(ji,jj) - 5.) & ! between 27 and 28 °C at 15N, -3 670 & / (53.5 * ( 1 + 11 / 53.5 * zcos_sais ) * 2.) ) ! and 13°C at 50N 53.5 + or - 11 671 ! ! = 1/4 period : 672 ! ! 64.5 in summer, 42.5 in winter 673 qt (ji,jj) = ztrp * ( tb(ji,jj,1) - t_star(ji,jj) ) 674 IF( gphit(ji,jj) >= 14.845 .AND. 37.2 >= gphit(ji,jj)) THEN 675 emp (ji,jj) = zemp_S * zconv & 676 & * SIN( rpi / 2 * (gphit(ji,jj) - 37.2) / (24.6 - 37.2) ) & ! zero at 37.8°, max at 24.6° 677 & * ( 1 - zemp_sais / zemp_S * zcos_sais) 678 emps (ji,jj) = emp (ji,jj) 679 ELSE 680 emp (ji,jj) = - zemp_N * zconv & 681 & * SIN( rpi / 2 * (gphit(ji,jj) - 37.2) / (46.8 - 37.2) ) & ! zero at 37.8°, max at 46.8° 682 & * ( 1 - zemp_sais / zemp_N * zcos_sais ) 683 emps (ji,jj) = emp (ji,jj) 684 ENDIF 685 qsr (ji,jj) = 230 * COS( 3.1415 * ( gphit(ji,jj) - 23.5 * zcos_sais ) / ( 0.9 * 180 ) ) ! 23.5° : tropics 686 END DO 669 DO ji = 1, jpi 670 ! domain from 15 deg to 50 deg between 27 and 28 degC at 15N, -3 671 ! and 13 degC at 50N 53.5 + or - 11 = 1/4 period : 672 ! 64.5 in summer, 42.5 in winter 673 t_star (ji,jj) = zTstar * ( 1 + 1. / 50. * zcos_sais ) & 674 & * COS( rpi * (gphit(ji,jj) - 5.) & 675 & / (53.5 * ( 1 + 11 / 53.5 * zcos_sais ) * 2.) ) 676 qt (ji,jj) = ztrp * ( tb(ji,jj,1) - t_star(ji,jj) ) 677 IF( gphit(ji,jj) >= 14.845 .AND. 37.2 >= gphit(ji,jj)) THEN 678 ! zero at 37.8 deg, max at 24.6 deg 679 emp (ji,jj) = zemp_S * zconv & 680 & * SIN( rpi / 2 * (gphit(ji,jj) - 37.2) / (24.6 - 37.2) ) & 681 & * ( 1 - zemp_sais / zemp_S * zcos_sais) 682 emps (ji,jj) = emp (ji,jj) 683 ELSE 684 ! zero at 37.8 deg, max at 46.8 deg 685 emp (ji,jj) = - zemp_N * zconv & 686 & * SIN( rpi / 2 * (gphit(ji,jj) - 37.2) / (46.8 - 37.2) ) & 687 & * ( 1 - zemp_sais / zemp_N * zcos_sais ) 688 emps (ji,jj) = emp (ji,jj) 689 ENDIF 690 ! 23.5 deg : tropics 691 qsr (ji,jj) = 230 * COS( 3.1415 * ( gphit(ji,jj) - 23.5 * zcos_sais ) / ( 0.9 * 180 ) ) 692 END DO 687 693 END DO 688 694 ! compute the emp flux such as its integration on the whole domain and at each time be zero -
trunk/NEMO/OPA_SRC/SBC/taumod.F90
r93 r167 169 169 ztimemax = zdaymax * 24 170 170 ! 21th december day in hours 171 ztimemin = ztimemax + rjjhh * raajj / 2 ! rjjhh * raajj / 4 = 1 seasonal cycle in hours 172 173 ztau = 0.105 / SQRT(2.) ! mean intensity at 0.105 ; srqt(2) because projected with 45° angle 174 ztau_sais = 0.015 ! seasonal oscillation intensity 171 ! rjjhh * raajj / 4 = 1 seasonal cycle in hours 172 ztimemin = ztimemax + rjjhh * raajj / 2 173 174 ! mean intensity at 0.105 ; srqt(2) because projected with 45deg angle 175 ztau = 0.105 / SQRT(2.) 176 ! seasonal oscillation intensity 177 ztau_sais = 0.015 175 178 ztaun = ztau - ztau_sais * COS( (ztime - ztimemax) / (ztimemin - ztimemax) * rpi ) 176 179 DO jj = 1, jpj 177 180 DO ji = 1, jpi 178 ! domain from 15° to 50° and 1/2 period along 14° so 5/4 of half period with seasonal cycle 181 ! domain from 15deg to 50deg and 1/2 period along 14deg 182 ! so 5/4 of half period with seasonal cycle 179 183 taux (ji,jj) = - ztaun * SIN( rpi * (gphiu(ji,jj) - 15.) / (29.-15.) ) 180 184 tauy (ji,jj) = ztaun * SIN( rpi * (gphiv(ji,jj) - 15.) / (29.-15.) )
Note: See TracChangeset
for help on using the changeset viewer.