Changeset 9427


Ignore:
Timestamp:
2018-03-26T15:16:45+02:00 (2 years ago)
Author:
clem
Message:

small bug in ice thickness redistribution

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/iceitd.F90

    r9421 r9427  
    449449               !           
    450450               !                                                     ! Ice age  
    451                ztrans          = oa_i_2d(ji,jl1) * pdaice(ji,jl)  !!clem: should be * zworka(ji) but it does not work ???? 
     451               ztrans          = oa_i_2d(ji,jl1) * zworka(ji) 
    452452               oa_i_2d(ji,jl1) = oa_i_2d(ji,jl1) - ztrans 
    453453               oa_i_2d(ji,jl2) = oa_i_2d(ji,jl2) + ztrans 
    454454               ! 
    455455               ztrans          = sv_i_2d(ji,jl1) * zworkv(ji)        ! Ice salinity 
    456                ! 
    457456               sv_i_2d(ji,jl1) = sv_i_2d(ji,jl1) - ztrans 
    458457               sv_i_2d(ji,jl2) = sv_i_2d(ji,jl2) + ztrans 
    459458               ! 
    460459               !                                                     ! Surface temperature 
    461                ztrans          = t_su_2d(ji,jl1) * pdaice(ji,jl)     !!clem: should be * zworka(ji) but it does not work ???? 
     460               ztrans          = zaTsfn(ji,jl1) * zworka(ji) 
    462461               zaTsfn(ji,jl1)  = zaTsfn(ji,jl1) - ztrans 
    463462               zaTsfn(ji,jl2)  = zaTsfn(ji,jl2) + ztrans 
     
    465464               IF ( ln_pnd_H12 ) THEN 
    466465                  !                                                  ! Pond fraction 
    467                   ztrans          = a_ip_2d(ji,jl1) * pdaice(ji,jl) !!clem: should be * zworka(ji) but it does not work 
     466                  ztrans          = a_ip_2d(ji,jl1) * zworka(ji) 
    468467                  a_ip_2d(ji,jl1) = a_ip_2d(ji,jl1) - ztrans 
    469468                  a_ip_2d(ji,jl2) = a_ip_2d(ji,jl2) + ztrans 
    470469                  !                                                  ! Pond volume (also proportional to da/a) 
    471                   ztrans          = v_ip_2d(ji,jl1) * pdaice(ji,jl) !!clem: should be * zworka(ji) but it does not work 
     470                  ztrans          = v_ip_2d(ji,jl1) * zworka(ji) 
    472471                  v_ip_2d(ji,jl1) = v_ip_2d(ji,jl1) - ztrans 
    473472                  v_ip_2d(ji,jl2) = v_ip_2d(ji,jl2) + ztrans 
Note: See TracChangeset for help on using the changeset viewer.