Changeset 3977 for branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO
- Timestamp:
- 2013-07-16T14:04:24+02:00 (11 years ago)
- Location:
- branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90
r2715 r3977 51 51 REAL(wp), PUBLIC :: ahi0 = 350._wp !: sea-ice hor. eddy diffusivity coeff. (m2/s) 52 52 REAL(wp), PUBLIC :: alphaevp = 1._wp !: coefficient for the solution of EVP int. stresses 53 REAL(wp), PUBLIC :: hminrhg = 0.05_wp !: clem : ice thickness (in m) below which ice velocity is set to ocean velocity 53 54 54 55 REAL(wp), PUBLIC :: usecc2 !: = 1.0 / ( ecc * ecc ) -
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/iceini.F90
r3938 r3977 85 85 CALL lim_itd_ini ! ice thickness distribution initialization 86 86 ! 87 CALL lim_sbc_init ! ice surface boundary condition88 89 87 90 88 ! ! Initial sea-ice state … … 102 100 ENDIF 103 101 ! 102 CALL lim_sbc_init ! ice surface boundary condition 103 ! 104 104 fr_i(:,:) = at_i(:,:) ! initialisation of sea-ice fraction 105 105 ! -
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90
r3964 r3977 54 54 REAL(wp) :: rzero = 0._wp 55 55 REAL(wp) :: rone = 1._wp 56 REAL(wp) :: r1_rau0 ! = 1 / rau0 56 57 57 58 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: utau_oce, vtau_oce ! air-ocean surface i- & j-stress [N/m2] … … 409 410 !!------------------------------------------------------------------- 410 411 REAL(wp) :: zsum, zarea 412 r1_rau0 = 1._wp / rau0 411 413 ! 412 414 IF(lwp) WRITE(numout,*) … … 436 438 snwice_mass (:,:) = tms(:,:) * ( rhosn * vt_s(:,:) + rhoic * vt_i(:,:) ) 437 439 snwice_mass_b(:,:) = snwice_mass(:,:) 438 IF( .NOT.ln_rstart ) THEN ! delete the initial ssh below sea-ice area 439 zarea = glob_sum( e1e2t(:,:) ) ! interior global domain surface 440 zsum = glob_sum( e1e2t(:,:) * ( snwice_mass(:,:) ) ) / zarea / rau0 441 sshn(:,:) = sshn(:,:) - zsum 442 sshb(:,:) = sshb(:,:) - zsum 440 IF( .NOT. ln_rstart ) THEN ! delete the initial ssh below sea-ice area 441 ! 442 zarea = glob_sum( e1e2t(:,:) ) ! interior global domain surface 443 zsum = glob_sum( e1e2t(:,:) * ( snwice_mass(:,:) ) ) / zarea * r1_rau0 444 sshn(:,:) = sshn(:,:) - zsum 445 sshb(:,:) = sshb(:,:) - zsum 443 446 ENDIF 444 447 ! -
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90
r3964 r3977 142 142 ! ! estimate from the previous year budget 143 143 ENDIF 144 IF(lwp) WRITE(numout,*) 'correction of fresh water budget :', fwfold 144 IF( MOD( kt, 5475 ) == 0 ) THEN 145 IF(lwp) WRITE(numout,*) 'correction of fresh water budget :', fwfold 146 ENDIF 145 147 ! 146 148 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN ! correct the freshwater fluxes
Note: See TracChangeset
for help on using the changeset viewer.