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 9675 for branches/UKMO/dev_merge_2017_CICE_interface/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90 – NEMO

Ignore:
Timestamp:
2018-05-29T15:47:05+02:00 (6 years ago)
Author:
dancopsey
Message:

Apply changes as suggested by Alex West - renaming coupling variables.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_merge_2017_CICE_interface/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r9631 r9675  
    367367         IF (ksbc == jp_flx) THEN 
    368368            DO jl=1,ncat 
    369                ztmpn(:,:,jl)=qla_ice(:,:,1)*a_i(:,:,jl) 
     369               ztmpn(:,:,jl)=zevap_ice(:,:,1)*a_i(:,:,jl) 
    370370            ENDDO 
    371371         ELSE IF (ksbc == jp_purecpl) THEN 
    372372            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) 
    374374            ENDDO 
    375375    ELSE 
    376            !In coupled mode - qla_ice calculated in sbc_cpl for each category 
    377            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) 
    378378         ENDIF 
    379379 
     
    384384!  Convert to GBM 
    385385            IF (ksbc == jp_flx .OR. ksbc == jp_purecpl) THEN 
    386                ztmp(:,:) = botmelt(:,:,jl)*a_i(:,:,jl) 
     386               ztmp(:,:) = qcn_ice(:,:,jl)*a_i(:,:,jl) 
    387387            ELSE 
    388                ztmp(:,:) = botmelt(:,:,jl) 
     388               ztmp(:,:) = qcn_ice(:,:,jl) 
    389389            ENDIF 
    390390            CALL nemo2cice(ztmp,fcondtopn_f(:,:,jl,:),'T', 1. ) 
     
    393393!  Convert to GBM 
    394394            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)  
    396396            ELSE 
    397                ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl)) 
     397               ztmp(:,:) = (qml_ice(:,:,jl)+qcn_ice(:,:,jl)) 
    398398            ENDIF 
    399399            CALL nemo2cice(ztmp,fsurfn_f(:,:,jl,:),'T', 1. ) 
     
    738738! Meltpond fraction and depth 
    739739      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. ) 
    742742      ENDDO 
    743743#endif 
     
    752752      IF (heat_capacity) THEN 
    753753         DO jl = 1,ncat 
    754             CALL cice2nemo(Tn_top(:,:,jl,:),tn_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. ) 
    756756         ENDDO 
    757757! Convert surface temperature to Kelvin 
    758          tn_ice(:,:,:)=tn_ice(:,:,:)+rt0 
     758         t1_ice(:,:,:)=t1_ice(:,:,:)+rt0 
    759759      ELSE 
    760          tn_ice(:,:,:) = 0.0 
    761          kn_ice(:,:,:) = 0.0 
     760         t1_ice(:,:,:) = 0.0 
     761         cnd_ice(:,:,:) = 0.0 
    762762      ENDIF        
    763763 
     
    877877      tprecip(:,:) = sf(jp_snow)%fnow(:,:,1)+sf(jp_rain)%fnow(:,:,1) 
    878878! 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) 
    890890 
    891891      ! control print (if less than 100 time-step asked) 
Note: See TracChangeset for help on using the changeset viewer.