Changeset 899 for codes/icosagcm/trunk/src/time
- Timestamp:
- 06/13/19 16:45:41 (5 years ago)
- Location:
- codes/icosagcm/trunk/src/time
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/time/euler_scheme.f90
r548 r899 36 36 REAL(rstd),POINTER :: hflux(:,:),wflux(:,:),hfluxt(:,:),wfluxt(:,:) 37 37 INTEGER :: ind 38 INTEGER :: i ,j,ij,l38 INTEGER :: ij,l 39 39 CALL trace_start("Euler_scheme") 40 40 … … 95 95 REAL(rstd), INTENT(IN) :: tau 96 96 LOGICAL, INTENT(INOUT) :: fluxt_zero 97 INTEGER :: l,i ,j,ij97 INTEGER :: l,ij 98 98 99 99 IF(fluxt_zero) THEN -
codes/icosagcm/trunk/src/time/explicit_scheme.f90
r548 r899 22 22 ! USE caldyn_gcm_mod, ONLY : req_ps, req_mass 23 23 24 REAL(rstd),POINTER :: q(:,:,:) 25 REAL(rstd),POINTER :: phis(:), ps(:) ,psm1(:), psm2(:), dps(:) 24 REAL(rstd),POINTER :: ps(:) ,psm1(:), psm2(:), dps(:) 26 25 REAL(rstd),POINTER :: u(:,:) , um1(:,:), um2(:,:), du(:,:) 27 REAL(rstd),POINTER :: rhodz(:,:), mass(:,:), massm1(:,:), massm2(:,:), dmass(:,:)26 REAL(rstd),POINTER :: mass(:,:), massm1(:,:), dmass(:,:) 28 27 REAL(rstd),POINTER :: theta_rhodz(:,:,:) , theta_rhodzm1(:,:,:), theta_rhodzm2(:,:,:), dtheta_rhodz(:,:,:) 29 28 REAL(rstd),POINTER :: hflux(:,:),wflux(:,:),hfluxt(:,:),wfluxt(:,:) … … 64 63 REAL(rstd), INTENT(IN) :: coef(:) 65 64 REAL(rstd) :: tau 66 INTEGER :: i ,j,ij,l65 INTEGER :: ij,l 67 66 68 67 CALL trace_start("RK_scheme") -
codes/icosagcm/trunk/src/time/time.f90
r706 r899 54 54 run_length=dt*itaumax 55 55 CALL getin('run_length',run_length) 56 itaumax= run_length/dt56 itaumax=INT(run_length/dt) 57 57 58 58 time_style='dcmip' -
codes/icosagcm/trunk/src/time/timeloop_gcm.f90
r895 r899 194 194 195 195 REAL(rstd) :: adv_over_out ! ratio itau_adv/itau_out 196 INTEGER :: ind, it, i,j,l,n, stage196 INTEGER :: ind, it,l 197 197 LOGICAL :: fluxt_zero(ndomain) ! set to .TRUE. to start accumulating mass fluxes in time 198 198 LOGICAL, PARAMETER :: check_rhodz=.FALSE. … … 401 401 elapsed = (stop_clock-start_clock)*1./rate_clock 402 402 per_step = elapsed/(it-itau0) 403 throughput = dt/per_step403 throughput = INT(dt/per_step) 404 404 total = per_step*itaumax 405 405 WRITE(*,'(A,I5,A,F6.2,A,I6)') 'Time spent (s):',INT(elapsed), &
Note: See TracChangeset
for help on using the changeset viewer.