- Timestamp:
- 2010-10-12T16:04:22+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/SBC/sbccpl.F90
r1970 r2224 40 40 USE restart ! 41 41 USE oce , ONLY : tn, un, vn 42 USE phycst, ONLY : rt0, rcp43 42 USE albedo ! 44 43 USE in_out_manager ! I/O manager … … 153 152 INTEGER , DIMENSION( jprcv) :: nrcvinfo ! OASIS info argument 154 153 154 #if ! defined key_lim2 && ! defined key_lim3 155 ! quick patch to be able to run the coupled model without sea-ice... 156 INTEGER, PARAMETER :: jpl = 1 157 REAL(wp), DIMENSION(jpi,jpj ) :: hicif, hsnif, u_ice, v_ice,fr1_i0,fr2_i0 158 REAL(wp), DIMENSION(jpi,jpj,jpl) :: tn_ice, alb_ice 159 REAL(wp) :: lfus 160 #endif 161 155 162 !! Substitution 156 163 # include "vectopt_loop_substitute.h90" … … 254 261 255 262 ! default definitions of srcv 256 srcv(:)%laction = .FALSE. ; srcv(:)%clgrid = 'T' ; srcv(:)%nsgn = 1 263 srcv(:)%laction = .FALSE. ; srcv(:)%clgrid = 'T' ; srcv(:)%nsgn = 1. 257 264 258 265 ! ! ------------------------- ! … … 439 446 440 447 ! default definitions of nsnd 441 ssnd(:)%laction = .FALSE. ; ssnd(:)%clgrid = 'T' ; ssnd(:)%nsgn = 1 448 ssnd(:)%laction = .FALSE. ; ssnd(:)%clgrid = 'T' ; ssnd(:)%nsgn = 1. 442 449 443 450 ! ! ------------------------- ! … … 492 499 ssnd(jps_ocz1)%clname = 'O_OCurz1' ; ssnd(jps_ivz1)%clname = 'O_IVelz1' 493 500 ! 494 ssnd(jps_ocx1:jps_ivz1)%nsgn = -1 501 ssnd(jps_ocx1:jps_ivz1)%nsgn = -1. ! vectors: change of the sign at the north fold 495 502 496 503 IF( cn_snd_crt(4) /= 'T' ) CALL ctl_stop( 'cn_snd_crt(4) must be equal to T' ) … … 714 721 ! ! non solar heat flux over the ocean (qns) 715 722 IF( srcv(jpr_qnsoce)%laction ) qns(:,:) = frcv(:,:,jpr_qnsoce) 716 IF( srcv(jpr_qnsmix)%laction ) qns(:,:) = frcv(:,:,jpr_qnsmix) 717 qns(:,:) = qns(:,:) - frcv(:,:,jpr_snow) * lfus ! add the latent heat of solid precip. melting 723 IF( srcv(jpr_qnsmix)%laction ) qns(:,:) = frcv(:,:,jpr_qnsmix) 724 ! add the latent heat of solid precip. melting 725 IF( srcv(jpr_snow )%laction ) qns(:,:) = qns(:,:) - frcv(:,:,jpr_snow) * lfus 718 726 719 727 ! ! solar flux over the ocean (qsr)
Note: See TracChangeset
for help on using the changeset viewer.