- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r5487 r6808 243 243 ! momentum fluxes (utau, vtau ) ! 244 244 !------------------------------------! 245 !CDIR COLLAPSE246 245 utau(:,:) = sf(jp_utau)%fnow(:,:,1) 247 !CDIR COLLAPSE248 246 vtau(:,:) = sf(jp_vtau)%fnow(:,:,1) 249 247 … … 251 249 ! wind stress module (taum ) ! 252 250 !------------------------------------! 253 !CDIR NOVERRCHK254 251 DO jj = 2, jpjm1 255 !CDIR NOVERRCHK256 252 DO ji = fs_2, fs_jpim1 ! vector opt. 257 253 ztx2 = utau(ji-1,jj ) + utau(ji,jj) … … 268 264 ! store the wind speed (wndm ) ! 269 265 !------------------------------------! 270 !CDIR COLLAPSE271 266 wndm(:,:) = sf(jp_wndm)%fnow(:,:,1) 272 267 wndm(:,:) = wndm(:,:) * tmask(:,:,1) … … 281 276 ! Other ocean fluxes ! 282 277 !------------------------! 283 !CDIR NOVERRCHK284 !CDIR COLLAPSE285 278 DO jj = 1, jpj 286 !CDIR NOVERRCHK287 279 DO ji = 1, jpi 288 280 ! … … 375 367 zcprec = rcp / rday ! convert prec ( mm/day ==> m/s) ==> W/m2 376 368 377 !CDIR COLLAPSE378 369 emp(:,:) = zqla(:,:) / cevap & ! freshwater flux 379 370 & - sf(jp_prec)%fnow(:,:,1) / rday * tmask(:,:,1) 380 371 ! 381 !CDIR COLLAPSE382 372 qns(:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:) & ! Downward Non Solar flux 383 373 & - zqla(:,:) * pst(:,:) * zcevap & ! remove evap. heat content at SST in Celcius … … 415 405 416 406 # if defined key_lim2 || defined key_lim3 407 417 408 SUBROUTINE blk_ice_clio_tau 418 409 !!--------------------------------------------------------------------------- … … 429 420 ! 430 421 IF( nn_timing == 1 ) CALL timing_start('blk_ice_clio_tau') 431 422 ! 432 423 SELECT CASE( cp_ice_msh ) 433 424 ! 434 425 CASE( 'C' ) ! C-grid ice dynamics 435 426 ! 436 427 zcoef = cai / cao ! Change from air-sea stress to air-ice stress 437 428 utau_ice(:,:) = zcoef * utau(:,:) 438 429 vtau_ice(:,:) = zcoef * vtau(:,:) 439 430 ! 440 431 CASE( 'I' ) ! I-grid ice dynamics: I-point (i.e. F-point lower-left corner) 441 432 ! 442 433 zcoef = 0.5_wp * cai / cao ! Change from air-sea stress to air-ice stress 443 434 DO jj = 2, jpj ! stress from ocean U- and V-points to ice U,V point … … 447 438 END DO 448 439 END DO 449 440 ! 450 441 CALL lbc_lnk( utau_ice(:,:), 'I', -1. ) ; CALL lbc_lnk( vtau_ice(:,:), 'I', -1. ) ! I-point 451 442 ! 452 443 END SELECT 453 444 ! 454 445 IF(ln_ctl) THEN 455 446 CALL prt_ctl(tab2d_1=utau_ice , clinfo1=' blk_ice_clio: utau_ice : ', tab2d_2=vtau_ice , clinfo2=' vtau_ice : ') 456 447 ENDIF 457 448 ! 458 449 IF( nn_timing == 1 ) CALL timing_stop('blk_ice_clio_tau') 459 450 ! 460 451 END SUBROUTINE blk_ice_clio_tau 452 461 453 #endif 462 454 463 455 # if defined key_lim2 || defined key_lim3 456 464 457 SUBROUTINE blk_ice_clio_flx( ptsu , palb_cs, palb_os, palb ) 465 458 !!--------------------------------------------------------------------------- … … 520 513 !-------------------------------------------------------------------------------- 521 514 522 !CDIR NOVERRCHK523 !CDIR COLLAPSE524 515 DO jj = 1, jpj 525 !CDIR NOVERRCHK526 516 DO ji = 1, jpi 527 517 ztatm (ji,jj) = sf(jp_tair)%fnow(ji,jj,1) ! air temperature in Kelvins … … 573 563 574 564 ! ! ========================== ! 575 DO jl = 1, jpl ! Loop over ice categories !565 DO jl = 1, jpl ! Loop over ice categories ! 576 566 ! ! ========================== ! 577 !CDIR NOVERRCHK578 !CDIR COLLAPSE579 567 DO jj = 1 , jpj 580 !CDIR NOVERRCHK581 568 DO ji = 1, jpi 582 569 !-------------------------------------------! … … 636 623 ! ----------------------------------------------------------------------------- ! 637 624 ! 638 !CDIR COLLAPSE639 625 qns_ice(:,:,:) = z_qlw (:,:,:) - z_qsb (:,:,:) - qla_ice (:,:,:) ! Downward Non Solar flux 640 !CDIR COLLAPSE641 626 tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) / rday ! total precipitation [kg/m2/s] 642 627 ! … … 644 629 ! Correct the OCEAN non solar flux with the existence of solid precipitation ! 645 630 ! ---------------=====--------------------------------------------------------- ! 646 !CDIR COLLAPSE647 631 qns(:,:) = qns(:,:) & ! update the non-solar heat flux with: 648 632 & - sprecip(:,:) * lfus & ! remove melting solid precip … … 782 766 ! Saturated water vapour and vapour pressure 783 767 ! ------------------------------------------ 784 !CDIR NOVERRCHK785 !CDIR COLLAPSE786 768 DO jj = 1, jpj 787 !CDIR NOVERRCHK788 769 DO ji = 1, jpi 789 770 ztamr = sf(jp_tair)%fnow(ji,jj,1) - rtt … … 814 795 zdaycor = 1.0 + 0.0013 * SIN( zdist ) + 0.0342 * COS( zdist ) 815 796 816 !CDIR NOVERRCHK817 797 DO jj = 1, jpj 818 !CDIR NOVERRCHK819 798 DO ji = 1, jpi 820 799 ! product of sine (cosine) of latitude and sine (cosine) of solar declination … … 837 816 838 817 ! compute and sum ocean qsr over the daylight (i.e. between sunrise and sunset) 839 !CDIR NOVERRCHK840 818 DO jt = 1, jp24 841 819 zcoef = FLOAT( jt ) - 0.5 842 !CDIR NOVERRCHK843 !CDIR COLLAPSE844 820 DO jj = 1, jpj 845 !CDIR NOVERRCHK846 821 DO ji = 1, jpi 847 822 zlha = COS( zlsrise(ji,jj) - zcoef * zdlha(ji,jj) ) ! local hour angle … … 862 837 ! Taking into account the ellipsity of the earth orbit, the clouds AND masked if sea-ice cover > 0% 863 838 zcoef1 = srgamma * zdaycor / ( 2. * rpi ) 864 !CDIR COLLAPSE865 839 DO jj = 1, jpj 866 840 DO ji = 1, jpi … … 920 894 ! Saturated water vapour and vapour pressure 921 895 ! ------------------------------------------ 922 !CDIR NOVERRCHK923 !CDIR COLLAPSE924 896 DO jj = 1, jpj 925 !CDIR NOVERRCHK926 897 DO ji = 1, jpi 927 898 ztamr = sf(jp_tair)%fnow(ji,jj,1) - rtt … … 952 923 zdaycor = 1.0 + 0.0013 * SIN( zdist ) + 0.0342 * COS( zdist ) 953 924 954 !CDIR NOVERRCHK955 925 DO jj = 1, jpj 956 !CDIR NOVERRCHK957 926 DO ji = 1, jpi 958 927 ! product of sine (cosine) of latitude and sine (cosine) of solar declination … … 979 948 DO jl = 1, ijpl ! loop over ice categories ! 980 949 ! !----------------------------! 981 !CDIR NOVERRCHK982 950 DO jt = 1, jp24 983 951 zcoef = FLOAT( jt ) - 0.5 984 !CDIR NOVERRCHK985 !CDIR COLLAPSE986 952 DO jj = 1, jpj 987 !CDIR NOVERRCHK988 953 DO ji = 1, jpi 989 954 zlha = COS( zlsrise(ji,jj) - zcoef * zdlha(ji,jj) ) ! local hour angle
Note: See TracChangeset
for help on using the changeset viewer.