Changeset 13286 for NEMO/trunk/src/OCE/ISF
- Timestamp:
- 2020-07-09T17:48:29+02:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 2 2 ^/utils/build/makenemo@HEAD makenemo 3 3 ^/utils/build/mk@HEAD mk 4 ^/utils/tools /@HEADtools4 ^/utils/tools@HEAD tools 5 5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM … … 8 8 9 9 # SETTE 10 ^/utils/CI/ sette@12931sette10 ^/utils/CI/r12931_sette_ticket2366@HEAD sette
-
- Property svn:externals
-
NEMO/trunk/src/OCE/ISF/isfcpl.F90
r13237 r13286 20 20 USE domqco , ONLY: dom_qco_zgr ! vertical scale factor interpolation 21 21 #endif 22 USE dom ngb, ONLY: dom_ngb ! find the closest grid point from a given lon/lat position22 USE domutl , ONLY: dom_ngb ! find the closest grid point from a given lon/lat position 23 23 ! 24 24 USE oce ! ocean dynamics and tracers … … 183 183 !!---------------------------------------------------------------------- 184 184 ! 185 CALL iom_get( numror, jpdom_auto glo, 'ssmask' , zssmask_b, ldxios = lrxios ) ! need to extrapolate T/S185 CALL iom_get( numror, jpdom_auto, 'ssmask' , zssmask_b, ldxios = lrxios ) ! need to extrapolate T/S 186 186 187 187 ! compute new ssh if we open a full water column … … 264 264 !!---------------------------------------------------------------------- 265 265 ! 266 CALL iom_get( numror, jpdom_auto glo, 'tmask' , ztmask_b, ldxios = lrxios ) ! need to extrapolate T/S267 !CALL iom_get( numror, jpdom_auto glo, 'wmask' , zwmask_b, ldxios = lrxios ) ! need to extrapolate T/S268 !CALL iom_get( numror, jpdom_auto glo, 'gdepw_n', zdepw_b(:,:,:), ldxios = lrxios ) ! need to interpol vertical profile (vvl)266 CALL iom_get( numror, jpdom_auto, 'tmask' , ztmask_b, ldxios = lrxios ) ! need to extrapolate T/S 267 !CALL iom_get( numror, jpdom_auto, 'wmask' , zwmask_b, ldxios = lrxios ) ! need to extrapolate T/S 268 !CALL iom_get( numror, jpdom_auto, 'gdepw_n', zdepw_b(:,:,:), ldxios = lrxios ) ! need to interpol vertical profile (vvl) 269 269 ! 270 270 ! … … 410 410 !!---------------------------------------------------------------------- 411 411 ! 412 CALL iom_get( numror, jpdom_auto glo, 'tmask' , ztmask_b, ldxios = lrxios )413 CALL iom_get( numror, jpdom_auto glo, 'e3u_n' , ze3u_b , ldxios = lrxios )414 CALL iom_get( numror, jpdom_auto glo, 'e3v_n' , ze3v_b , ldxios = lrxios )412 CALL iom_get( numror, jpdom_auto, 'tmask' , ztmask_b, ldxios = lrxios ) 413 CALL iom_get( numror, jpdom_auto, 'e3u_n' , ze3u_b , ldxios = lrxios ) 414 CALL iom_get( numror, jpdom_auto, 'e3v_n' , ze3v_b , ldxios = lrxios ) 415 415 ! 416 416 ! 1.0: compute horizontal volume flux divergence difference before-after coupling … … 520 520 521 521 ! get restart variable 522 CALL iom_get( numror, jpdom_auto glo, 'tmask' , ztmask_b(:,:,:), ldxios = lrxios ) ! need to extrapolate T/S523 CALL iom_get( numror, jpdom_auto glo, 'e3t_n' , ze3t_b(:,:,:) , ldxios = lrxios )524 CALL iom_get( numror, jpdom_auto glo, 'tn' , zt_b(:,:,:) , ldxios = lrxios )525 CALL iom_get( numror, jpdom_auto glo, 'sn' , zs_b(:,:,:) , ldxios = lrxios )522 CALL iom_get( numror, jpdom_auto, 'tmask' , ztmask_b(:,:,:), ldxios = lrxios ) ! need to extrapolate T/S 523 CALL iom_get( numror, jpdom_auto, 'e3t_n' , ze3t_b(:,:,:) , ldxios = lrxios ) 524 CALL iom_get( numror, jpdom_auto, 'tn' , zt_b(:,:,:) , ldxios = lrxios ) 525 CALL iom_get( numror, jpdom_auto, 'sn' , zs_b(:,:,:) , ldxios = lrxios ) 526 526 527 527 ! compute run length … … 544 544 545 545 DO jk = 1,jpk-1 546 DO jj = nldj,nlej547 DO ji = nldi,nlei546 DO jj = Njs0,Nje0 547 DO ji = Nis0,Nie0 548 548 549 549 ! volume diff … … 578 578 nisfl(:)=0 579 579 DO jk = 1,jpk-1 580 DO jj = nldj,nlej581 DO ji = nldi,nlei580 DO jj = Njs0,Nje0 581 DO ji = Nis0,Nie0 582 582 jip1=MIN(ji+1,jpi) ; jim1=MAX(ji-1,1) ; jjp1=MIN(jj+1,jpj) ; jjm1=MAX(jj-1,1) ; 583 583 IF ( tmask(ji,jj,jk) == 0._wp .AND. ztmask_b(ji,jj,jk) == 1._wp ) THEN … … 600 600 jisf = 0 601 601 DO jk = 1,jpk-1 602 DO jj = nldj,nlej603 DO ji = nldi,nlei602 DO jj = Njs0,Nje0 603 DO ji = Nis0,Nie0 604 604 IF ( tmask(ji,jj,jk) == 0._wp .AND. ztmask_b(ji,jj,jk) == 1._wp ) THEN 605 605 -
NEMO/trunk/src/OCE/ISF/isfrst.F90
r11931 r13286 53 53 IF( iom_varid( numror, cfwf_b, ldstop = .FALSE. ) > 0 ) THEN 54 54 IF(lwp) WRITE(numout,*) ' nit000-1 isf tracer content forcing fields read in the restart file' 55 CALL iom_get( numror, jpdom_auto glo, cfwf_b, pfwf_b(:,:) , ldxios = lrxios ) ! before ice shelf melt56 CALL iom_get( numror, jpdom_auto glo, chc_b , ptsc_b (:,:,jp_tem), ldxios = lrxios ) ! before ice shelf heat flux57 CALL iom_get( numror, jpdom_auto glo, csc_b , ptsc_b (:,:,jp_sal), ldxios = lrxios ) ! before ice shelf heat flux55 CALL iom_get( numror, jpdom_auto, cfwf_b, pfwf_b(:,:) , ldxios = lrxios ) ! before ice shelf melt 56 CALL iom_get( numror, jpdom_auto, chc_b , ptsc_b (:,:,jp_tem), ldxios = lrxios ) ! before ice shelf heat flux 57 CALL iom_get( numror, jpdom_auto, csc_b , ptsc_b (:,:,jp_sal), ldxios = lrxios ) ! before ice shelf heat flux 58 58 ELSE 59 59 pfwf_b(:,:) = pfwf(:,:) -
NEMO/trunk/src/OCE/ISF/isfutils.F90
r12271 r13286 12 12 !!---------------------------------------------------------------------- 13 13 14 USE iom , ONLY: iom_open, iom_get, iom_close, jpdom_ data! read input file14 USE iom , ONLY: iom_open, iom_get, iom_close, jpdom_global ! read input file 15 15 USE lib_fortran , ONLY: glob_sum, glob_min, glob_max ! compute global value 16 USE par_oce , ONLY: jpi,jpj,jpk, jpnij 17 USE dom_oce , ONLY: n ldi, nlei, nldj, nlej, narea, tmask_h, tmask_i! local domain16 USE par_oce , ONLY: jpi,jpj,jpk, jpnij, Nis0, Nie0, Njs0, Nje0 ! domain size 17 USE dom_oce , ONLY: narea, tmask_h, tmask_i ! local domain 18 18 USE in_out_manager, ONLY: i8, wp, lwp, numout ! miscelenious 19 19 USE lib_mpp … … 47 47 48 48 CALL iom_open( TRIM(cdfile), inum ) 49 CALL iom_get( inum, jpdom_ data, TRIM(cdvar), pvar)49 CALL iom_get( inum, jpdom_global, TRIM(cdvar), pvar) 50 50 CALL iom_close(inum) 51 51 … … 84 84 ! 85 85 ! local MOD sum 86 DO jj= nldj,nlej87 DO ji= nldi,nlei86 DO jj=Njs0,Nje0 87 DO ji=Nis0,Nie0 88 88 idums = ABS(MOD(TRANSFER(pvar(ji,jj), ip),imodd)) 89 89 itmps(narea) = MOD(itmps(narea) + idums, imods) … … 138 138 ! local MOD sum 139 139 DO jk=1,jpk 140 DO jj= nldj,nlej141 DO ji= nldi,nlei140 DO jj=Njs0,Nje0 141 DO ji=Nis0,Nie0 142 142 idums = ABS(MOD(TRANSFER(pvar(ji,jj,jk), ip),imodd)) 143 143 itmps(narea) = MOD(itmps(narea) + idums, imods)
Note: See TracChangeset
for help on using the changeset viewer.