Changeset 3977


Ignore:
Timestamp:
2013-07-16T14:04:24+02:00 (7 years ago)
Author:
flavoni
Message:

add print of fwb correction value, and hminrhg for LIM2, in CNRS LIM3 branch

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  
    5151   REAL(wp), PUBLIC ::   ahi0   = 350._wp      !: sea-ice hor. eddy diffusivity coeff. (m2/s) 
    5252   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 
    5354 
    5455   REAL(wp), PUBLIC ::   usecc2                !:  = 1.0 / ( ecc * ecc ) 
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/iceini.F90

    r3938 r3977  
    8585      CALL lim_itd_ini                 ! ice thickness distribution initialization 
    8686      ! 
    87       CALL lim_sbc_init                ! ice surface boundary condition    
    88  
    8987 
    9088      !                                ! Initial sea-ice state 
     
    102100      ENDIF 
    103101      ! 
     102      CALL lim_sbc_init                ! ice surface boundary condition    
     103      ! 
    104104      fr_i(:,:) = at_i(:,:)           ! initialisation of sea-ice fraction 
    105105      ! 
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

    r3964 r3977  
    5454   REAL(wp)  ::   rzero  = 0._wp     
    5555   REAL(wp)  ::   rone   = 1._wp 
     56   REAL(wp) ::   r1_rau0   ! = 1 / rau0 
    5657 
    5758   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   utau_oce, vtau_oce   ! air-ocean surface i- & j-stress     [N/m2] 
     
    409410      !!------------------------------------------------------------------- 
    410411      REAL(wp) :: zsum, zarea 
     412      r1_rau0 = 1._wp / rau0 
    411413      ! 
    412414      IF(lwp) WRITE(numout,*) 
     
    436438         snwice_mass  (:,:) = tms(:,:) * ( rhosn * vt_s(:,:) + rhoic * vt_i(:,:)  ) 
    437439         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 
    443446         ENDIF 
    444447      ! 
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90

    r3964 r3977  
    142142            !                                         ! estimate from the previous year budget 
    143143         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  
    145147         !  
    146148         IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN      ! correct the freshwater fluxes 
Note: See TracChangeset for help on using the changeset viewer.