Changeset 2816 for branches/2011/dev_r2802_UKMO8_sbccpl
- Timestamp:
- 2011-07-29T10:31:06+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2802_UKMO8_sbccpl/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r2813 r2816 522 522 ! ! ------------------------- ! 523 523 ssnd(jps_fice)%clname = 'OIceFrc' 524 ssnd(jps_hice)%clname = 'OIceTck' 524 525 ssnd(jps_hsnw)%clname = 'OSnwTck' 525 ssnd(jps_hice)%clname = 'OIceTck'526 526 IF( k_ice /= 0 ) THEN 527 527 ssnd(jps_fice)%laction = .TRUE. ! if ice treated in the ocean (even in climato case) … … 532 532 SELECT CASE ( TRIM( sn_snd_thick%cldes ) ) 533 533 CASE ( 'ice and snow' ) 534 ssnd(jps_h snw:jps_hice)%laction = .TRUE.535 ssnd(jps_h snw:jps_hice)%nct = jpl534 ssnd(jps_hice:jps_hsnw)%laction = .TRUE. 535 ssnd(jps_hice:jps_hsnw)%nct = jpl 536 536 CASE ( 'weighted ice and snow' ) 537 ssnd(jps_h snw:jps_hice)%laction = .TRUE.538 IF ( TRIM( sn_snd_thick%clcat ) == 'yes' ) ssnd(jps_h snw:jps_hice)%nct = jpl537 ssnd(jps_hice:jps_hsnw)%laction = .TRUE. 538 IF ( TRIM( sn_snd_thick%clcat ) == 'yes' ) ssnd(jps_hice:jps_hsnw)%nct = jpl 539 539 CASE default ; CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_thick%cldes' ) 540 540 END SELECT … … 1333 1333 IF( ssnd(jps_fice)%laction ) CALL cpl_prism_snd( jps_fice, isec, ztmp3, info ) 1334 1334 1335 ! Send snow and icethickness field1335 ! Send ice and snow thickness field 1336 1336 SELECT CASE( sn_snd_thick%cldes) 1337 1337 CASE( 'weighted ice and snow' ) 1338 1338 SELECT CASE( sn_snd_thick%clcat ) 1339 1339 CASE( 'yes' ) 1340 ztmp3(:,:,1:jpl) = ht_ s(:,:,1:jpl) * a_i(:,:,1:jpl)1341 ztmp4(:,:,1:jpl) = ht_ i(:,:,1:jpl) * a_i(:,:,1:jpl)1340 ztmp3(:,:,1:jpl) = ht_i(:,:,1:jpl) * a_i(:,:,1:jpl) 1341 ztmp4(:,:,1:jpl) = ht_s(:,:,1:jpl) * a_i(:,:,1:jpl) 1342 1342 CASE( 'no' ) 1343 1343 ztmp3(:,:,:) = 0.0 ; ztmp4(:,:,:) = 0.0 1344 1344 DO jl=1,jpl 1345 ztmp3(:,:,1) = ztmp3(:,:,1) + ht_ s(:,:,jl) * a_i(:,:,jl)1346 ztmp4(:,:,1) = ztmp4(:,:,1) + ht_ i(:,:,jl) * a_i(:,:,jl)1345 ztmp3(:,:,1) = ztmp3(:,:,1) + ht_i(:,:,jl) * a_i(:,:,jl) 1346 ztmp4(:,:,1) = ztmp4(:,:,1) + ht_s(:,:,jl) * a_i(:,:,jl) 1347 1347 ENDDO 1348 1348 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_thick%clcat' ) 1349 1349 END SELECT 1350 1350 CASE( 'ice and snow' ) 1351 ztmp3(:,:,1:jpl) = ht_ s(:,:,1:jpl)1352 ztmp4(:,:,1:jpl) = ht_ i(:,:,1:jpl)1351 ztmp3(:,:,1:jpl) = ht_i(:,:,1:jpl) 1352 ztmp4(:,:,1:jpl) = ht_s(:,:,1:jpl) 1353 1353 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_thick%cldes' ) 1354 1354 END SELECT 1355 IF( ssnd(jps_h snw)%laction ) CALL cpl_prism_snd( jps_hsnw, isec, ztmp3, info )1356 IF( ssnd(jps_h ice)%laction ) CALL cpl_prism_snd( jps_hice, isec, ztmp4, info )1355 IF( ssnd(jps_hice)%laction ) CALL cpl_prism_snd( jps_hice, isec, ztmp3, info ) 1356 IF( ssnd(jps_hsnw)%laction ) CALL cpl_prism_snd( jps_hsnw, isec, ztmp4, info ) 1357 1357 ! 1358 1358 #if defined key_cpl_carbon_cycle
Note: See TracChangeset
for help on using the changeset viewer.