Changeset 6711


Ignore:
Timestamp:
2016-06-14T15:24:49+02:00 (4 years ago)
Author:
clem
Message:

correct a bug in coupled mode in case calving is activated

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r6489 r6711  
    14211421      CALL lim_thd_snwblow( p_frld, zsnw ) 
    14221422       
    1423       ! --- evaporation (kg/m2/s) --- ! 
     1423      ! --- evaporation (used later in sbccpl) --- ! 
     1424      zevap(:,:) = zemp_tot(:,:) + ztprecip(:,:) 
     1425 
     1426      ! --- evaporation over ice (kg/m2/s) --- ! 
    14241427      zevap_ice(:,:) = frcv(jpr_ievp)%z3(:,:,1) 
    14251428      ! since the sensitivity of evap to temperature (devap/dT) is not prescribed by the atmosphere, we set it to 0 
     
    14291432      ! --- evaporation minus precipitation corrected for the effect of wind blowing on snow --- ! 
    14301433      zemp_oce(:,:) = zemp_tot(:,:) - zemp_ice(:,:) - zsprecip * (1._wp - zsnw) 
    1431       zemp_ice(:,:) = zemp_ice(:,:) + zsprecip * (1._wp - zsnw)           
    1432  
    1433       ! Sublimation over sea-ice (cell average) 
    1434       IF( iom_use('subl_ai_cea') )  CALL iom_put( 'subl_ai_cea', zevap_ice(:,:) * zicefr(:,:) ) 
    1435       ! runoffs and calving (put in emp_tot) 
     1434      zemp_ice(:,:) = zemp_ice(:,:) + zsprecip * (1._wp - zsnw)     
     1435 
     1436      ! --- runoffs (included in emp later on) --- ! 
    14361437      IF( srcv(jpr_rnf)%laction )   rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1) 
     1438 
     1439      ! --- calving (put in emp_tot and emp_oce) --- ! 
    14371440      IF( srcv(jpr_cal)%laction ) THEN  
    14381441         zemp_tot(:,:) = zemp_tot(:,:) - frcv(jpr_cal)%z3(:,:,1) 
     1442         zemp_oce(:,:) = zemp_oce(:,:) - frcv(jpr_cal)%z3(:,:,1) 
    14391443         CALL iom_put( 'calving_cea', frcv(jpr_cal)%z3(:,:,1) ) 
    14401444      ENDIF 
    1441  
     1445       
    14421446      IF( ln_mixcpl ) THEN 
    14431447         emp_tot(:,:) = emp_tot(:,:) * xcplmask(:,:,0) + zemp_tot(:,:) * zmsk(:,:) 
     
    14611465         ENDDO 
    14621466      ENDIF 
    1463  
     1467       
     1468      IF( iom_use('subl_ai_cea') )   CALL iom_put( 'subl_ai_cea', zevap_ice(:,:) * zicefr(:,:)    )  ! Sublimation over sea-ice (cell average) 
    14641469                                     CALL iom_put( 'snowpre'    , sprecip                         )  ! Snow 
    14651470      IF( iom_use('snow_ao_cea') )   CALL iom_put( 'snow_ao_cea', sprecip(:,:) * ( 1._wp - zsnw ) )  ! Snow over ice-free ocean  (cell average) 
     
    15491554 
    15501555#if defined key_lim3       
    1551       ! --- evaporation --- ! 
    1552       zevap(:,:) = zemp_tot(:,:) + ztprecip(:,:) ! evaporation over ocean 
    1553  
    15541556      ! --- non solar flux over ocean --- ! 
    15551557      !         note: p_frld cannot be = 0 since we limit the ice concentration to amax 
Note: See TracChangeset for help on using the changeset viewer.