- Timestamp:
- 2015-06-19T17:18:00+02:00 (9 years ago)
- Location:
- branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OFF_SRC
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OFF_SRC/domain.F90
r5234 r5443 116 116 USE ioipsl 117 117 INTEGER :: ios ! Local integer output status for namelist read 118 NAMELIST/namrun/ nn_no , cn_exp , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl, & 118 NAMELIST/namrun/ cn_ocerst_indir, cn_ocerst_outdir, nn_stocklist, ln_rst_list, & 119 & nn_no , cn_exp , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl, & 119 120 & nn_it000, nn_itend , nn_date0 , nn_leapy , nn_istate , nn_stock , & 120 & nn_write, ln_dimgnnn, ln_mskland , ln_c lobber, nn_chunksz, nn_euler121 & nn_write, ln_dimgnnn, ln_mskland , ln_cfmeta , ln_clobber, nn_chunksz, nn_euler 121 122 NAMELIST/namdom/ nn_bathy , rn_bathy, rn_e3zps_min, rn_e3zps_rat, nn_msh , rn_hmin, & 122 123 & nn_acc , rn_atfp , rn_rdt , rn_rdtmin , & … … 159 160 WRITE(numout,*) ' multi file dimgout ln_dimgnnn = ', ln_dimgnnn 160 161 WRITE(numout,*) ' mask land points ln_mskland = ', ln_mskland 162 WRITE(numout,*) ' additional CF standard metadata ln_cfmeta = ', ln_cfmeta 161 163 WRITE(numout,*) ' overwrite an existing file ln_clobber = ', ln_clobber 162 164 WRITE(numout,*) ' NetCDF chunksize (bytes) nn_chunksz = ', nn_chunksz … … 171 173 ninist = nn_istate 172 174 nstock = nn_stock 175 nstocklist = nn_stocklist 173 176 nwrite = nn_write 174 177 -
branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OFF_SRC/dommsk.F90
r5234 r5443 15 15 USE in_out_manager ! I/O manager 16 16 USE wrk_nemo 17 USE lbclnk 17 18 18 19 IMPLICIT NONE … … 48 49 !!---------------------------------------------------------------------- 49 50 ! 50 INTEGER :: ji, j k ! dummy loop indices51 INTEGER :: ji, jj, jk ! dummy loop indices 51 52 INTEGER :: iif, iil, ijf, ijl ! local integers 52 53 INTEGER, POINTER, DIMENSION(:,:) :: imsk … … 83 84 ENDIF 84 85 ! 86 ! (ISF) MIN(1,SUM(umask)) is here to check if you have effectively at 87 ! least 1 wet u point 88 DO jj = 1, jpjm1 89 DO ji = 1, fs_jpim1 ! vector loop 90 umask_i(ji,jj) = ssmask(ji,jj) * ssmask(ji+1,jj ) * MIN(1._wp,SUM(umask(ji,jj,:))) 91 vmask_i(ji,jj) = ssmask(ji,jj) * ssmask(ji ,jj+1) * MIN(1._wp,SUM(vmask(ji,jj,:))) 92 END DO 93 DO ji = 1, jpim1 ! NO vector opt. 94 fmask_i(ji,jj) = ssmask(ji,jj ) * ssmask(ji+1,jj ) & 95 & * ssmask(ji,jj+1) * ssmask(ji+1,jj+1) * MIN(1._wp,SUM(fmask(ji,jj,:))) 96 END DO 97 END DO 98 CALL lbc_lnk( umask_i, 'U', 1._wp ) ! Lateral boundary conditions 99 CALL lbc_lnk( vmask_i, 'V', 1._wp ) 100 CALL lbc_lnk( fmask_i, 'F', 1._wp ) 101 102 ! 3. Ocean/land mask at wu-, wv- and w points 103 !---------------------------------------------- 104 wmask (:,:,1) = tmask(:,:,1) ! ???????? 105 wumask(:,:,1) = umask(:,:,1) ! ???????? 106 wvmask(:,:,1) = vmask(:,:,1) ! ???????? 107 DO jk=2,jpk 108 wmask (:,:,jk)=tmask(:,:,jk) * tmask(:,:,jk-1) 109 wumask(:,:,jk)=umask(:,:,jk) * umask(:,:,jk-1) 110 wvmask(:,:,jk)=vmask(:,:,jk) * vmask(:,:,jk-1) 111 END DO 112 ! 85 113 IF( nprint == 1 .AND. lwp ) THEN ! Control print 86 114 imsk(:,:) = INT( tmask_i(:,:) ) -
branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90
r5234 r5443 245 245 tsn(:,:,:,jp_sal) = sf_dyn(jf_sal)%fnow(:,:,:) * tmask(:,:,:) ! salinity 246 246 ! 247 CALL eos ( tsn, rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop 247 ! 248 CALL eos ( tsn, rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop 249 CALL eos_rab( tsn, rab_n ) ! now local thermal/haline expension ratio at T-points 250 CALL bn2 ( tsn, rab_n, rn2 ) ! before Brunt-Vaisala frequency need for zdfmxl 251 252 rn2b(:,:,:) = rn2(:,:,:) ! need for zdfmxl 248 253 CALL zdf_mxl( kt ) ! In any case, we need mxl 249 254 ! … … 259 264 fr_i(:,:) = sf_dyn(jf_ice)%fnow(:,:,1) * tmask(:,:,1) ! Sea-ice fraction 260 265 qsr (:,:) = sf_dyn(jf_qsr)%fnow(:,:,1) * tmask(:,:,1) ! solar radiation 261 IF 266 IF( ln_dynrnf ) & 262 267 rnf (:,:) = sf_dyn(jf_rnf)%fnow(:,:,1) * tmask(:,:,1) ! river runoffs 263 268 … … 383 388 384 389 ! 385 IF 390 IF( ln_dynrnf ) THEN 386 391 jf_rnf = jfld + 1 ; jfld = jf_rnf 387 392 slf_d(jf_rnf) = sn_rnf … … 535 540 !!--------------------------------------------------------------------- 536 541 #if defined key_ldfslp && ! defined key_c1d 542 CALL eos ( pts, rhd, rhop, gdept_0(:,:,:) ) 537 543 CALL eos_rab( pts, rab_n ) ! now local thermal/haline expension ratio at T-points 538 544 CALL bn2 ( pts, rab_n, rn2 ) ! now Brunt-Vaisala 539 IF( ln_zps ) & ! Partial steps: before Horizontal DErivative 540 & CALL zps_hde( kt, jpts, pts, gtsu, gtsv, & ! Partial steps: before horizontal gradient 541 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & ! 542 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level 543 ! only gtsu, gtsv, rhd, gru , grv are used 544 545 546 ! ! of t, s, rd at the bottom ocean level 545 546 ! Partial steps: before Horizontal DErivative 547 IF( ln_zps .AND. .NOT. ln_isfcav) & 548 & CALL zps_hde ( kt, jpts, pts, gtsu, gtsv, & ! Partial steps: before horizontal gradient 549 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 550 IF( ln_zps .AND. ln_isfcav) & 551 & CALL zps_hde_isf( kt, jpts, pts, gtsu, gtsv, & ! Partial steps for top cell (ISF) 552 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 553 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the first ocean level 554 555 rn2b(:,:,:) = rn2(:,:,:) ! need for zdfmxl 547 556 CALL zdf_mxl( kt ) ! mixed layer depth 548 557 CALL ldf_slp( kt, rhd, rn2 ) ! slopes -
branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OFF_SRC/nemogcm.F90
- Property svn:keywords set to Id
r5234 r5443 149 149 & nn_bench, nn_timing 150 150 NAMELIST/namcfg/ cp_cfg, cp_cfz, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, & 151 & jpizoom, jpjzoom, jperio 151 & jpizoom, jpjzoom, jperio, ln_use_jattr 152 152 !!---------------------------------------------------------------------- 153 153 cltxt = '' 154 cxios_context = 'nemo' 154 155 ! 155 156 ! ! Open reference namelist and configuration namelist files … … 182 183 #if defined key_iomput 183 184 CALL xios_initialize( "nemo",return_comm=ilocal_comm ) 184 narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm ) ! Nodes selection185 narea = mynode( cltxt, 'output.namelist.dyn', numnam_ref, numnam_cfg, numond , nstop, ilocal_comm ) ! Nodes selection 185 186 #else 186 187 ilocal_comm = 0 187 narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop ) ! Nodes selection (control print return in cltxt)188 narea = mynode( cltxt, 'output.namelist.dyn', numnam_ref, numnam_cfg, numond , nstop ) ! Nodes selection (control print return in cltxt) 188 189 #endif 189 190 … … 233 234 WRITE(numout,*) ' NEMO team' 234 235 WRITE(numout,*) ' Ocean General Circulation Model' 235 WRITE(numout,*) ' version 3. 5 (2012) '236 WRITE(numout,*) ' version 3.6 (2015) ' 236 237 WRITE(numout,*) 237 238 WRITE(numout,*) … … 359 360 WRITE(numout,*) ' left bottom j index of the zoom (in data domain) jpizoom = ', jpjzoom 360 361 WRITE(numout,*) ' lateral cond. type (between 0 and 6) jperio = ', jperio 362 WRITE(numout,*) ' use file attribute if exists as i/p j-start ln_use_jattr = ', ln_use_jattr 361 363 ENDIF 362 364 ! ! Parameter control
Note: See TracChangeset
for help on using the changeset viewer.