- Timestamp:
- 2022-02-04T17:43:35+01:00 (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.3_GO8_package/src/OCE/DIA/diaprod.F90
r13590 r15695 65 65 ! 66 66 67 IF( iom_use("urhop") .OR. iom_use("vrhop") .OR. iom_use("wrhop") & 68 #if ! defined key_diaar5 69 & .OR. iom_use("rhop") & 70 #endif 71 & ) THEN 67 IF( iom_use("urhop") .OR. iom_use("vrhop") .OR. iom_use("wrhop") ) THEN 72 68 CALL eos( tsn, z3d, zrhop ) ! now in situ and potential density 73 69 zrhop(:,:,:) = zrhop(:,:,:)-1000.e0 ! reference potential density to 1000 to avoid precision issues in rhop2 calculation 74 70 zrhop(:,:,jpk) = 0._wp 75 #if ! defined key_diaar576 CALL iom_put( 'rhop', zrhop )77 #else78 ! If key_diaar5 set then there is already an iom_put call to output rhop.79 ! Really should be a standard diagnostics option?80 #endif81 71 ENDIF 82 72
Note: See TracChangeset
for help on using the changeset viewer.