Changeset 4292 for branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_2
- Timestamp:
- 2013-11-20T17:28:04+01:00 (10 years ago)
- Location:
- branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_2
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90
r4148 r4292 24 24 USE phycst ! physical constants 25 25 USE dom_oce ! ocean domain 26 USE domvvl ! ocean vertical scale factors 26 27 USE dom_ice_2 ! LIM-2: ice domain 27 28 USE ice_2 ! LIM-2: ice variables … … 59 60 !! * Substitutions 60 61 # include "vectopt_loop_substitute.h90" 62 # include "domzgr_substitute.h90" 61 63 !!---------------------------------------------------------------------- 62 64 !! NEMO/LIM2 4.0 , UCL - NEMO Consortium (2011) … … 446 448 !!------------------------------------------------------------------- 447 449 ! 450 INTEGER :: jk ! local integer 451 ! 448 452 IF(lwp) WRITE(numout,*) 449 453 IF(lwp) WRITE(numout,*) 'lim_sbc_init_2 : LIM-2 sea-ice - surface boundary condition' … … 472 476 snwice_mass (:,:) = 0.e0 ! no mass exchanges 473 477 snwice_mass_b(:,:) = 0.e0 ! no mass exchanges 478 snwice_fmass (:,:) = 0.e0 ! no mass exchanges 474 479 ENDIF 475 480 IF( nn_ice_embd == 2 .AND. & ! full embedment (case 2) & no restart : … … 477 482 sshn(:,:) = sshn(:,:) - snwice_mass(:,:) * r1_rau0 478 483 sshb(:,:) = sshb(:,:) - snwice_mass(:,:) * r1_rau0 484 do jk = 1,jpkm1 ! adjust initial vertical scale factors 485 fse3t_n(:,:,jk) = e3t_0(:,:,jk)*( 1._wp + sshn(:,:)*tmask(:,:,1)/(ht_0(:,:) + 1.0 - tmask(:,:,1)) ) 486 fse3t_b(:,:,jk) = e3t_0(:,:,jk)*( 1._wp + sshb(:,:)*tmask(:,:,1)/(ht_0(:,:) + 1.0 - tmask(:,:,1)) ) 487 end do 488 fse3t_a(:,:,:) = fse3t_b(:,:,:) 489 ! Reconstruction of all vertical scale factors at now and before time steps 490 ! ============================================================================= 491 ! Horizontal scale factor interpolations 492 ! -------------------------------------- 493 CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3u_b(:,:,:), 'U' ) 494 CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3v_b(:,:,:), 'V' ) 495 CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3u_n(:,:,:), 'U' ) 496 CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3v_n(:,:,:), 'V' ) 497 CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n(:,:,:), 'F' ) 498 ! Vertical scale factor interpolations 499 ! ------------------------------------ 500 CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n (:,:,:), 'W' ) 501 CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3uw_n(:,:,:), 'UW' ) 502 CALL dom_vvl_interpol( fse3v_n(:,:,:), fse3vw_n(:,:,:), 'VW' ) 503 CALL dom_vvl_interpol( fse3u_b(:,:,:), fse3uw_b(:,:,:), 'UW' ) 504 CALL dom_vvl_interpol( fse3v_b(:,:,:), fse3vw_b(:,:,:), 'VW' ) 505 ! t- and w- points depth 506 ! ---------------------- 507 fsdept_n(:,:,1) = 0.5_wp * fse3w_n(:,:,1) 508 fsdepw_n(:,:,1) = 0.0_wp 509 fsde3w_n(:,:,1) = fsdept_n(:,:,1) - sshn(:,:) 510 DO jk = 2, jpk 511 fsdept_n(:,:,jk) = fsdept_n(:,:,jk-1) + fse3w_n(:,:,jk) 512 fsdepw_n(:,:,jk) = fsdepw_n(:,:,jk-1) + fse3t_n(:,:,jk-1) 513 fsde3w_n(:,:,jk) = fsdept_n(:,:,jk ) - sshn (:,:) 514 END DO 479 515 ENDIF 480 516 ! -
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90
r4147 r4292 237 237 238 238 ! energy needed to bring ocean surface layer until its freezing 239 qcmif (ji,jj) = rau0 * rcp * fse3t_m(ji,jj ,1) * ( tfu(ji,jj) - sst_m(ji,jj) - rt0 ) * ( 1 - zinda )239 qcmif (ji,jj) = rau0 * rcp * fse3t_m(ji,jj) * ( tfu(ji,jj) - sst_m(ji,jj) - rt0 ) * ( 1 - zinda ) 240 240 241 241 ! calculate oceanic heat flux.
Note: See TracChangeset
for help on using the changeset viewer.