Changeset 3964


Ignore:
Timestamp:
2013-07-09T17:42:43+02:00 (7 years ago)
Author:
flavoni
Message:

correction of fresh water budget for LIM3, see ticket 1116

Location:
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

    r3962 r3964  
    286286      !   mass of snow and ice per unit area          ! 
    287287      !-----------------------------------------------! 
    288       IF( nn_ice_embd /= 0 ) THEN                               ! embedded sea-ice (mass required) 
     288     !SF IF( nn_ice_embd /= 0 ) THEN                               ! embedded sea-ice (mass required) 
    289289         snwice_mass_b(:,:) = snwice_mass(:,:)                  ! save mass from the previous ice time step 
    290290         !                                                      ! new mass per unit area 
     
    292292         !                                                      ! time evolution of snow+ice mass 
    293293         snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) * r1_rdtice 
    294       ENDIF 
     294      !SF ENDIF 
    295295 
    296296      !-----------------------------------------------! 
     
    408408      !! ** input   : Namelist namicedia 
    409409      !!------------------------------------------------------------------- 
     410      REAL(wp) :: zsum, zarea 
    410411      ! 
    411412      IF(lwp) WRITE(numout,*) 
     
    435436         snwice_mass  (:,:) = tms(:,:) * ( rhosn * vt_s(:,:) + rhoic * vt_i(:,:)  ) 
    436437         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 
     443         ENDIF 
    437444      ! 
    438445      ! 
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90

    r3962 r3964  
    3636   REAL(wp) ::   fwfold    ! fwfold to be suppressed 
    3737   REAL(wp) ::   area      ! global mean ocean surface (interior domain) 
     38   REAL(wp) ::   r1_rau0   ! = 1 / rau0  
    3839 
    3940   !! * Substitutions 
     
    130131         ikty = 365 * 86400 / rdttra(1)    !!bug  use of 365 days leap year or 360d year !!!!!!! 
    131132         IF( MOD( kt, ikty ) == 0 ) THEN 
     133            ! 
     134            r1_rau0 = 1._wp / rau0 
     135            ! 
    132136            a_fwb_b = a_fwb                           ! mean sea level taking into account the ice+snow 
    133137                                                      ! sum over the global domain 
     
    138142            !                                         ! estimate from the previous year budget 
    139143         ENDIF 
     144         IF(lwp) WRITE(numout,*) 'correction of fresh water budget :', fwfold 
    140145         !  
    141146         IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN      ! correct the freshwater fluxes 
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/oce.F90

    r3962 r3964  
    8585         &     spgu  (jpi,jpj)   , spgv(jpi,jpj)   ,                       & 
    8686         &     gtsu(jpi,jpj,jpts), gtsv(jpi,jpj,jpts),                     & 
    87          &     gru(jpi,jpj)      , grv(jpi,jpj)                      , STAT=ierr(2) ) 
     87         &     gru(jpi,jpj)      , grv(jpi,jpj),                           & 
     88         &     snwice_mass(jpi,jpj), snwice_mass_b(jpi,jpj), snwice_fmass(jpi,jpj)    , STAT=ierr(2) ) 
    8889         ! 
    8990      oce_alloc = MAXVAL( ierr ) 
Note: See TracChangeset for help on using the changeset viewer.