- Timestamp:
- 2011-03-15T16:27:46+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r2636 r2690 31 31 USE dom_oce ! ocean space and time domain 32 32 USE phycst ! physical constants 33 USE in_out_manager ! I/O manager34 33 USE trdmod ! ocean dynamics trends 35 34 USE trdmod_oce ! ocean variables trends 35 USE in_out_manager ! I/O manager 36 36 USE prtctl ! Print control 37 37 USE lbclnk ! lateral boundary condition … … 77 77 !! - Save the trend (l_trddyn=T) 78 78 !!---------------------------------------------------------------------- 79 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released80 USE wrk_nemo, ONLY: ztrdu => wrk_3d_1, ztrdv => wrk_3d_2 ! 3D workspace79 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 80 USE wrk_nemo, ONLY: ztrdu => wrk_3d_1 , ztrdv => wrk_3d_2 ! 3D workspace 81 81 !! 82 82 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 193 193 !! ** Action : - Update (ua,va) with the now hydrastatic pressure trend 194 194 !!---------------------------------------------------------------------- 195 USE oce, ONLY : zhpi => ta ! use ta as 3D workspace 196 USE oce, ONLY : zhpj => sa ! use sa as 3D workspace 195 USE oce, ONLY: zhpi => ta , zhpj => sa ! (ta,sa) used as 3D workspace 197 196 !! 198 197 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 208 207 ENDIF 209 208 210 ! Local constant initialization 211 zcoef0 = - grav * 0.5_wp 209 zcoef0 = - grav * 0.5_wp ! Local constant initialization 212 210 213 211 ! Surface value … … 255 253 !! ** Action : - Update (ua,va) with the now hydrastatic pressure trend 256 254 !!---------------------------------------------------------------------- 257 USE oce, ONLY : zhpi => ta ! use ta as 3D workspace 258 USE oce, ONLY : zhpj => sa ! use sa as 3D workspace 255 USE oce, ONLY: zhpi => ta , zhpj => sa ! (ta,sa) used as 3D workspace 259 256 !! 260 257 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 357 354 !! ** Action : - Update (ua,va) with the now hydrastatic pressure trend 358 355 !!---------------------------------------------------------------------- 359 USE oce, ONLY : zhpi => ta ! use ta as 3D workspace 360 USE oce, ONLY : zhpj => sa ! use sa as 3D workspace 356 USE oce, ONLY: zhpi => ta , zhpj => sa ! (ta,sa) used as 3D workspace 361 357 !! 362 358 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 443 439 !! - Save the trend (l_trddyn=T) 444 440 !!---------------------------------------------------------------------- 445 USE oce, ONLY : zhpi => ta ! use ta as 3D workspace 446 USE oce, ONLY : zhpj => sa ! use sa as 3D workspace 441 USE oce, ONLY: zhpi => ta , zhpj => sa ! (ta,sa) used as 3D workspace 447 442 !! 448 443 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 520 515 !! Reference : Song, Mon. Wea. Rev., 126, 3213-3230, 1998. 521 516 !!---------------------------------------------------------------------- 522 USE oce, ONLY : zhpi => ta ! use ta as 3D workspace 523 USE oce, ONLY : zhpj => sa ! use sa as 3D workspace 517 USE oce, ONLY: zhpi => ta , zhpj => sa ! (ta,sa) used as 3D workspace 524 518 !! 525 519 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 600 594 !! Reference: Shchepetkin and McWilliams, J. Geophys. Res., 108(C3), 3090, 2003 601 595 !!---------------------------------------------------------------------- 602 USE oce, ONLY : zhpi => ta ! use ta as 3D workspace 603 USE oce, ONLY : zhpj => sa ! use sa as 3D workspace 604 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 605 USE wrk_nemo, ONLY: drhox => wrk_3d_1 , dzx => wrk_3d_2 606 USE wrk_nemo, ONLY: drhou => wrk_3d_3 , dzu => wrk_3d_4 , rho_i => wrk_3d_5 607 USE wrk_nemo, ONLY: drhoy => wrk_3d_6 , dzy => wrk_3d_7 608 USE wrk_nemo, ONLY: drhov => wrk_3d_8 , dzv => wrk_3d_9 , rho_j => wrk_3d_10 609 USE wrk_nemo, ONLY: drhoz => wrk_3d_11 , dzz => wrk_3d_12 610 USE wrk_nemo, ONLY: drhow => wrk_3d_13 , dzw => wrk_3d_14 611 USE wrk_nemo, ONLY: rho_k => wrk_3d_15 596 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 597 USE oce , ONLY: zhpi => ta , zhpj => sa ! (ta,sa) used as 3D workspace 598 USE wrk_nemo, ONLY: drhox => wrk_3d_1 , dzx => wrk_3d_2 599 USE wrk_nemo, ONLY: drhou => wrk_3d_3 , dzu => wrk_3d_4 , rho_i => wrk_3d_5 600 USE wrk_nemo, ONLY: drhoy => wrk_3d_6 , dzy => wrk_3d_7 601 USE wrk_nemo, ONLY: drhov => wrk_3d_8 , dzv => wrk_3d_9 , rho_j => wrk_3d_10 602 USE wrk_nemo, ONLY: drhoz => wrk_3d_11 , dzz => wrk_3d_12 603 USE wrk_nemo, ONLY: drhow => wrk_3d_13 , dzw => wrk_3d_14 604 USE wrk_nemo, ONLY: rho_k => wrk_3d_15 612 605 !! 613 606 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 620 613 621 614 IF( wrk_in_use(3, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15) ) THEN 622 CALL ctl_stop('dyn:hpg_djc 615 CALL ctl_stop('dyn:hpg_djc: requested workspace arrays unavailable') ; RETURN 623 616 ENDIF 624 617 … … 628 621 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, density Jacobian with cubic polynomial scheme' 629 622 ENDIF 630 631 623 632 624 ! Local constant initialization … … 820 812 ! 821 813 IF( wrk_not_released(3, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15) ) & 822 CALL ctl_stop('dyn:hpg_djc 814 CALL ctl_stop('dyn:hpg_djc: failed to release workspace arrays') 823 815 ! 824 816 END SUBROUTINE hpg_djc … … 833 825 !! Reference: Thiem & Berntsen, Ocean Modelling, In press, 2005. 834 826 !!---------------------------------------------------------------------- 835 USE oce, ONLY : zhpi => ta ! use ta as 3D workspace 836 USE oce, ONLY : zhpj => sa ! use sa as 3D workspace 837 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 838 USE wrk_nemo, ONLY: zdistr => wrk_2d_1 , zsina => wrk_2d_2 , zcosa => wrk_2d_3 839 USE wrk_nemo, ONLY: zhpiorg => wrk_3d_1 , zhpirot => wrk_3d_2 840 USE wrk_nemo, ONLY: zhpitra => wrk_3d_3 , zhpine => wrk_3d_4 841 USE wrk_nemo, ONLY: zhpjorg => wrk_3d_5 , zhpjrot => wrk_3d_6 842 USE wrk_nemo, ONLY: zhpjtra => wrk_3d_7 , zhpjne => wrk_3d_8 827 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 828 USE oce , ONLY: zhpi => ta , zhpj => sa ! (ta,sa) used as 3D workspace 829 USE wrk_nemo, ONLY: zdistr => wrk_2d_1 , zsina => wrk_2d_2 , zcosa => wrk_2d_3 830 USE wrk_nemo, ONLY: zhpiorg => wrk_3d_1 , zhpirot => wrk_3d_2 831 USE wrk_nemo, ONLY: zhpitra => wrk_3d_3 , zhpine => wrk_3d_4 832 USE wrk_nemo, ONLY: zhpjorg => wrk_3d_5 , zhpjrot => wrk_3d_6 833 USE wrk_nemo, ONLY: zhpjtra => wrk_3d_7 , zhpjne => wrk_3d_8 843 834 !! 844 835 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 849 840 !!---------------------------------------------------------------------- 850 841 851 IF( wrk_in_use(2, 1,2,3) .OR.&842 IF( wrk_in_use(2, 1,2,3) .OR. & 852 843 wrk_in_use(3, 1,2,3,4,5,6,7,8) ) THEN 853 CALL ctl_stop('dyn:hpg_rot 854 END 844 CALL ctl_stop('dyn:hpg_rot: requested workspace arrays unavailable') ; RETURN 845 ENDIF 855 846 856 847 IF( kt == nit000 ) THEN … … 1009 1000 END DO 1010 1001 ! 1011 IF( wrk_not_released(2, 1,2,3) .OR.&1012 wrk_not_released(3, 1,2,3,4,5,6,7,8) ) CALL ctl_stop('dyn:hpg_rot 1002 IF( wrk_not_released(2, 1,2,3) .OR. & 1003 wrk_not_released(3, 1,2,3,4,5,6,7,8) ) CALL ctl_stop('dyn:hpg_rot: failed to release workspace arrays') 1013 1004 ! 1014 1005 END SUBROUTINE hpg_rot
Note: See TracChangeset
for help on using the changeset viewer.