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 9944 for NEMO/trunk/tests/ICEDYN/MY_SRC/usrdef_sbc.F90 – NEMO

Ignore:
Timestamp:
2018-07-13T16:33:14+02:00 (6 years ago)
Author:
mathiot
Message:

wrong variable name in ICEDYN userdef

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/ICEDYN/MY_SRC/usrdef_sbc.F90

    r9789 r9944  
    132132      emp_oce  (:,:)   = emp_oce(:,:) - sprecip(:,:) * (1._wp - zsnw(:,:) ) 
    133133      qevap_ice(:,:,:) =   0._wp 
    134       qprec_ice(:,:)   =   rhosn * ( sst_m(:,:) * cpic - lfus ) * tmask(:,:,1) !  in J/m3 
     134      qprec_ice(:,:)   =   rhos * ( sst_m(:,:) * rcpi - rLfus ) * tmask(:,:,1) !  in J/m3 
    135135      qemp_oce (:,:)   = - emp_oce(:,:) * sst_m(:,:) * rcp 
    136       qemp_ice (:,:)   =   sprecip(:,:) * zsnw * ( sst_m(:,:) * cpic - lfus ) * tmask(:,:,1) ! solid precip (only) 
     136      qemp_ice (:,:)   =   sprecip(:,:) * zsnw * ( sst_m(:,:) * rcpi - rLfus ) * tmask(:,:,1) ! solid precip (only) 
    137137 
    138138      ! total fluxes 
     
    146146      ! 
    147147      WHERE    ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) <  0.1_wp )       ! linear decrease from hi=0 to 10cm   
    148          qsr_ice_tr(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 
     148         qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 
    149149      ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp )       ! constant (zfr1) when hi>10cm 
    150          qsr_ice_tr(:,:,:) = qsr_ice(:,:,:) * zfr1 
     150         qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * zfr1 
    151151      ELSEWHERE                                                         ! zero when hs>0 
    152          qsr_ice_tr(:,:,:) = 0._wp  
     152         qtr_ice_top(:,:,:) = 0._wp  
    153153      END WHERE 
    154154           
Note: See TracChangeset for help on using the changeset viewer.