New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3905 for trunk/NEMOGCM/NEMO/OPA_SRC – NEMO

Ignore:
Timestamp:
2013-05-24T15:40:39+02:00 (11 years ago)
Author:
cetlod
Message:

bugfix : store the flux due to freezing/melting, needed at least for PISCES, see ticket #1100

Location:
trunk/NEMOGCM/NEMO/OPA_SRC/SBC
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

    r3680 r3905  
    7070   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   sfx    , sfx_b    !: salt flux                                    [PSU/m2/s] 
    7171   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   emp_tot           !: total E-P over ocean and ice                 [Kg/m2/s] 
     72   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   fmmflx            !: freshwater budget: freezing/melting          [Kg/m2/s] 
    7273   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   rnf    , rnf_b    !: river runoff   [Kg/m2/s]   
    7374   !! 
    … …  
    115116         &      qsr_tot(jpi,jpj) , qsr  (jpi,jpj) ,                        & 
    116117         &      emp    (jpi,jpj) , emp_b(jpi,jpj) ,                        & 
    117          &      sfx    (jpi,jpj) , sfx_b(jpi,jpj) , emp_tot(jpi,jpj) , STAT=ierr(2) ) 
     118         &      sfx    (jpi,jpj) , sfx_b(jpi,jpj) , emp_tot(jpi,jpj), fmmflx(jpi,jpj), STAT=ierr(2) ) 
    118119         ! 
    119120      ALLOCATE( rnf  (jpi,jpj) , sbc_tsc  (jpi,jpj,jpts) , qsr_hc  (jpi,jpj,jpk) ,     & 

  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r3764 r3905  
    146146      sfx(:,:) = 0.0_wp                            ! the salt flux due to freezing/melting will be computed (i.e. will be non-zero)  
    147147                                                   ! only if sea-ice is present 
     148  
     149      fmmflx(:,:) = 0.0_wp                        ! freezing-melting array initialisation 
    148150 
    149151      !                                            ! restartability    
    … …  
    362364                                                                ! (includes virtual salt flux beneath ice  
    363365                                                                ! in linear free surface case) 
     366         CALL iom_put( "fmmflx", fmmflx  )                      ! Freezing-melting water flux 
    364367         CALL iom_put( "qt"    , qns  + qsr )                   ! total heat flux  
    365368         CALL iom_put( "qns"   , qns        )                   ! solar heat flux 
Note: See TracChangeset for help on using the changeset viewer.