Changeset 11027 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/SAS
- Timestamp:
- 2019-05-21T17:33:54+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/SAS
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/SAS/diawri.F90
r10425 r11027 78 78 79 79 80 SUBROUTINE dia_wri( kt )80 SUBROUTINE dia_wri( kt, Kmm ) 81 81 !!--------------------------------------------------------------------- 82 82 !! *** ROUTINE dia_wri *** … … 90 90 !! 91 91 INTEGER, INTENT( in ) :: kt ! ocean time-step index 92 INTEGER, INTENT( in ) :: Kmm ! ocean time levelindex 92 93 !!---------------------------------------------------------------------- 93 94 ! 94 95 ! Output the initial state and forcings 95 96 IF( ninist == 1 ) THEN 96 CALL dia_wri_state( 'output.init' )97 CALL dia_wri_state( 'output.init', Kmm ) 97 98 ninist = 0 98 99 ENDIF … … 330 331 #endif 331 332 332 SUBROUTINE dia_wri_state( cdfile_name )333 SUBROUTINE dia_wri_state( cdfile_name, Kmm ) 333 334 !!--------------------------------------------------------------------- 334 335 !! *** ROUTINE dia_wri_state *** … … 344 345 !!---------------------------------------------------------------------- 345 346 CHARACTER (len=* ), INTENT( in ) :: cdfile_name ! name of the file created 347 INTEGER , INTENT( in ) :: Kmm ! ocean time levelindex 346 348 !! 347 349 INTEGER :: inum … … 359 361 #endif 360 362 361 CALL iom_rstput( 0, 0, inum, 'votemper', ts n(:,:,:,jp_tem) ) ! now temperature362 CALL iom_rstput( 0, 0, inum, 'vosaline', ts n(:,:,:,jp_sal) ) ! now salinity363 CALL iom_rstput( 0, 0, inum, 'sossheig', ssh n) ! sea surface height364 CALL iom_rstput( 0, 0, inum, 'vozocrtx', u n) ! now i-velocity365 CALL iom_rstput( 0, 0, inum, 'vomecrty', v n) ! now j-velocity366 CALL iom_rstput( 0, 0, inum, 'vovecrtz', w n) ! now k-velocity367 CALL iom_rstput( 0, 0, inum, 'sowaflup', emp - rnf ) ! freshwater budget368 CALL iom_rstput( 0, 0, inum, 'sohefldo', qsr + qns ) ! total heat flux369 CALL iom_rstput( 0, 0, inum, 'soshfldo', qsr ) ! solar heat flux370 CALL iom_rstput( 0, 0, inum, 'soicecov', fr_i ) ! ice fraction371 CALL iom_rstput( 0, 0, inum, 'sozotaux', utau ) ! i-wind stress372 CALL iom_rstput( 0, 0, inum, 'sometauy', vtau ) ! j-wind stress363 CALL iom_rstput( 0, 0, inum, 'votemper', ts (:,:,:,jp_tem,Kmm) ) ! now temperature 364 CALL iom_rstput( 0, 0, inum, 'vosaline', ts (:,:,:,jp_sal,Kmm) ) ! now salinity 365 CALL iom_rstput( 0, 0, inum, 'sossheig', ssh(:,:, Kmm) ) ! sea surface height 366 CALL iom_rstput( 0, 0, inum, 'vozocrtx', uu (:,:,:, Kmm) ) ! now i-velocity 367 CALL iom_rstput( 0, 0, inum, 'vomecrty', vv (:,:,:, Kmm) ) ! now j-velocity 368 CALL iom_rstput( 0, 0, inum, 'vovecrtz', ww ) ! now k-velocity 369 CALL iom_rstput( 0, 0, inum, 'sowaflup', emp - rnf ) ! freshwater budget 370 CALL iom_rstput( 0, 0, inum, 'sohefldo', qsr + qns ) ! total heat flux 371 CALL iom_rstput( 0, 0, inum, 'soshfldo', qsr ) ! solar heat flux 372 CALL iom_rstput( 0, 0, inum, 'soicecov', fr_i ) ! ice fraction 373 CALL iom_rstput( 0, 0, inum, 'sozotaux', utau ) ! i-wind stress 374 CALL iom_rstput( 0, 0, inum, 'sometauy', vtau ) ! j-wind stress 373 375 374 376 #if defined key_si3 -
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/SAS/nemogcm.F90
r10998 r11027 355 355 Nbb = 1; Nnn = 2; Naa = 3; Nrhs = Naa 356 356 357 ! Initialisation of temporary pointers (to be deleted after development finished)358 CALL update_pointers( Nbb, Nnn, Naa )359 357 ! !-------------------------------! 360 358 ! ! NEMO general initialization ! -
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/SAS/sbcssm.F90
r10922 r11027 121 121 IF( .NOT. ln_linssh ) e3t_m(:,:) = e3t_0(:,:,1) !clem: necessary at least for sas2D 122 122 frq_m(:,:) = 1._wp ! - - 123 ssh n (:,:) = 0._wp ! - -123 ssh (:,:,Kmm) = 0._wp ! - - 124 124 ENDIF 125 125 126 126 IF ( nn_ice == 1 ) THEN 127 ts n(:,:,1,jp_tem) = sst_m(:,:)128 ts n(:,:,1,jp_sal) = sss_m(:,:)129 ts b(:,:,1,jp_tem) = sst_m(:,:)130 ts b(:,:,1,jp_sal) = sss_m(:,:)131 ENDIF 132 u b (:,:,1) = ssu_m(:,:)133 v b (:,:,1) = ssv_m(:,:)127 ts(:,:,1,jp_tem,Kmm) = sst_m(:,:) 128 ts(:,:,1,jp_sal,Kmm) = sss_m(:,:) 129 ts(:,:,1,jp_tem,Kbb) = sst_m(:,:) 130 ts(:,:,1,jp_sal,Kbb) = sss_m(:,:) 131 ENDIF 132 uu (:,:,1,Kbb) = ssu_m(:,:) 133 vv (:,:,1,Kbb) = ssv_m(:,:) 134 134 135 135 IF(ln_ctl) THEN ! print control -
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/SAS/step.F90
r10975 r11027 47 47 48 48 PUBLIC stp ! called by nemogcm.F90 49 PUBLIC update_pointers ! called by nemo_init50 49 51 50 !!---------------------------------------------------------------------- … … 105 104 CALL sbc ( kstp, Nbb, Nnn ) ! Sea Boundary Condition (including sea-ice) 106 105 107 CALL dia_wri( kstp )! ocean model: outputs106 CALL dia_wri( kstp, Nnn ) ! ocean model: outputs 108 107 109 108 #if defined key_agrif … … 126 125 IF( indic < 0 ) THEN 127 126 CALL ctl_stop( 'step: indic < 0' ) 128 CALL dia_wri_state( 'output.abort' )127 CALL dia_wri_state( 'output.abort', Nnn ) 129 128 ENDIF 130 IF( kstp == nit000 ) CALL iom_close( numror ) ! close input ocean restart file129 IF( kstp == nit000 ) CALL iom_close( numror ) ! close input ocean restart file 131 130 132 131 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 133 132 ! Coupled mode 134 133 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 135 IF( lk_oasis ) CALL sbc_cpl_snd( kstp, N nn ) ! coupled mode : field exchanges if OASIS-coupled ice134 IF( lk_oasis ) CALL sbc_cpl_snd( kstp, Nbb, Nnn ) ! coupled mode : field exchanges if OASIS-coupled ice 136 135 137 136 #if defined key_iomput … … 153 152 END SUBROUTINE stp 154 153 155 SUBROUTINE update_pointers( Kbb, Kmm, Kaa )156 !!----------------------------------------------------------------------157 !! *** ROUTINE update_pointers ***158 !!159 !! ** Purpose : Associate temporary pointer arrays.160 !! For IMMERSE development phase only - to be deleted161 !!162 !! ** Method :163 !!----------------------------------------------------------------------164 INTEGER, INTENT( in ) :: Kbb, Kmm, Kaa ! time level indices165 166 ub => uu(:,:,:,Kbb); un => uu(:,:,:,Kmm); ua => uu(:,:,:,Kaa)167 vb => vv(:,:,:,Kbb); vn => vv(:,:,:,Kmm); va => vv(:,:,:,Kaa)168 wn => ww(:,:,:)169 hdivn => hdiv(:,:,:)170 171 sshb => ssh(:,:,Kbb); sshn => ssh(:,:,Kmm); ssha => ssh(:,:,Kaa)172 ub_b => uu_b(:,:,Kbb); un_b => uu_b(:,:,Kmm); ua_b => uu_b(:,:,Kaa)173 vb_b => vv_b(:,:,Kbb); vn_b => vv_b(:,:,Kmm); va_b => vv_b(:,:,Kaa)174 175 tsb => ts(:,:,:,:,Kbb); tsn => ts(:,:,:,:,Kmm); tsa => ts(:,:,:,:,Kaa)176 177 e3t_b => e3t(:,:,:,Kbb); e3t_n => e3t(:,:,:,Kmm); e3t_a => e3t(:,:,:,Kaa)178 e3u_b => e3u(:,:,:,Kbb); e3u_n => e3u(:,:,:,Kmm); e3u_a => e3u(:,:,:,Kaa)179 e3v_b => e3v(:,:,:,Kbb); e3v_n => e3v(:,:,:,Kmm); e3v_a => e3v(:,:,:,Kaa)180 181 e3f_n => e3f(:,:,:)182 183 e3w_b => e3w (:,:,:,Kbb); e3w_n => e3w (:,:,:,Kmm)184 e3uw_b => e3uw(:,:,:,Kbb); e3uw_n => e3uw(:,:,:,Kmm)185 e3vw_b => e3vw(:,:,:,Kbb); e3vw_n => e3vw(:,:,:,Kmm)186 187 gdept_b => gdept(:,:,:,Kbb); gdept_n => gdept(:,:,:,Kmm)188 gdepw_b => gdepw(:,:,:,Kbb); gdepw_n => gdepw(:,:,:,Kmm)189 gde3w_n => gde3w(:,:,:)190 191 END SUBROUTINE update_pointers192 193 154 !!====================================================================== 194 155 END MODULE step
Note: See TracChangeset
for help on using the changeset viewer.