Ignore:
Timestamp:
2020-04-23T10:11:38+02:00 (7 months ago)
Author:
gsamson
Message:

move rhoa computation from sbcabl to sbcblk (to avoid restartability issue over sea-ice); correct prt_ctl calls in ablmod and comment a ctl_stop in iom to allow using 3D atm forcing with bulk mode

Location:
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/OCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/OCE/IOM/iom.F90

    r12489 r12799  
    11931193                           &         'we accept this case, even if there is a possible mix-up between z and time dimension' )    
    11941194                     idmspc = idmspc - 1 
    1195                   ELSE 
    1196                      CALL ctl_stop( TRIM(clinfo), 'To keep iom lisibility, when reading a '//clrankpv//'D array,'         ,   & 
    1197                         &                         'we do not accept data with '//cldmspc//' spatial dimensions',   & 
    1198                         &                         'Use ncwa -a to suppress the unnecessary dimensions' ) 
     1195                  !!GS: possibility to read 3D ABL atmopsheric forcing and use 1st level to force BULK simulation 
     1196                  !ELSE 
     1197                  !   CALL ctl_stop( TRIM(clinfo), 'To keep iom lisibility, when reading a '//clrankpv//'D array,',   & 
     1198                  !      &                         'we do not accept data with '//cldmspc//' spatial dimensions'  ,   & 
     1199                  !      &                         'Use ncwa -a to suppress the unnecessary dimensions' ) 
    11991200                  ENDIF 
    12001201            ENDIF 
  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/OCE/SBC/sbcblk.F90

    r12588 r12799  
    658658 
    659659      IF( ln_abl ) THEN         !==  ABL formulation  ==!   multiplication by rho_air and turbulent fluxes computation done in ablstp 
    660          !! FL do we need this multiplication by tmask ... ??? 
    661660         DO_2D_11_11 
    662             zztmp = zU_zu(ji,jj) !* tmask(ji,jj,1) 
     661            zztmp = zU_zu(ji,jj) 
    663662            wndm(ji,jj)   = zztmp                   ! Store zU_zu in wndm to compute ustar2 in ablmod 
    664663            pcd_du(ji,jj) = zztmp * zcd_oce(ji,jj) 
    665664            psen(ji,jj)   = zztmp * zch_oce(ji,jj) 
    666665            pevp(ji,jj)   = zztmp * zce_oce(ji,jj) 
     666            rhoa(ji,jj)   = rho_air( ptair(ji,jj), phumi(ji,jj), pslp(ji,jj) ) 
    667667         END_2D 
    668668      ELSE                      !==  BLK formulation  ==!   turbulent fluxes computation 
Note: See TracChangeset for help on using the changeset viewer.