Changeset 9675 for branches/UKMO/dev_merge_2017_CICE_interface/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Timestamp:
- 2018-05-29T15:47:05+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_merge_2017_CICE_interface/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r9631 r9675 367 367 IF (ksbc == jp_flx) THEN 368 368 DO jl=1,ncat 369 ztmpn(:,:,jl)= qla_ice(:,:,1)*a_i(:,:,jl)369 ztmpn(:,:,jl)=zevap_ice(:,:,1)*a_i(:,:,jl) 370 370 ENDDO 371 371 ELSE IF (ksbc == jp_purecpl) THEN 372 372 DO jl=1,ncat 373 ztmpn(:,:,jl)= qla_ice(:,:,jl)*a_i(:,:,jl)373 ztmpn(:,:,jl)=zevap_ice(:,:,jl)*a_i(:,:,jl) * (-1.0*Lsub) 374 374 ENDDO 375 375 ELSE 376 !In coupled mode - qla_ice calculated in sbc_cpl for each category377 ztmpn(:,:,1:ncat)= qla_ice(:,:,1:ncat)376 !In coupled mode - zevap_ice calculated in sbc_cpl for each category 377 ztmpn(:,:,1:ncat)=zevap_ice(:,:,1:ncat) 378 378 ENDIF 379 379 … … 384 384 ! Convert to GBM 385 385 IF (ksbc == jp_flx .OR. ksbc == jp_purecpl) THEN 386 ztmp(:,:) = botmelt(:,:,jl)*a_i(:,:,jl)386 ztmp(:,:) = qcn_ice(:,:,jl)*a_i(:,:,jl) 387 387 ELSE 388 ztmp(:,:) = botmelt(:,:,jl)388 ztmp(:,:) = qcn_ice(:,:,jl) 389 389 ENDIF 390 390 CALL nemo2cice(ztmp,fcondtopn_f(:,:,jl,:),'T', 1. ) … … 393 393 ! Convert to GBM 394 394 IF (ksbc == jp_flx .OR. ksbc == jp_purecpl) THEN 395 ztmp(:,:) = ( topmelt(:,:,jl)+botmelt(:,:,jl))*a_i(:,:,jl)395 ztmp(:,:) = (qml_ice(:,:,jl)+qcn_ice(:,:,jl))*a_i(:,:,jl) 396 396 ELSE 397 ztmp(:,:) = ( topmelt(:,:,jl)+botmelt(:,:,jl))397 ztmp(:,:) = (qml_ice(:,:,jl)+qcn_ice(:,:,jl)) 398 398 ENDIF 399 399 CALL nemo2cice(ztmp,fsurfn_f(:,:,jl,:),'T', 1. ) … … 738 738 ! Meltpond fraction and depth 739 739 DO jl = 1,ncat 740 CALL cice2nemo(apeffn(:,:,jl,:),a_ p(:,:,jl),'T', 1. )741 CALL cice2nemo(trcrn(:,:,nt_hpnd,jl,:), ht_p(:,:,jl),'T', 1. )740 CALL cice2nemo(apeffn(:,:,jl,:),a_ip(:,:,jl),'T', 1. ) 741 CALL cice2nemo(trcrn(:,:,nt_hpnd,jl,:),v_ip(:,:,jl),'T', 1. ) 742 742 ENDDO 743 743 #endif … … 752 752 IF (heat_capacity) THEN 753 753 DO jl = 1,ncat 754 CALL cice2nemo(Tn_top(:,:,jl,:),t n_ice(:,:,jl),'T', 1. )755 CALL cice2nemo(keffn_top(:,:,jl,:), kn_ice(:,:,jl),'T', 1. )754 CALL cice2nemo(Tn_top(:,:,jl,:),t1_ice(:,:,jl),'T', 1. ) 755 CALL cice2nemo(keffn_top(:,:,jl,:),cnd_ice(:,:,jl),'T', 1. ) 756 756 ENDDO 757 757 ! Convert surface temperature to Kelvin 758 t n_ice(:,:,:)=tn_ice(:,:,:)+rt0758 t1_ice(:,:,:)=t1_ice(:,:,:)+rt0 759 759 ELSE 760 t n_ice(:,:,:) = 0.0761 kn_ice(:,:,:) = 0.0760 t1_ice(:,:,:) = 0.0 761 cnd_ice(:,:,:) = 0.0 762 762 ENDIF 763 763 … … 877 877 tprecip(:,:) = sf(jp_snow)%fnow(:,:,1)+sf(jp_rain)%fnow(:,:,1) 878 878 ! May be better to do this conversion somewhere else 879 qla_ice(:,:,1) = -Lsub*sf(jp_sblm)%fnow(:,:,1)880 topmelt(:,:,1) = sf(jp_top1)%fnow(:,:,1)881 topmelt(:,:,2) = sf(jp_top2)%fnow(:,:,1)882 topmelt(:,:,3) = sf(jp_top3)%fnow(:,:,1)883 topmelt(:,:,4) = sf(jp_top4)%fnow(:,:,1)884 topmelt(:,:,5) = sf(jp_top5)%fnow(:,:,1)885 botmelt(:,:,1) = sf(jp_bot1)%fnow(:,:,1)886 botmelt(:,:,2) = sf(jp_bot2)%fnow(:,:,1)887 botmelt(:,:,3) = sf(jp_bot3)%fnow(:,:,1)888 botmelt(:,:,4) = sf(jp_bot4)%fnow(:,:,1)889 botmelt(:,:,5) = sf(jp_bot5)%fnow(:,:,1)879 zevap_ice(:,:,1) = -Lsub*sf(jp_sblm)%fnow(:,:,1) 880 qml_ice(:,:,1) = sf(jp_top1)%fnow(:,:,1) 881 qml_ice(:,:,2) = sf(jp_top2)%fnow(:,:,1) 882 qml_ice(:,:,3) = sf(jp_top3)%fnow(:,:,1) 883 qml_ice(:,:,4) = sf(jp_top4)%fnow(:,:,1) 884 qml_ice(:,:,5) = sf(jp_top5)%fnow(:,:,1) 885 qcn_ice(:,:,1) = sf(jp_bot1)%fnow(:,:,1) 886 qcn_ice(:,:,2) = sf(jp_bot2)%fnow(:,:,1) 887 qcn_ice(:,:,3) = sf(jp_bot3)%fnow(:,:,1) 888 qcn_ice(:,:,4) = sf(jp_bot4)%fnow(:,:,1) 889 qcn_ice(:,:,5) = sf(jp_bot5)%fnow(:,:,1) 890 890 891 891 ! control print (if less than 100 time-step asked)
Note: See TracChangeset
for help on using the changeset viewer.