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 14908 – NEMO

Changeset 14908


Ignore:
Timestamp:
2021-05-26T17:43:13+02:00 (3 years ago)
Author:
mathiot
Message:

fix consistency between cav and par cases and sign issues

Location:
NEMO/branches/2021/ticket2669_isf_fluxes/src/OCE/ISF
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/ticket2669_isf_fluxes/src/OCE/ISF/isfpar.F90

    r14907 r14908  
    8888      ! 
    8989      ! set temperature content 
    90       ptsc(:,:,jp_tem) = zqh(:,:) * r1_rho0_rcp 
     90      ptsc(:,:,jp_tem) = - zqh(:,:) * r1_rho0_rcp 
    9191      ! 
    9292      ! write restart variables (qoceisf, qhcisf, fwfisf for now and before) 

  • NEMO/branches/2021/ticket2669_isf_fluxes/src/OCE/ISF/isfparmlt.F90

    r14907 r14908  
    112112      CALL isf_tbl(Kmm, ztfrz3d, ztfrz, 'T', misfkt_par, rhisf_tbl_par, misfkb_par, rfrac_tbl_par ) 
    113113      ! 
    114       pqfwf(:,:) = - sf_isfpar_fwf(1)%fnow(:,:,1)      ! fresh water flux from the isf (fwfisf <0 mean melting)  
    115       pqoce(:,:) =   pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux 
    116       pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux  
     114      pqfwf(:,:) = - sf_isfpar_fwf(1)%fnow(:,:,1)      ! fresh water flux from the isf (fwfisf <0 mean melting)       ( > 0 out ) 
     115      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux  ( > 0 out ) 
     116      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux                                            ( > 0 out ) 
    117117      ! 
    118118      CALL iom_put('isftfrz_par', ztfrz ) 
    … …  
    157157      ! 
    158158      ! 2. ------------Net heat flux and fresh water flux due to the ice shelf 
    159       pqoce(:,:) =   rho0 * rcp * rn_gammat0 * risfLeff(:,:) * e1t(:,:) * ( ztavg(:,:) - ztfrz(:,:) ) * r1_e1e2t(:,:) 
    160       pqfwf(:,:) = - pqoce(:,:) / rLfusisf             ! derived from the latent heat flux 
    161       pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux  
     159      pqfwf(:,:) = - rho0 * rcp * rn_gammat0 * risfLeff(:,:) * e1t(:,:) * (ztavg(:,:) - ztfrz(:,:) ) * r1_e1e2t(:,:) / rLfusisf  ! ( > 0 out ) 
     160      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux  ( > 0 out ) 
     161      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux                                            ( > 0 out ) 
    162162      ! 
    163163      ! 3. ------------BG03 output 
Note: See TracChangeset for help on using the changeset viewer.