diff --git a/src/arome/aux/mode_budget.F90 b/src/arome/aux/mode_budget_phy.F90
similarity index 99%
rename from src/arome/aux/mode_budget.F90
rename to src/arome/aux/mode_budget_phy.F90
index 3a258b2443dee2db8d5da4990db60e41534e6ed8..d2881a9679bcbf83eaa50a3fc3503bf7f5cfaf50 100644
--- a/src/arome/aux/mode_budget.F90
+++ b/src/arome/aux/mode_budget_phy.F90
@@ -1,4 +1,4 @@
-MODULE MODE_BUDGET
+MODULE MODE_BUDGET_PHY
USE MODD_BUDGET, ONLY : TBUDGETDATA
IMPLICIT NONE
CONTAINS
@@ -291,5 +291,4 @@ ENDIF
IF (LHOOK) CALL DR_HOOK('BUDGET_DDH',1,ZHOOK_HANDLE)
END SUBROUTINE BUDGET_DDH
-END MODULE MODE_BUDGET
-
+END MODULE MODE_BUDGET_PHY
diff --git a/src/arome/ext/aro_adjust_lima.F90 b/src/arome/ext/aro_adjust_lima.F90
index b7854d832bc629ca86ec37bb54e7266f26dafa66..018e5aa73e1d5aa935663ea77478e97974acf1a3 100644
--- a/src/arome/ext/aro_adjust_lima.F90
+++ b/src/arome/ext/aro_adjust_lima.F90
@@ -91,7 +91,7 @@ USE MODI_LIMA_ADJUST_SPLIT
USE MODE_SET_CONC_LIMA
USE MODE_SET_CONC_LIMA_LBC
USE MODE_FILL_DIMPHYEX, ONLY: FILL_DIMPHYEX
-!USE MODE_BUDGET, ONLY: BUDGET_DDH
+!USE MODE_BUDGET_PHY, ONLY: BUDGET_DDH
!
USE DDH_MIX, ONLY : TYP_DDH
USE YOMLDDH, ONLY : TLDDH
diff --git a/src/arome/ext/aro_convbu.F90 b/src/arome/ext/aro_convbu.F90
index 6cb03579814bc14aebf704684fc84775b686ace3..06d195d23d961020f056e99f94ab0be57c7b91de 100644
--- a/src/arome/ext/aro_convbu.F90
+++ b/src/arome/ext/aro_convbu.F90
@@ -3,7 +3,7 @@ USE PARKIND1, ONLY : JPRB
USE YOMHOOK , ONLY : LHOOK, DR_HOOK
USE MODD_BUDGET,ONLY: LBU_ENABLE, LBUDGET_TH, LBUDGET_RV, LBUDGET_RC, LBUDGET_RI
-USE MODE_BUDGET, ONLY: BUDGET_DDH
+USE MODE_BUDGET_PHY, ONLY: BUDGET_DDH
USE DDH_MIX, ONLY : TYP_DDH
USE YOMLDDH, ONLY : TLDDH
USE YOMMDDH, ONLY : TMDDH
diff --git a/src/arome/ext/aro_lima.F90 b/src/arome/ext/aro_lima.F90
index edef5e82779c4dfc8784d9ed5ca29077fd226722..6c8b4bb4089757a029990dc9aa02547e76811a54 100644
--- a/src/arome/ext/aro_lima.F90
+++ b/src/arome/ext/aro_lima.F90
@@ -51,7 +51,7 @@ USE MODD_PARAM_LIMA
USE MODD_NSV
!
USE MODD_BUDGET
-USE MODE_BUDGET, ONLY: BUDGET_DDH
+USE MODE_BUDGET_PHY, ONLY: BUDGET_DDH
USE MODE_FILL_DIMPHYEX, ONLY: FILL_DIMPHYEX
!
USE MODI_LIMA
diff --git a/src/arome/ext/aro_rain_ice.F90 b/src/arome/ext/aro_rain_ice.F90
index 9fbb932871d0eb76d26a2654e050177db459a235..8fd95fc4d9f092f209e61990836d98d4c7ad97c2 100644
--- a/src/arome/ext/aro_rain_ice.F90
+++ b/src/arome/ext/aro_rain_ice.F90
@@ -94,7 +94,7 @@ USE MODD_PARAM_ICE, ONLY: PARAM_ICE
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
USE MODD_BUDGET, ONLY: TBUDGETDATA, NBUDGET_RH, TBUCONF
-USE MODE_BUDGET, ONLY: BUDGET_DDH
+USE MODE_BUDGET_PHY, ONLY: BUDGET_DDH
USE MODE_FILL_DIMPHYEX, ONLY: FILL_DIMPHYEX
!
USE MODI_RAIN_ICE
diff --git a/src/arome/ext/aro_startbu.F90 b/src/arome/ext/aro_startbu.F90
index 8eebeb0fca08d209f47737ac8847c262c14e73c0..c9b10fcdeb1cf54c83f74f8466f8e16a4efdf7d4 100644
--- a/src/arome/ext/aro_startbu.F90
+++ b/src/arome/ext/aro_startbu.F90
@@ -44,7 +44,7 @@ USE YOMHOOK , ONLY : LHOOK, DR_HOOK
! ------------
!
USE MODD_BUDGET
-USE MODE_BUDGET, ONLY: BUDGET_DDH
+USE MODE_BUDGET_PHY, ONLY: BUDGET_DDH
USE DDH_MIX , ONLY : TYP_DDH
USE YOMLDDH , ONLY : TLDDH
USE YOMMDDH , ONLY : TMDDH
diff --git a/src/arome/micro/rain_ice_old.F90 b/src/arome/micro/rain_ice_old.F90
index e7b5994e23cda700e473725afa0ac7ea66e9fbcc..e291da334a50bfc82c41c8b99d4d1f35f0d5a6d5 100644
--- a/src/arome/micro/rain_ice_old.F90
+++ b/src/arome/micro/rain_ice_old.F90
@@ -173,7 +173,7 @@ USE MODD_RAIN_ICE_PARAM
USE MODD_PARAM_ICE
USE MODD_BUDGET
USE MODD_LES
-USE MODE_BUDGET, ONLY: BUDGET_DDH
+USE MODE_BUDGET_PHY, ONLY: BUDGET_DDH
USE MODI_GAMMA
USE MODE_TIWMX
USE MODE_ICECLOUD, ONLY : ICECLOUD
diff --git a/src/common/aux/mode_budget.F90 b/src/common/aux/mode_budget_phy.F90
similarity index 98%
rename from src/common/aux/mode_budget.F90
rename to src/common/aux/mode_budget_phy.F90
index e50f2185b75c75ceb45d12d9de4de0c746b02456..30db2f33b4fc10006c42562750ef1a5c7ddb131b 100644
--- a/src/common/aux/mode_budget.F90
+++ b/src/common/aux/mode_budget_phy.F90
@@ -1,4 +1,4 @@
-MODULE MODE_BUDGET
+MODULE MODE_BUDGET_PHY
USE MODD_BUDGET, ONLY : TBUDGETDATA
IMPLICIT NONE
CONTAINS
@@ -48,4 +48,4 @@ SUBROUTINE BUDGET_STORE_ADD(TPBUDGET, HSOURCE, PVARS)
REAL, DIMENSION(:,:,:), INTENT(IN) :: PVARS ! Current value to be stored
REAL, DIMENSION(SIZE(PVARS, 1), SIZE(PVARS, 2), SIZE(PVARS, 3)) :: ZVARS
END SUBROUTINE BUDGET_STORE_ADD
-END MODULE MODE_BUDGET
+END MODULE MODE_BUDGET_PHY
diff --git a/src/common/aux/gradient_m_phy.F90 b/src/common/aux/mode_gradient_m_phy.F90
similarity index 99%
rename from src/common/aux/gradient_m_phy.F90
rename to src/common/aux/mode_gradient_m_phy.F90
index 348dea22e2e0f2301f5b7e15a5ad8dd6bb750b48..6f3c34f01c872849d90154f8b73a90b649854bd3 100644
--- a/src/common/aux/gradient_m_phy.F90
+++ b/src/common/aux/mode_gradient_m_phy.F90
@@ -166,7 +166,7 @@ SUBROUTINE GX_M_M_PHY(D,OFLAT,PA,PDXX,PDZZ,PDZX,PGX_M_M)
!
!
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
-USE SHUMAN_PHY, ONLY: DXF_PHY, MZF_PHY, DZM_PHY, MXF_PHY, MXM_PHY
+USE MODE_SHUMAN_PHY, ONLY: DXF_PHY, MZF_PHY, DZM_PHY, MXF_PHY, MXM_PHY
!
IMPLICIT NONE
!
@@ -290,7 +290,7 @@ END SUBROUTINE GX_M_M_PHY
!
!
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
-USE SHUMAN_PHY, ONLY: DYF_PHY, MZF_PHY, DZM_PHY, MYF_PHY, MYM_PHY
+USE MODE_SHUMAN_PHY, ONLY: DYF_PHY, MZF_PHY, DZM_PHY, MYF_PHY, MYM_PHY
!
IMPLICIT NONE
!
diff --git a/src/common/aux/gradient_u_phy.F90 b/src/common/aux/mode_gradient_u_phy.F90
similarity index 98%
rename from src/common/aux/gradient_u_phy.F90
rename to src/common/aux/mode_gradient_u_phy.F90
index ff685a0c1319532998fddd9a0051fd17dd1fc4eb..f66dfeff79fac30aeabb45fe0d91866866caa170 100644
--- a/src/common/aux/gradient_u_phy.F90
+++ b/src/common/aux/mode_gradient_u_phy.F90
@@ -53,7 +53,7 @@ CONTAINS
!* 0. DECLARATIONS
!
!
-USE SHUMAN_PHY, ONLY : DZM_PHY, MXM_PHY
+USE MODE_SHUMAN_PHY, ONLY : DZM_PHY, MXM_PHY
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
@@ -158,7 +158,7 @@ END SUBROUTINE GZ_U_UW_PHY
!* 0. DECLARATIONS
!
!
-USE SHUMAN_PHY, ONLY : DZM_PHY, DXF_PHY, MXF_PHY, MZF_PHY
+USE MODE_SHUMAN_PHY, ONLY : DZM_PHY, DXF_PHY, MXF_PHY, MZF_PHY
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
diff --git a/src/common/aux/gradient_v_phy.F90 b/src/common/aux/mode_gradient_v_phy.F90
similarity index 98%
rename from src/common/aux/gradient_v_phy.F90
rename to src/common/aux/mode_gradient_v_phy.F90
index 66ec0b4ca7708bd3dffd385752e336a49910f9d6..0e29c0064d67e01b07396f7c7e544ae9968592b6 100644
--- a/src/common/aux/gradient_v_phy.F90
+++ b/src/common/aux/mode_gradient_v_phy.F90
@@ -54,7 +54,7 @@ CONTAINS
!* 0. DECLARATIONS
!
!
-USE SHUMAN_PHY, ONLY : DZM_PHY, MYM_PHY
+USE MODE_SHUMAN_PHY, ONLY : DZM_PHY, MYM_PHY
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
@@ -154,7 +154,7 @@ END SUBROUTINE GZ_V_VW_PHY
!* 0. DECLARATIONS
!
!
-USE SHUMAN_PHY, ONLY : DZM_PHY, DYF_PHY, MYF_PHY, MZF_PHY
+USE MODE_SHUMAN_PHY, ONLY : DZM_PHY, DYF_PHY, MYF_PHY, MZF_PHY
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
diff --git a/src/common/aux/gradient_w_phy.F90 b/src/common/aux/mode_gradient_w_phy.F90
similarity index 98%
rename from src/common/aux/gradient_w_phy.F90
rename to src/common/aux/mode_gradient_w_phy.F90
index 0377f0991ceffd5f90533f5af1091b9bc653dc36..0f8db721bfbc5be4d986dd23373fb00eb0f36b2f 100644
--- a/src/common/aux/gradient_w_phy.F90
+++ b/src/common/aux/mode_gradient_w_phy.F90
@@ -50,7 +50,7 @@ CONTAINS
!* 0. DECLARATIONS
!
!
-USE SHUMAN_PHY, ONLY: MZF_PHY, DZF_PHY, MXM_PHY, DXM_PHY, MZM_PHY, DZM_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZF_PHY, DZF_PHY, MXM_PHY, DXM_PHY, MZM_PHY, DZM_PHY
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
@@ -162,7 +162,7 @@ END SUBROUTINE GX_W_UW_PHY
!* 0. DECLARATIONS
!
!
-USE SHUMAN_PHY, ONLY: MZF_PHY, DZF_PHY, MYM_PHY, DYM_PHY, MZM_PHY, DZM_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZF_PHY, DZF_PHY, MYM_PHY, DYM_PHY, MZM_PHY, DZM_PHY
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
@@ -280,7 +280,7 @@ END SUBROUTINE GY_W_VW_PHY
!* 0. DECLARATIONS
!
!
-USE SHUMAN_PHY, ONLY: MZF_PHY, DZF_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZF_PHY, DZF_PHY
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
diff --git a/src/common/aux/mode_io_field_write.F90 b/src/common/aux/mode_io_field_write.F90
index aa739877b4202d245da26e70bc48c21513390a57..fe30f03ccb7931b0a87ce185671c3f86fdad4988 100644
--- a/src/common/aux/mode_io_field_write.F90
+++ b/src/common/aux/mode_io_field_write.F90
@@ -4,6 +4,8 @@ USE MODD_FIELD, ONLY: TFIELDMETADATA
CONTAINS
SUBROUTINE IO_FIELD_WRITE(TPFILE,TZFIELD,PFIELD)
!
+ ! THIS ROUTINE IS TEMPORARY AND IS ONLY NEEDED BY MODE_TURB_HOR* WHICH HAVE NOT YET
+ ! BEEN TRANSFORMED FOR GPU ADAPTATION
!* 0.1 Declarations of arguments
!
TYPE(TFILEDATA), INTENT(IN) :: TPFILE
@@ -13,18 +15,4 @@ SUBROUTINE IO_FIELD_WRITE(TPFILE,TZFIELD,PFIELD)
CALL ABORT
END SUBROUTINE IO_FIELD_WRITE
!
-SUBROUTINE IO_FIELD_WRITE_PHY(D,TPFILE,TZFIELD,PFIELD)
- USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(DIMPHYEX_t), INTENT(IN) :: D
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- TYPE(TFIELDMETADATA), INTENT(IN) :: TZFIELD
- REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PFIELD ! array containing the data field
- !
- CALL ABORT
-END SUBROUTINE IO_FIELD_WRITE_PHY
-!
END MODULE MODE_IO_FIELD_WRITE
-
diff --git a/src/common/aux/mode_io_field_write_phy.F90 b/src/common/aux/mode_io_field_write_phy.F90
new file mode 100644
index 0000000000000000000000000000000000000000..84a66b4746bef3c5a0f35e6b072108a19cd32538
--- /dev/null
+++ b/src/common/aux/mode_io_field_write_phy.F90
@@ -0,0 +1,30 @@
+MODULE MODE_IO_FIELD_WRITE_PHY
+USE MODD_IO, ONLY: TFILEDATA
+USE MODD_FIELD, ONLY: TFIELDMETADATA
+CONTAINS
+SUBROUTINE IO_FIELD_WRITE(TPFILE,TZFIELD,PFIELD)
+ !
+ !* 0.1 Declarations of arguments
+ !
+ TYPE(TFILEDATA), INTENT(IN) :: TPFILE
+ TYPE(TFIELDMETADATA), INTENT(IN) :: TZFIELD
+ REAL, DIMENSION(:,:,:), INTENT(IN) :: PFIELD ! array containing the data field
+ !
+ CALL ABORT
+END SUBROUTINE IO_FIELD_WRITE
+!
+SUBROUTINE IO_FIELD_WRITE_PHY(D,TPFILE,TZFIELD,PFIELD)
+ USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
+ !
+ !* 0.1 Declarations of arguments
+ !
+ TYPE(DIMPHYEX_t), INTENT(IN) :: D
+ TYPE(TFILEDATA), INTENT(IN) :: TPFILE
+ TYPE(TFIELDMETADATA), INTENT(IN) :: TZFIELD
+ REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PFIELD ! array containing the data field
+ !
+ CALL ABORT
+END SUBROUTINE IO_FIELD_WRITE_PHY
+!
+END MODULE MODE_IO_FIELD_WRITE_PHY
+
diff --git a/src/common/aux/shuman_phy.F90 b/src/common/aux/shuman_phy.F90
index 112da4864b78e285820b35b5d9d0776cbc387db8..ab02447333cc790df851676ff903c851228c3d98 100644
--- a/src/common/aux/shuman_phy.F90
+++ b/src/common/aux/shuman_phy.F90
@@ -1,4 +1,4 @@
-MODULE SHUMAN_PHY
+MODULE MODE_SHUMAN_PHY
IMPLICIT NONE
CONTAINS
! ###############################
@@ -1331,4 +1331,4 @@ CALL ABORT ! AROME SHOULD NOT CALLED HORIZONTAL FINITE DIFFERENCE
IF (LHOOK) CALL DR_HOOK('DYF',1,ZHOOK_HANDLE)
END SUBROUTINE DYF_PHY
!
-END MODULE SHUMAN_PHY
+END MODULE MODE_SHUMAN_PHY
diff --git a/src/common/micro/ice_adjust.F90 b/src/common/micro/ice_adjust.F90
index deec98d8ef92cda88971c6eb49b4b497ea1caac5..e2981ff0107fa5d87b62c408022608aa35e95ca8 100644
--- a/src/common/micro/ice_adjust.F90
+++ b/src/common/micro/ice_adjust.F90
@@ -118,7 +118,7 @@ USE MODD_TURB_n, ONLY: TURB_t
USE MODD_BUDGET, ONLY: TBUDGETDATA, TBUDGETCONF_t, NBUDGET_TH, NBUDGET_RV, NBUDGET_RC, NBUDGET_RI
USE MODD_RAIN_ICE_PARAM, ONLY : RAIN_ICE_PARAM_t
!
-USE MODE_BUDGET, ONLY: BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
+USE MODE_BUDGET_PHY, ONLY: BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
!
USE MODI_CONDENSATION
!
diff --git a/src/common/micro/mode_ice4_budgets.F90 b/src/common/micro/mode_ice4_budgets.F90
index 99cf34f3123f146272562530c881fec1efb33a29..03c550d31736c25e59c0ebe77506f0a9ad3c4d74 100644
--- a/src/common/micro/mode_ice4_budgets.F90
+++ b/src/common/micro/mode_ice4_budgets.F90
@@ -25,7 +25,7 @@ USE MODD_PARAM_ICE, ONLY: PARAM_ICE_t
!
USE MODD_FIELDS_ADDRESS ! index number for prognostic (theta and mixing ratios) and budgets
!
-USE MODE_BUDGET, ONLY: BUDGET_STORE_ADD_PHY, BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
+USE MODE_BUDGET_PHY, ONLY: BUDGET_STORE_ADD_PHY, BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
!
!
IMPLICIT NONE
diff --git a/src/common/micro/mode_ice4_sedimentation.F90 b/src/common/micro/mode_ice4_sedimentation.F90
index f6fc795edff66ad05a393e199992dfb0fb4955d5..c110a81f1cdf5e8616d6b5d1fd52d14ab7948ddc 100644
--- a/src/common/micro/mode_ice4_sedimentation.F90
+++ b/src/common/micro/mode_ice4_sedimentation.F90
@@ -42,7 +42,7 @@ USE MODD_RAIN_ICE_PARAM, ONLY: RAIN_ICE_PARAM_t
USE MODD_PARAM_ICE, ONLY: PARAM_ICE_t
!
USE MODE_MSG, ONLY: PRINT_MSG, NVERB_FATAL
-USE MODE_BUDGET, ONLY: BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
+USE MODE_BUDGET_PHY, ONLY: BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
!
USE MODE_ICE4_SEDIMENTATION_STAT, ONLY: ICE4_SEDIMENTATION_STAT
USE MODE_ICE4_SEDIMENTATION_SPLIT, ONLY: ICE4_SEDIMENTATION_SPLIT
diff --git a/src/common/micro/rain_ice.F90 b/src/common/micro/rain_ice.F90
index 2008e76f6b6552c424e91d78e78e8a310eed4a79..c17b5c3cd39bd3e426f2e41208e16ee11ec1550c 100644
--- a/src/common/micro/rain_ice.F90
+++ b/src/common/micro/rain_ice.F90
@@ -197,7 +197,7 @@ USE MODD_FIELDS_ADDRESS, ONLY : & ! common fields adress
& IRG, & ! Graupel
& IRH ! Hail
-USE MODE_BUDGET, ONLY: BUDGET_STORE_ADD_PHY, BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
+USE MODE_BUDGET_PHY, ONLY: BUDGET_STORE_ADD_PHY, BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
USE MODE_MSG, ONLY: PRINT_MSG, NVERB_FATAL
USE MODE_ICE4_RAINFR_VERT, ONLY: ICE4_RAINFR_VERT
diff --git a/src/common/turb/mode_prandtl.F90 b/src/common/turb/mode_prandtl.F90
index 9a0d77c6dd45218aecc31b142635e444bbb1d2c4..f911cba340f766ae9cb41483fdb0fbed6072bb3b 100644
--- a/src/common/turb/mode_prandtl.F90
+++ b/src/common/turb/mode_prandtl.F90
@@ -17,7 +17,7 @@ USE MODD_CTURB, ONLY : CSTURB_t
USE MODD_DIMPHYEX, ONLY : DIMPHYEX_t
USE MODD_PARAMETERS, ONLY : JPVEXT_TURB
!
-USE SHUMAN_PHY, ONLY: MZM_PHY,MZF_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZM_PHY,MZF_PHY
USE MODE_GRADIENT_M_PHY
IMPLICIT NONE
!----------------------------------------------------------------------------
@@ -154,7 +154,7 @@ USE MODD_PARAMETERS, ONLY: JPVEXT_TURB
USE MODE_EMOIST, ONLY: EMOIST
USE MODE_ETHETA, ONLY: ETHETA
USE MODE_GRADIENT_M_PHY, ONLY: GX_M_M_PHY, GY_M_M_PHY
-USE MODE_IO_FIELD_WRITE, ONLY: IO_FIELD_WRITE_PHY
+USE MODE_IO_FIELD_WRITE_PHY, ONLY: IO_FIELD_WRITE_PHY
!
IMPLICIT NONE
!
diff --git a/src/common/turb/mode_rmc01.F90 b/src/common/turb/mode_rmc01.F90
index 5a980a92be0abf8fd509eba3bf3a572b9df28d89..59f254a0d4af720a86f3692144d7fcad68e4d806 100644
--- a/src/common/turb/mode_rmc01.F90
+++ b/src/common/turb/mode_rmc01.F90
@@ -53,7 +53,7 @@ USE MODD_CTURB, ONLY: CSTURB_t
USE MODE_UPDATE_IIJU_PHY, ONLY: UPDATE_IIJU_PHY
USE MODE_SBL_PHY, ONLY: BUSINGER_PHIM, BUSINGER_PHIE
!
-USE SHUMAN_PHY, ONLY: MZF_PHY, MYF_PHY, MXF_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZF_PHY, MYF_PHY, MXF_PHY
!
IMPLICIT NONE
!
diff --git a/src/common/turb/mode_tke_eps_sources.F90 b/src/common/turb/mode_tke_eps_sources.F90
index f7f6a5082ae0c68b332cc8b4fed22347fbd2bfd1..333b8d6415f6542685e59989cd21175ee72ed1fd 100644
--- a/src/common/turb/mode_tke_eps_sources.F90
+++ b/src/common/turb/mode_tke_eps_sources.F90
@@ -128,7 +128,7 @@ CONTAINS
! ------------
!
USE PARKIND1, ONLY: JPRB
-USE SHUMAN_PHY, ONLY: MZM_PHY, MZF_PHY, DZF_PHY, DZM_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZM_PHY, MZF_PHY, DZF_PHY, DZM_PHY
USE YOMHOOK, ONLY: LHOOK, DR_HOOK
!
USE MODD_ARGSLIST_ll, ONLY: LIST_ll
@@ -142,8 +142,8 @@ USE MODD_LES, ONLY: TLES_t
USE MODD_PARAMETERS, ONLY: JPVEXT_TURB
USE MODD_TURB_n, ONLY: TURB_t
!
-USE MODE_BUDGET, ONLY: BUDGET_STORE_ADD_PHY, BUDGET_STORE_END_PHY, BUDGET_STORE_INIT_PHY
-USE MODE_IO_FIELD_WRITE, ONLY: IO_FIELD_WRITE_PHY
+USE MODE_BUDGET_PHY, ONLY: BUDGET_STORE_ADD_PHY, BUDGET_STORE_END_PHY, BUDGET_STORE_INIT_PHY
+USE MODE_IO_FIELD_WRITE_PHY, ONLY: IO_FIELD_WRITE_PHY
USE MODE_ll
!
USE MODI_GET_HALO
diff --git a/src/common/turb/mode_tridiag_thermo.F90 b/src/common/turb/mode_tridiag_thermo.F90
index fcef93a9789f5513e9b7a213a0c24796dcc516a8..23d959b9bdb5714bdc06ecfe4592456c2a0a5751 100644
--- a/src/common/turb/mode_tridiag_thermo.F90
+++ b/src/common/turb/mode_tridiag_thermo.F90
@@ -122,7 +122,7 @@ USE MODD_DIMPHYEX, ONLY : DIMPHYEX_t
USE MODD_PARAMETERS, ONLY : JPVEXT_TURB
!
USE MODI_SHUMAN, ONLY : MZM
-USE SHUMAN_PHY, ONLY: MZM_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZM_PHY
!
IMPLICIT NONE
!
diff --git a/src/common/turb/mode_turb_ver.F90 b/src/common/turb/mode_turb_ver.F90
index b014aa191d35db68f45b60b86ade811c2159ec3e..6847c703e900270816207ad37a90eec4e25386e7 100644
--- a/src/common/turb/mode_turb_ver.F90
+++ b/src/common/turb/mode_turb_ver.F90
@@ -225,7 +225,7 @@ USE MODD_TURB_n, ONLY: TURB_t
USE MODE_EMOIST, ONLY: EMOIST
USE MODE_ETHETA, ONLY: ETHETA
USE MODE_GRADIENT_M_PHY, ONLY: GZ_M_W_PHY
-USE MODE_IO_FIELD_WRITE, ONLY: IO_FIELD_WRITE_PHY
+USE MODE_IO_FIELD_WRITE_PHY, ONLY: IO_FIELD_WRITE_PHY
USE MODE_PRANDTL, ONLY: PSI_SV, PSI3, PHI3, PRANDTL
USE MODE_SBL_DEPTH, ONLY: SBL_DEPTH
USE MODE_TURB_VER_DYN_FLUX, ONLY: TURB_VER_DYN_FLUX
diff --git a/src/common/turb/mode_turb_ver_dyn_flux.F90 b/src/common/turb/mode_turb_ver_dyn_flux.F90
index 17f8d55009e8e4a0c251bc3f79d9b715a4ebcf6c..acda739bf432db377a546504ce3cb145282f4738 100644
--- a/src/common/turb/mode_turb_ver_dyn_flux.F90
+++ b/src/common/turb/mode_turb_ver_dyn_flux.F90
@@ -205,7 +205,7 @@ SUBROUTINE TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,TLES,KSV,O2D,OFLAT, &
! ------------
!
USE PARKIND1, ONLY: JPRB
-USE SHUMAN_PHY
+USE MODE_SHUMAN_PHY
USE YOMHOOK, ONLY: LHOOK, DR_HOOK
!
USE MODD_CST, ONLY: CST_t
@@ -221,7 +221,7 @@ USE MODE_GRADIENT_U_PHY, ONLY : GZ_U_UW_PHY, GX_U_M_PHY
USE MODE_GRADIENT_V_PHY, ONLY : GZ_V_VW_PHY, GY_V_M_PHY
USE MODE_GRADIENT_W_PHY, ONLY : GX_W_UW_PHY, GY_W_VW_PHY, GZ_W_M_PHY
USE MODE_GRADIENT_M_PHY, ONLY : GX_M_U_PHY, GY_M_V_PHY
-USE MODE_IO_FIELD_WRITE, only: IO_FIELD_WRITE_PHY
+USE MODE_IO_FIELD_WRITE_PHY, only: IO_FIELD_WRITE_PHY
USE MODE_ll
USE MODE_TRIDIAG_WIND, ONLY: TRIDIAG_WIND
!
diff --git a/src/common/turb/mode_turb_ver_sv_corr.F90 b/src/common/turb/mode_turb_ver_sv_corr.F90
index 2a1915b221dabe682ce5f2e6cf85e163a2941ff9..64acbc7f83418a05b71b816b7d6cf9cb86b452fa 100644
--- a/src/common/turb/mode_turb_ver_sv_corr.F90
+++ b/src/common/turb/mode_turb_ver_sv_corr.F90
@@ -62,7 +62,7 @@ USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
USE MODD_PARAMETERS, ONLY: JPVEXT_TURB
USE MODD_LES, ONLY: TLES_t
!
-USE SHUMAN_PHY, ONLY: MZF_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZF_PHY
USE MODE_GRADIENT_M_PHY, ONLY : GZ_M_W_PHY
USE MODE_EMOIST, ONLY: EMOIST
USE MODE_ETHETA, ONLY: ETHETA
diff --git a/src/common/turb/mode_turb_ver_sv_flux.F90 b/src/common/turb/mode_turb_ver_sv_flux.F90
index ab552090a450c7e48e04baf5d99f1b14dc4c2c6f..40a52e4ce77ea68335688a57bdc43324495a7207 100644
--- a/src/common/turb/mode_turb_ver_sv_flux.F90
+++ b/src/common/turb/mode_turb_ver_sv_flux.F90
@@ -210,7 +210,7 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, &
! ------------
!
USE PARKIND1, ONLY: JPRB
-USE SHUMAN_PHY, ONLY: DZM_PHY, MZM_PHY, MZF_PHY
+USE MODE_SHUMAN_PHY, ONLY: DZM_PHY, MZM_PHY, MZF_PHY
USE YOMHOOK, ONLY: LHOOK, DR_HOOK
!
USE MODD_CST, ONLY: CST_t
@@ -226,7 +226,7 @@ USE MODE_EMOIST, ONLY: EMOIST
USE MODE_ETHETA, ONLY: ETHETA
USE MODE_GRADIENT_W_PHY, ONLY: GZ_W_M_PHY
USE MODE_GRADIENT_M_PHY, ONLY: GZ_M_W_PHY
-USE MODE_IO_FIELD_WRITE, ONLY: IO_FIELD_WRITE_PHY
+USE MODE_IO_FIELD_WRITE_PHY, ONLY: IO_FIELD_WRITE_PHY
USE MODE_TRIDIAG, ONLY: TRIDIAG
!
USE MODI_LES_MEAN_SUBGRID_PHY
diff --git a/src/common/turb/mode_turb_ver_thermo_corr.F90 b/src/common/turb/mode_turb_ver_thermo_corr.F90
index 56c4e01d7c54ac8f2e02d8f7a272f356846a767a..8a53e4917a552ea1eda3b790e8ee7dcc90a5887a 100644
--- a/src/common/turb/mode_turb_ver_thermo_corr.F90
+++ b/src/common/turb/mode_turb_ver_thermo_corr.F90
@@ -205,7 +205,7 @@ SUBROUTINE TURB_VER_THERMO_CORR(D,CST,CSTURB,TURBN,TLES, &
! ------------
!
USE PARKIND1, ONLY: JPRB
-USE SHUMAN_PHY, ONLY: MZM_PHY, MZF_PHY, DZM_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZM_PHY, MZF_PHY, DZM_PHY
USE YOMHOOK, ONLY: LHOOK, DR_HOOK
!
USE MODD_CST, ONLY: CST_t
@@ -217,7 +217,7 @@ USE MODD_LES, ONLY: TLES_t
USE MODD_PARAMETERS, ONLY: JPVEXT_TURB
USE MODD_TURB_n, ONLY: TURB_t
!
-USE MODE_IO_FIELD_WRITE, ONLY: IO_FIELD_WRITE_PHY
+USE MODE_IO_FIELD_WRITE_PHY, ONLY: IO_FIELD_WRITE_PHY
USE MODE_PRANDTL
!
USE MODI_LES_MEAN_SUBGRID_PHY
diff --git a/src/common/turb/mode_turb_ver_thermo_flux.F90 b/src/common/turb/mode_turb_ver_thermo_flux.F90
index 618fb377a19214004ad7dcd906efc3e6fbad40b5..f041feff7ba15e721fe247aab52b85c4cd5c6d79 100644
--- a/src/common/turb/mode_turb_ver_thermo_flux.F90
+++ b/src/common/turb/mode_turb_ver_thermo_flux.F90
@@ -230,7 +230,7 @@ SUBROUTINE TURB_VER_THERMO_FLUX(D,CST,CSTURB,TURBN,TLES, &
! ------------
!
USE PARKIND1, ONLY: JPRB
-USE SHUMAN_PHY, ONLY: DZF_PHY, DZM_PHY, MXF_PHY, MYF_PHY, MZF_PHY, MZM_PHY
+USE MODE_SHUMAN_PHY, ONLY: DZF_PHY, DZM_PHY, MXF_PHY, MYF_PHY, MZF_PHY, MZM_PHY
USE YOMHOOK, ONLY: LHOOK, DR_HOOK
!
USE MODD_CST, ONLY: CST_t
@@ -243,7 +243,7 @@ USE MODD_PARAMETERS, ONLY: JPVEXT_TURB, JPHEXT, XUNDEF
USE MODD_TURB_n, ONLY: TURB_t
!
USE MODE_GRADIENT_W_PHY, ONLY: GZ_W_M_PHY
-USE MODE_IO_FIELD_WRITE, ONLY: IO_FIELD_WRITE_PHY
+USE MODE_IO_FIELD_WRITE_PHY, ONLY: IO_FIELD_WRITE_PHY
USE MODE_PRANDTL
USE MODE_TM06_H, ONLY: TM06_H
USE MODE_TRIDIAG_THERMO, ONLY: TRIDIAG_THERMO
diff --git a/src/common/turb/shallow_mf.F90 b/src/common/turb/shallow_mf.F90
index 898883076890bc04c938b2339d1c956b653fff5c..0eaa8b4056af502c876fde0cc385a7c7f277c8ec 100644
--- a/src/common/turb/shallow_mf.F90
+++ b/src/common/turb/shallow_mf.F90
@@ -73,7 +73,7 @@
! ------------
!
USE PARKIND1, ONLY: JPRB
-USE SHUMAN_PHY, ONLY: MXM_PHY, MYM_PHY
+USE MODE_SHUMAN_PHY, ONLY: MXM_PHY, MYM_PHY
USE YOMHOOK, ONLY: LHOOK, DR_HOOK
!
USE MODD_BUDGET, ONLY: TBUDGETCONF_t, TBUDGETDATA, NBUDGET_U, NBUDGET_V, &
@@ -86,7 +86,7 @@ USE MODD_TURB_n, ONLY: TURB_t
USE MODD_CTURB, ONLY: CSTURB_t
USE MODD_PARAMETERS, ONLY: JPSVMAX
!
-USE MODE_BUDGET, ONLY: BUDGET_STORE_ADD_PHY
+USE MODE_BUDGET_PHY, ONLY: BUDGET_STORE_ADD_PHY
USE MODE_COMPUTE_MF_CLOUD, ONLY: COMPUTE_MF_CLOUD
USE MODE_COMPUTE_UPDRAFT, ONLY: COMPUTE_UPDRAFT
USE MODE_COMPUTE_UPDRAFT_RAHA, ONLY: COMPUTE_UPDRAFT_RAHA
diff --git a/src/common/turb/turb.F90 b/src/common/turb/turb.F90
index 0cfafc17ae355cf3decb1037417d15745a7df402..3cf7afa999a48d7087cbd67d7c483183f9ed3804 100644
--- a/src/common/turb/turb.F90
+++ b/src/common/turb/turb.F90
@@ -240,7 +240,7 @@
! ------------
!
USE PARKIND1, ONLY: JPRB
-USE SHUMAN_PHY, ONLY: MZF_PHY,MXF_PHY,MYF_PHY
+USE MODE_SHUMAN_PHY, ONLY: MZF_PHY,MXF_PHY,MYF_PHY
USE YOMHOOK , ONLY: LHOOK, DR_HOOK
!
USE MODD_BUDGET, ONLY: NBUDGET_U, NBUDGET_V, NBUDGET_W, NBUDGET_TH, NBUDGET_RV, NBUDGET_RC, &
@@ -256,7 +256,7 @@ USE MODD_PARAMETERS, ONLY: JPVEXT_TURB, XUNDEF
USE MODD_TURB_n, ONLY: TURB_t
!
USE MODE_BL89, ONLY: BL89
-USE MODE_BUDGET, ONLY: BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
+USE MODE_BUDGET_PHY, ONLY: BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
USE MODE_EMOIST, ONLY: EMOIST
USE MODE_ETHETA, ONLY: ETHETA
USE MODE_GRADIENT_U_PHY, ONLY: GZ_U_UW_PHY
@@ -264,7 +264,7 @@ USE MODE_GRADIENT_V_PHY, ONLY: GZ_V_VW_PHY
USE MODE_GRADIENT_W_PHY, ONLY: GZ_W_M_PHY
USE MODE_GRADIENT_M_PHY, ONLY: GZ_M_W_PHY
USE MODE_IBM_MIXINGLENGTH, ONLY: IBM_MIXINGLENGTH
-USE MODE_IO_FIELD_WRITE, ONLY: IO_FIELD_WRITE_PHY
+USE MODE_IO_FIELD_WRITE_PHY, ONLY: IO_FIELD_WRITE_PHY
USE MODE_RMC01, ONLY: RMC01
USE MODE_ROTATE_WIND, ONLY: ROTATE_WIND, UPDATE_ROTATE_WIND
USE MODE_SBL_PHY, ONLY: LMO
diff --git a/src/mesonh/aux/mode_budget.f90 b/src/mesonh/aux/mode_budget.f90
deleted file mode 100644
index 89aad3a6eb956fdeddfcf552472c115ccc794a02..0000000000000000000000000000000000000000
--- a/src/mesonh/aux/mode_budget.f90
+++ /dev/null
@@ -1,349 +0,0 @@
-!MNH_LIC Copyright 1994-2021 CNRS, Meteo-France and Universite Paul Sabatier
-!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
-!MNH_LIC for details. version 1.
-!-----------------------------------------------------------------
-! Modifications
-! P. Wautelet 28/01/2020: new subroutines: Budget_store_init, Budget_store_end and Budget_source_id_find in new module mode_budget
-! P. Wautelet 17/08/2020: treat LES budgets correctly
-! P. Wautelet 05/03/2021: measure cpu_time for budgets
-!-----------------------------------------------------------------
-
-!#################
-module mode_budget
-!#################
-
-use modd_budget, only: cbutype, nbutime, tbudgetdata, xtime_bu, xtime_bu_process
-use modd_les_budget, only: xtime_les_bu, xtime_les_bu_process
-
-use modi_cart_compress, only: Cart_compress
-use modi_mask_compress, only: Mask_compress
-use modi_second_mnh, only: Second_mnh
-
-use mode_msg
-
-implicit none
-
-private
-
-public :: Budget_store_init,Budget_store_init_phy
-public :: Budget_store_end, Budget_store_end_phy
-public :: Budget_store_add, Budget_store_add_phy
-
-real :: ztime1, ztime2
-
-contains
-
-subroutine Budget_store_init_phy(D, tpbudget, hsource, pvars)
- USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
- TYPE(DIMPHYEX_t), INTENT(IN) :: D
- type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
- character(len=*), intent(in) :: hsource ! Name of the source term
- real, dimension(D%NIT,D%NJT,D%NKT), intent(in) :: pvars ! Current value to be stored
-!
- call Budget_store_init(tpbudget, hsource, pvars)
-!
-end subroutine Budget_store_init_phy
-!
-subroutine Budget_store_end_phy(D, tpbudget, hsource, pvars)
- USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
- TYPE(DIMPHYEX_t), INTENT(IN) :: D
- type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
- character(len=*), intent(in) :: hsource ! Name of the source term
- real, dimension(D%NIT,D%NJT,D%NKT), intent(in) :: pvars ! Current value to be stored
-!
- call Budget_store_end(tpbudget, hsource, pvars)
-!
-end subroutine Budget_store_end_phy
-!
-subroutine Budget_store_add_phy(D, tpbudget, hsource, pvars)
- USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
- TYPE(DIMPHYEX_t), INTENT(IN) :: D
- type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
- character(len=*), intent(in) :: hsource ! Name of the source term
- real, dimension(D%NIT,D%NJT,D%NKT), intent(in) :: pvars ! Current value to be stored
-!
- call Budget_store_add(tpbudget, hsource, pvars)
-!
-end subroutine Budget_store_add_phy
-!
-subroutine Budget_store_init( tpbudget, hsource, pvars )
- use modd_les, only: lles_call
-
- type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
- character(len=*), intent(in) :: hsource ! Name of the source term
- real, dimension(:,:,:), intent(in) :: pvars ! Current value to be stored
-
- integer :: iid ! Reference number of the current source term
-
- call Print_msg( NVERB_DEBUG, 'BUD', 'Budget_store_init', trim( tpbudget%cname )//':'//trim( hsource ) )
-
- if ( lles_call ) then
- call Second_mnh( ztime1 )
-
- if ( allocated( tpbudget%xtmplesstore ) ) then
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_init', 'xtmplesstore already allocated' )
- else
- allocate( tpbudget%xtmplesstore( Size( pvars, 1 ), Size( pvars, 2 ), Size ( pvars, 3 ) ) )
- end if
- tpbudget%xtmplesstore(:, :, :) = pvars(:, :, :)
-
- tpbudget%clessource = hsource
-
- call Second_mnh( ztime2 )
- xtime_les_bu = xtime_les_bu + ztime2 - ztime1
- xtime_les_bu_process = xtime_les_bu_process + ztime2 - ztime1
- end if
-
- ! Nothing else to do if budgets are not enabled
- if ( .not. tpbudget%lenabled ) return
-
- call Second_mnh( ztime1 )
-
- call Budget_source_id_find( tpbudget, hsource, iid )
-
- if ( tpbudget%ntmpstoresource /= 0 ) then
- cmnhmsg(1) = 'ntmpstoresource already set (previous call to '//'Budget_store_end missing?)'
- cmnhmsg(2) = 'Set for: ' // Trim( tpbudget%cname ) // ':' // Trim( tpbudget%tsources(tpbudget%ntmpstoresource)%cmnhname )
- cmnhmsg(3) = 'Working on: ' // Trim( tpbudget%cname ) // ':' // Trim( hsource )
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_init' )
- end if
-
- if ( tpbudget%tsources(iid)%ldonotinit ) then
- ! If ldonotinit is set, this subroutine should not be called
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_init', 'should not be called for ' &
- //trim( tpbudget%cname )//':'//trim( hsource ) )
- return
- end if
-
- if ( tpbudget%tsources(iid)%lenabled ) then
- if ( tpbudget%ntmpstoresource /= 0 ) then
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_init', 'xtmpstore already used by ' &
- //trim( tpbudget%tsources(tpbudget%ntmpstoresource)%cmnhname ) )
- return
- end if
-
- tpbudget%ntmpstoresource = iid
-
- !Store data into the budget temporary array
- !This value will be subtracted from the next one (in Budget_store_end) to get the evolution of the array between the 2 calls
- if ( cbutype == 'CART' ) then
- tpbudget%xtmpstore(:, :, :) = Cart_compress( pvars(:, :, :) )
- else if ( cbutype == 'MASK' ) then
- tpbudget%xtmpstore(:, nbutime, :) = Mask_compress( pvars(:, :, :) )
- else
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_init', 'unknown cbutype: '//trim( cbutype ) )
- end if
- end if
-
- call Second_mnh( ztime2 )
- xtime_bu = xtime_bu + ztime2 - ztime1
- xtime_bu_process = xtime_bu_process + ztime2 - ztime1
-
- end subroutine Budget_store_init
-
-
-subroutine Budget_store_end( tpbudget, hsource, pvars )
- use modd_les, only: lles_call
-
- use modi_les_budget, only: Les_budget
-
- type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
- character(len=*), intent(in) :: hsource ! Name of the source term
- real, dimension(:,:,:), intent(in) :: pvars ! Current value to be stored
-
- integer :: iid ! Reference number of the current source term
- integer :: igroup ! Number of the group where to store the source term
- real, dimension(:,:,:), allocatable :: zvars_add
-
- call Print_msg( NVERB_DEBUG, 'BUD', 'Budget_store_end', trim( tpbudget%cname )//':'//trim( hsource ) )
-
- if ( lles_call ) then
- if ( hsource /= tpbudget%clessource ) &
- call Print_msg( NVERB_FATAL, 'BUD', 'Budget_store_end', 'hsource not the same as in Budget_store_init (' &
- // Trim( hsource ) // ' / ' // Trim( tpbudget%clessource ) // ')' )
-
- tpbudget%clessource = 'reset'
-
- if ( allocated( tpbudget%xtmplesstore ) ) then
- ! Do the call to Les_budget with oadd=.true.
- ! This is necessary when the call to Budget_store_init was done with pvars not strictly
- ! equal to the source term
- Allocate( zvars_add( Size( pvars, 1 ), Size( pvars, 2 ), Size ( pvars, 3 ) ) )
- zvars_add(:, :, :) = pvars(:, :, :) - tpbudget%xtmplesstore(:, :, :)
- call Les_budget( zvars_add, tpbudget%nid, hsource, oadd = .true. )
- Deallocate( zvars_add )
- Deallocate( tpbudget%xtmplesstore )
- else
- call Les_budget( pvars, tpbudget%nid, hsource, oadd = .false. )
- end if
- end if
-
- ! Nothing to do if budgets are not enabled
- if ( .not. tpbudget%lenabled ) return
-
- call Second_mnh( ztime1 )
-
- call Budget_source_id_find( tpbudget, hsource, iid )
-
- if ( tpbudget%tsources(iid)%lenabled ) then
- if ( iid /= tpbudget%ntmpstoresource .and. .not.tpbudget%tsources(iid)%ldonotinit ) then
- if ( tpbudget%ntmpstoresource == 0 ) then
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_end', 'ntmpstoresource not set for ' &
- //trim( tpbudget%tsources(iid)%cmnhname ) )
- else
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_end', 'xtmpstore used by an other source: ' &
- //trim( tpbudget%tsources(tpbudget%ntmpstoresource)%cmnhname )//', expected: ' &
- //trim( tpbudget%tsources(iid)%cmnhname ) )
- end if
- end if
-
- !Store data into the budget array
- !The values are computed by the difference between the values stored in the temporary array (filled in Budget_store_init)
- !and the current values added to the already stored ones.
- !Except if ldonotinit is true. In that case, overwrite the array.
- igroup = tpbudget%tsources(iid)%ngroup
- if ( cbutype == 'CART' ) then
- if ( tpbudget%tsources(iid)%ldonotinit ) then
- if ( tpbudget%tsources(iid)%loverwrite ) then
- tpbudget%tgroups(igroup)%xdata(:, :, :) = Cart_compress( pvars(:, :, :) )
- else
- tpbudget%tgroups(igroup)%xdata(:, :, :) = tpbudget%tgroups(igroup)%xdata(:, :, :) &
- + Cart_compress( pvars(:, :, :) )
- end if
- else
- if ( tpbudget%tsources(iid)%loverwrite ) then
- tpbudget%tgroups(igroup)%xdata(:, :, :) = Cart_compress( pvars(:, :, :) ) &
- - tpbudget%xtmpstore(:, :, :)
- else
- tpbudget%tgroups(igroup)%xdata(:, :, :) = tpbudget%tgroups(igroup)%xdata(:, :, :) &
- + Cart_compress( pvars(:, :, :) ) &
- - tpbudget%xtmpstore(:, :, :)
- end if
- end if
- else if ( cbutype == 'MASK' ) then
- if ( tpbudget%tsources(iid)%ldonotinit ) then
- if ( tpbudget%tsources(iid)%loverwrite ) then
- tpbudget%tgroups(igroup)%xdata(:, nbutime, :) = Mask_compress( pvars(:, :, :) )
- else
- tpbudget%tgroups(igroup)%xdata(:, nbutime, :) = tpbudget%tgroups(igroup)%xdata(:, nbutime, :) &
- + Mask_compress( pvars(:, :, :) )
- end if
- else
- if ( tpbudget%tsources(iid)%loverwrite ) then
- tpbudget%tgroups(igroup)%xdata(:, nbutime, :) = Mask_compress( pvars(:, :, :) ) &
- - tpbudget%xtmpstore(:, nbutime, :)
- else
- tpbudget%tgroups(igroup)%xdata(:, nbutime, :) = tpbudget%tgroups(igroup)%xdata(:, nbutime, :) &
- + Mask_compress( pvars(:, :, :) ) &
- - tpbudget%xtmpstore(:, nbutime, :)
- end if
- end if
- else
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_end', 'unknown cbutype: '//trim( cbutype ) )
- end if
-
- ! Release the budget temporary array
- tpbudget%ntmpstoresource = 0
- end if
-
- call Second_mnh( ztime2 )
- xtime_bu = xtime_bu + ztime2 - ztime1
- xtime_bu_process = xtime_bu_process + ztime2 - ztime1
-
-end subroutine Budget_store_end
-
-
-subroutine Budget_store_add( tpbudget, hsource, pvars )
- use modd_les, only: lles_call
-
- use modi_les_budget, only: Les_budget
-
- type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
- character(len=*), intent(in) :: hsource ! Name of the source term
- real, dimension(:,:,:), intent(in) :: pvars ! Current value to be stored
-
- integer :: iid ! Reference number of the current source term
- integer :: igroup ! Number of the group where to store the source term
-
- call Print_msg( NVERB_DEBUG, 'BUD', 'Budget_store_add', trim( tpbudget%cname )//':'//trim( hsource ) )
-
- if ( tpbudget%ntmpstoresource /= 0 ) &
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_add', 'inside a Budget_store_init/Budget_store_end zone' )
-
- if ( lles_call ) call Les_budget( pvars, tpbudget%nid, hsource, oadd = .true. )
-
- ! Nothing to do if budgets are not enabled
- if ( .not. tpbudget%lenabled ) return
-
- call Second_mnh( ztime1 )
-
- call Budget_source_id_find( tpbudget, hsource, iid )
-
- if ( tpbudget%tsources(iid)%lenabled ) then
- if ( tpbudget%tsources(iid)%loverwrite ) &
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_add', 'loverwrite=.true. is not allowed' )
-
- !Store data into the budget array
- igroup = tpbudget%tsources(iid)%ngroup
- if ( cbutype == 'CART' ) then
- tpbudget%tgroups(igroup)%xdata(:, :, :) = tpbudget%tgroups(igroup)%xdata(:, :, :) &
- + Cart_compress( pvars(:, :, :) )
- else if ( cbutype == 'MASK' ) then
- tpbudget%tgroups(igroup)%xdata(:, nbutime, :) = tpbudget%tgroups(igroup)%xdata(:, nbutime, :) &
- + Mask_compress( pvars(:, :, :) )
- else
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_store_add', 'unknown cbutype: '//trim( cbutype ) )
- end if
- end if
-
- call Second_mnh( ztime2 )
- xtime_bu = xtime_bu + ztime2 - ztime1
- xtime_bu_process = xtime_bu_process + ztime2 - ztime1
-
-end subroutine Budget_store_add
-
-
-subroutine Budget_source_id_find( tpbudget, hsource, kid )
- type(tbudgetdata), intent(in) :: tpbudget ! Budget datastructure
- character(len=*), intent(in) :: hsource ! Name of the source term
- integer, intent(out) :: kid ! Reference number of the current source term
-
- integer :: iid
- integer :: ji
-
- call Print_msg( NVERB_DEBUG, 'BUD', 'Budget_source_id_find', trim( tpbudget%cname )//':'//trim( hsource ) )
-
- iid = 0
- do ji = 1, tpbudget%nsources
- if ( trim( hsource ) == trim( tpbudget%tsources(ji)%cmnhname ) ) then
- iid = ji
- exit
- end if
- end do
-
- if ( iid > 0 ) then
- call Print_msg( NVERB_DEBUG, 'BUD', 'Budget_source_id_find', trim( tpbudget%cname )//':'//trim( hsource )//' found' )
- else
- !Search also in the non-available source term list
- do ji = tpbudget%nsources + 1, tpbudget%nsourcesmax
- if ( trim( hsource ) == trim( tpbudget%tsources(ji)%cmnhname ) ) then
- iid = ji
- exit
- end if
- end do
-
- if ( iid == 0 ) then
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_source_id_find', trim( tpbudget%cname )//':'//trim( hsource )//' not found' )
- else
- cmnhmsg(1) = Trim( tpbudget%cname ) // ':' // Trim( hsource ) // ' found'
- cmnhmsg(2) = 'in non-available source term list.'
- cmnhmsg(3) = 'Check availability condition in Ini_budget.'
- call Print_msg( NVERB_ERROR, 'BUD', 'Budget_source_id_find' )
- end if
- end if
-
- kid = iid
-end subroutine Budget_source_id_find
-
-end module mode_budget
diff --git a/src/mesonh/aux/mode_budget_phy.f90 b/src/mesonh/aux/mode_budget_phy.f90
new file mode 100644
index 0000000000000000000000000000000000000000..04a610d2066ab9f38f9879b3edc31c62c38bfdb7
--- /dev/null
+++ b/src/mesonh/aux/mode_budget_phy.f90
@@ -0,0 +1,62 @@
+!MNH_LIC Copyright 1994-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
+!MNH_LIC for details. version 1.
+!-----------------------------------------------------------------
+! Modifications
+! P. Wautelet 28/01/2020: new SUBROUTINEs: Budget_store_init, Budget_store_end and Budget_source_id_find in new module mode_budget
+! P. Wautelet 17/08/2020: treat LES budgets correctly
+! P. Wautelet 05/03/2021: measure cpu_time for budgets
+!-----------------------------------------------------------------
+
+!#################
+MODULE MODE_BUDGET_PHY
+!#################
+
+USE MODD_BUDGET, ONLY: TBUDGETDATA
+USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
+
+IMPLICIT NONE
+
+PRIVATE
+
+PUBLIC :: Budget_store_init_phy
+PUBLIC :: Budget_store_end_phy
+PUBLIC :: Budget_store_add_phy
+
+CONTAINS
+
+SUBROUTINE Budget_store_init_phy(D, tpbudget, hsource, pvars)
+ USE MODE_BUDGET, ONLY: BUDGET_STORE_INIT
+ TYPE(DIMPHYEX_t), INTENT(IN) :: D
+ type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
+ character(len=*), intent(in) :: hsource ! Name of the source term
+ real, dimension(D%NIT,D%NJT,D%NKT), intent(in) :: pvars ! Current value to be stored
+!
+ CALL Budget_store_init(tpbudget, hsource, pvars)
+!
+END SUBROUTINE Budget_store_init_phy
+!
+SUBROUTINE Budget_store_end_phy(D, tpbudget, hsource, pvars)
+ USE MODE_BUDGET, ONLY: BUDGET_STORE_END
+ TYPE(DIMPHYEX_t), INTENT(IN) :: D
+ type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
+ character(len=*), intent(in) :: hsource ! Name of the source term
+ real, dimension(D%NIT,D%NJT,D%NKT), intent(in) :: pvars ! Current value to be stored
+!
+ CALL Budget_store_end(tpbudget, hsource, pvars)
+!
+END SUBROUTINE Budget_store_end_phy
+!
+SUBROUTINE Budget_store_add_phy(D, tpbudget, hsource, pvars)
+ USE MODE_BUDGET, ONLY: BUDGET_STORE_ADD
+ TYPE(DIMPHYEX_t), INTENT(IN) :: D
+ type(tbudgetdata), intent(inout) :: tpbudget ! Budget datastructure
+ character(len=*), intent(in) :: hsource ! Name of the source term
+ real, dimension(D%NIT,D%NJT,D%NKT), intent(in) :: pvars ! Current value to be stored
+!
+ CALL Budget_store_add(tpbudget, hsource, pvars)
+!
+END SUBROUTINE Budget_store_add_phy
+!
+END MODULE MODE_BUDGET_PHY
diff --git a/src/mesonh/aux/mode_io_field_write.f90 b/src/mesonh/aux/mode_io_field_write.f90
deleted file mode 100644
index e513093ca60097b15caa167b738b12576ebc5a8e..0000000000000000000000000000000000000000
--- a/src/mesonh/aux/mode_io_field_write.f90
+++ /dev/null
@@ -1,4400 +0,0 @@
-!MNH_LIC Copyright 1994-2021 CNRS, Meteo-France and Universite Paul Sabatier
-!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
-!MNH_LIC for details. version 1.
-!-----------------------------------------------------------------
-! Modifications:
-! J. Escobar 15/09/2015: WENO5 & JPHEXT <> 1
-! P. Wautelet 05/2016-04/2018: new data structures and calls for I/O
-! P. Wautelet 10/01/2019: do not write scalars in Z-split files
-! P. Wautelet 10/01/2019: write header also for Z-split files
-! P. Wautelet 05/03/2019: rename IO subroutines and modules
-! P. Wautelet 10/04/2019: replace ABORT and STOP calls by Print_msg
-! P. Wautelet 12/04/2019: added pointers for C1D, L1D, N1D, X5D and X6D structures in TFIELDDATA
-! P. Wautelet 12/04/2019: use MNHTIME for time measurement variables
-! P. Wautelet 12/07/2019: add support for 1D array of dates
-! J. Escobar 11/02/2020: for GA & // IO, add sync, & mpi_allreduce for error handling in // IO
-! P. Wautelet 22/09/2020: use ldimreduced to allow reduction in the number of dimensions of fields (used by 2D simulations)
-! P. Wautelet 30/09/2020: add IO_Field_write_box_byfield_X3 and IO_Field_write_error_check subroutines
-! P. Wautelet 04/12/2020: add IO_Field_create and IO_Ndimlist_reduce subroutines
-! P. Wautelet 07/12/2020: add support for partial write of fields (optional argument: koffset, not all subroutines, no LFI spport)
-! P. Wautelet 14/01/2021: add IO_Field_write_byname_N4 and IO_Field_write_byfield_N4 subroutines
-!-----------------------------------------------------------------
-
-#define MNH_SCALARS_IN_SPLITFILES 0
-
-MODULE MODE_IO_FIELD_WRITE
-
- use modd_field, only: tfieldlist, tfieldmetadata, tfieldmetadata_base, &
- TYPECHAR, TYPEDATE, TYPEINT, TYPELOG, TYPEREAL
- USE MODD_IO, ONLY: TFILEDATA, TOUTBAK
- USE MODD_MPIF
- use modd_parameters, only: NMNHNAMELGTMAX
- use modd_precision, only: MNHINT_MPI, MNHREAL_MPI, MNHTIME
-
- use mode_field, only: Find_field_id_from_mnhname
- USE MODE_IO_WRITE_LFI
-#ifdef MNH_IOCDF4
- USE MODE_IO_WRITE_NC4
-#endif
- use mode_msg
-
- IMPLICIT NONE
-
- PRIVATE
-
- public :: IO_Field_write, IO_Field_write_box, IO_Field_write_lb
- public :: IO_Field_write_phy
- public :: IO_Fieldlist_write, IO_Field_user_write
- public :: IO_Header_write, IO_Field_create
-
- INTERFACE IO_Field_write
- MODULE PROCEDURE IO_Field_write_byname_X0, IO_Field_write_byname_X1, &
- IO_Field_write_byname_X2, IO_Field_write_byname_X3, &
- IO_Field_write_byname_X4, IO_Field_write_byname_X5, &
- IO_Field_write_byname_X6, &
- IO_Field_write_byname_N0, IO_Field_write_byname_N1, &
- IO_Field_write_byname_N2, IO_Field_write_byname_N3, &
- IO_Field_write_byname_N4, &
- IO_Field_write_byname_L0, IO_Field_write_byname_L1, &
- IO_Field_write_byname_C0, IO_Field_write_byname_C1, &
- IO_Field_write_byname_T0, IO_Field_write_byname_T1, &
- IO_Field_write_byfield_X0,IO_Field_write_byfield_X1, &
- IO_Field_write_byfield_X2,IO_Field_write_byfield_X3, &
- IO_Field_write_byfield_X4,IO_Field_write_byfield_X5, &
- IO_Field_write_byfield_X6, &
- IO_Field_write_byfield_N0,IO_Field_write_byfield_N1, &
- IO_Field_write_byfield_N2,IO_Field_write_byfield_N3, &
- IO_Field_write_byfield_N4, &
- IO_Field_write_byfield_L0,IO_Field_write_byfield_L1, &
- IO_Field_write_byfield_C0,IO_Field_write_byfield_C1, &
- IO_Field_write_byfield_T0,IO_Field_write_byfield_T1
- END INTERFACE
-
- INTERFACE IO_Field_write_phy
- MODULE PROCEDURE IO_Field_write_phy_byfield_X2, IO_Field_write_phy_byfield_X1
- END INTERFACE
-
- INTERFACE IO_Field_write_box
- MODULE PROCEDURE IO_Field_write_box_byfield_X2, IO_Field_write_box_byfield_X3, &
- IO_Field_write_box_byfield_X4, IO_Field_write_box_byfield_X5
- END INTERFACE
-
- INTERFACE IO_Field_write_lb
- MODULE PROCEDURE IO_Field_write_byname_lb, IO_Field_write_byfield_lb
- END INTERFACE
-
-CONTAINS
-
- SUBROUTINE IO_Field_metadata_check(TPFIELD,KTYPE,KDIMS,HCALLER)
- CLASS(tfieldmetadata_base), INTENT(IN) :: TPFIELD ! Field to check
- INTEGER, INTENT(IN) :: KTYPE ! Expected datatype
- INTEGER, INTENT(IN) :: KDIMS ! Expected number of dimensions
- CHARACTER(LEN=*), INTENT(IN) :: HCALLER ! name of the calling subroutine
- !
- CHARACTER(LEN=2) :: YDIMOK,YDIMKO
- CHARACTER(LEN=8) :: YTYPEOK,YTYPEKO
- !
- IF (TPFIELD%NGRID<0 .OR. TPFIELD%NGRID>8) THEN
- CALL PRINT_MSG(NVERB_WARNING,'IO',HCALLER,'TPFIELD%NGRID is invalid for '//TRIM(TPFIELD%CMNHNAME))
- END IF
- IF (TPFIELD%NTYPE/=KTYPE) THEN
- CALL TYPE_WRITE(KTYPE,YTYPEOK)
- CALL TYPE_WRITE(TPFIELD%NTYPE,YTYPEKO)
- CALL PRINT_MSG(NVERB_WARNING,'IO',HCALLER,&
- 'TPFIELD%NTYPE should be '//YTYPEOK//' instead of '//YTYPEKO//' for '//TRIM(TPFIELD%CMNHNAME))
- END IF
- IF (TPFIELD%NDIMS/=KDIMS) THEN
- WRITE (YDIMOK,'(I2)') KDIMS
- WRITE (YDIMKO,'(I2)') TPFIELD%NDIMS
- CALL PRINT_MSG(NVERB_WARNING,'IO',HCALLER,&
- 'TPFIELD%NDIMS should be '//YDIMOK//' instead of '//YDIMKO//' for '//TRIM(TPFIELD%CMNHNAME))
- END IF
- !
- CONTAINS
- SUBROUTINE TYPE_WRITE(KTYPEINT,HTYPE)
- INTEGER, INTENT(IN) :: KTYPEINT
- CHARACTER(LEN=8),INTENT(OUT) :: HTYPE
- !
- SELECT CASE(KTYPEINT)
- CASE(TYPEINT)
- HTYPE = 'TYPEINT'
- CASE(TYPELOG)
- HTYPE = 'TYPELOG'
- CASE(TYPEREAL)
- HTYPE = 'TYPEREAL'
- CASE(TYPECHAR)
- HTYPE = 'TYPECHAR'
- CASE(TYPEDATE)
- HTYPE = 'TYPEDATE'
- CASE DEFAULT
- HTYPE = 'UNKNOWN'
- END SELECT
- !
- END SUBROUTINE TYPE_WRITE
- END SUBROUTINE IO_Field_metadata_check
-
-
- SUBROUTINE IO_File_write_check(TPFILE,HSUBR,KRESP)
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HSUBR
- INTEGER, INTENT(OUT) :: KRESP
- !
- KRESP = 0
- !
- !Check if file is opened
- IF (.NOT.TPFILE%LOPENED) THEN
- CALL PRINT_MSG(NVERB_ERROR,'IO',HSUBR,TRIM(TPFILE%CNAME)//' is not opened')
- KRESP = -201
- RETURN
- END IF
- !
- !Check if file is in the right opening mode
- IF (TPFILE%CMODE/='WRITE') THEN
- CALL PRINT_MSG(NVERB_WARNING,'IO',HSUBR,&
- TRIM(TPFILE%CNAME)//': writing in a file opened in '//TRIM(TPFILE%CMODE)//' mode')
- END IF
- !
- !Check fileformat
- IF (TPFILE%CFORMAT/='NETCDF4' .AND. TPFILE%CFORMAT/='LFI' .AND. TPFILE%CFORMAT/='LFICDF4') THEN
- CALL PRINT_MSG(NVERB_FATAL,'IO',HSUBR,&
- TRIM(TPFILE%CNAME)//': invalid fileformat ('//TRIM(TPFILE%CFORMAT)//')')
- KRESP = -202
- RETURN
- END IF
- !
- END SUBROUTINE IO_File_write_check
-
-
- SUBROUTINE IO_Format_write_select(TPFILE,OLFI,ONC4)
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! File structure
- LOGICAL, INTENT(OUT) :: OLFI ! Write in LFI format?
- LOGICAL, INTENT(OUT) :: ONC4 ! Write in netCDF format?
-
- OLFI = .FALSE.
- ONC4 = .FALSE.
- IF (TPFILE%CFORMAT=='LFI' .OR. TPFILE%CFORMAT=='LFICDF4') OLFI = .TRUE.
-#ifdef MNH_IOCDF4
- IF (TPFILE%CFORMAT=='NETCDF4' .OR. TPFILE%CFORMAT=='LFICDF4') ONC4 = .TRUE.
-#endif
- END SUBROUTINE IO_Format_write_select
-
-
- subroutine IO_Field_write_error_check( tpfile, tpfield, hsubr, kresp_in, kresp_lfi, kresp_nc4, kresp_out )
- use modd_io, only: gsmonoproc
-
- type(tfiledata), intent(in) :: tpfile
- class(tfieldmetadata_base), intent(in) :: tpfield
- character(len=*), intent(in) :: hsubr
- integer, intent(in) :: kresp_in
- integer, intent(in) :: kresp_lfi
- integer, intent(in) :: kresp_nc4
- integer, intent(out) :: kresp_out
-
- character(len=:), allocatable :: ymsg
- character(len=6) :: yresp
- integer :: ierr_mpi
- integer, dimension(3) :: iresps
-
- kresp_out = 0
-
- iresps(1) = kresp_in
- iresps(2) = kresp_lfi
- iresps(3) = kresp_nc4
-
- if ( .not. gsmonoproc ) call MPI_BCAST( iresps, 3, MNHINT_MPI, tpfile%nmaster_rank - 1, tpfile%nmpicomm, ierr_mpi )
-
- if ( iresps(1) /= 0 ) then
- write(yresp, '( i6 )') iresps(1)
- ymsg = Trim( tpfile%cname ) // ': resp=' // yresp // ' when writing ' // Trim( tpfield%cmnhname )
- call Print_msg( NVERB_ERROR, 'IO', hsubr, ymsg )
- kresp_out = iresps(1)
- end if
-
-#ifdef MNH_IOLFI
- if ( iresps(2) /= 0 ) then
- write(yresp, '( i6 )') iresps(2)
- ymsg = Trim( tpfile%cname ) // ': LFI: resp=' // yresp // ' when writing ' // Trim( tpfield%cmnhname )
- call Print_msg( NVERB_ERROR, 'IO', hsubr, ymsg )
- kresp_out = iresps(2)
- end if
-#endif
-
-#ifdef MNH_IOCDF4
- if ( iresps(3) /= 0 ) then
- write(yresp, '( i6 )') iresps(3)
- ymsg = Trim( tpfile%cname ) // ': netCDF: resp=' // yresp // ' when writing ' // Trim( tpfield%cmnhname )
- call Print_msg( NVERB_ERROR, 'IO', hsubr, ymsg )
- kresp_out = iresps(3)
- end if
-#endif
-
- end subroutine IO_Field_write_error_check
-
-
- SUBROUTINE IO_Header_write(TPFILE,HDAD_NAME)
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! File structure
- CHARACTER(LEN=*),OPTIONAL,INTENT(IN) :: HDAD_NAME
-
- integer :: ifile
-
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Header_write_FILE','called for file '//TRIM(TPFILE%CNAME))
-
- CALL IO_Header_onefile_write(TPFILE,HDAD_NAME)
-
- !Write header also for the Z-split files
- DO IFILE=1,TPFILE%NSUBFILES_IOZ
- CALL IO_Header_onefile_write(TPFILE%TFILES_IOZ(IFILE)%TFILE,HDAD_NAME)
- END DO
- END SUBROUTINE IO_Header_write
-
-
- SUBROUTINE IO_Header_onefile_write(TPFILE,HDAD_NAME)
- !
- USE MODD_CONF
- USE MODD_CONF_n, ONLY: CSTORAGE_TYPE
- USE MODD_PARAMETERS, ONLY: NFILENAMELGTMAXLFI
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! File structure
- CHARACTER(LEN=*),OPTIONAL,INTENT(IN) :: HDAD_NAME
- !
- CHARACTER(LEN=:),ALLOCATABLE :: YDAD_NAME
- INTEGER :: ILEN,ILEN2
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Header_onefile_write','called for file '//TRIM(TPFILE%CNAME))
- !
- IF ( ASSOCIATED(TPFILE%TDADFILE) .AND. PRESENT(HDAD_NAME) ) THEN
- IF ( TRIM(TPFILE%TDADFILE%CNAME) /= TRIM(HDAD_NAME) ) THEN
- CALL PRINT_MSG(NVERB_ERROR,'IO','IO_Header_onefile_write','TPFILE%TDADFILE%CNAME /= HDAD_NAME')
- END IF
- END IF
- !
- CALL IO_Header_write_nc4(TPFILE)
- !
- CALL IO_Field_write(TPFILE,'MNHVERSION', NMNHVERSION)
- CALL IO_Field_write(TPFILE,'MASDEV', NMASDEV)
- CALL IO_Field_write(TPFILE,'BUGFIX', NBUGFIX)
- CALL IO_Field_write(TPFILE,'BIBUSER', CBIBUSER)
- CALL IO_Field_write(TPFILE,'PROGRAM', CPROGRAM)
- CALL IO_Field_write(TPFILE,'STORAGE_TYPE',CSTORAGE_TYPE)
- CALL IO_Field_write(TPFILE,'MY_NAME', TPFILE%CNAME)
- !
- IF ( ASSOCIATED(TPFILE%TDADFILE) ) THEN
- ILEN = LEN_TRIM(TPFILE%TDADFILE%CNAME)
- ILEN2 = MAX(NFILENAMELGTMAXLFI,ILEN)
- ALLOCATE(CHARACTER(LEN=ILEN2) :: YDAD_NAME)
- IF(ILEN>0) THEN
- YDAD_NAME(1:ILEN) = TPFILE%TDADFILE%CNAME(1:ILEN)
- YDAD_NAME(ILEN+1:ILEN2) = ' '
- ELSE
- YDAD_NAME(:) = ' '
- END IF
- ELSE IF (PRESENT(HDAD_NAME)) THEN
- ILEN = LEN_TRIM(HDAD_NAME)
- ILEN2 = MAX(NFILENAMELGTMAXLFI,ILEN)
- ALLOCATE(CHARACTER(LEN=ILEN2) :: YDAD_NAME)
- IF(ILEN>0) THEN
- YDAD_NAME(1:ILEN) = HDAD_NAME(1:ILEN)
- YDAD_NAME(ILEN+1:ILEN2) = ' '
- ELSE
- YDAD_NAME(:) = ' '
- END IF
- ELSE
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Header_onefile_write',TRIM(TPFILE%CNAME)// &
- ': TPFILE%TDADFILE not associated and HDAD_NAME not provided')
- ALLOCATE(CHARACTER(LEN=NFILENAMELGTMAXLFI) :: YDAD_NAME)
- YDAD_NAME(:) = ' '
- ENDIF
- CALL IO_Field_write(TPFILE,'DAD_NAME',YDAD_NAME)
- DEALLOCATE(YDAD_NAME)
- !
- END SUBROUTINE IO_Header_onefile_write
-
-
-subroutine IO_Field_create( tpfile, tpfield )
- ! Subroutine to create a variable in a file and write its metadata without writing its data
- ! LFI files are not supported
- use modd_field
- use modd_io, only: gsmonoproc, isp
-
- type(tfiledata), intent(in) :: tpfile
- class(tfieldmetadata), intent(in) :: tpfield
-
- integer :: ik_file
- integer :: iresp
- logical :: glfi, gnc4
- class(tfieldmetadata), allocatable :: tzfield
- type(tfiledata), pointer :: tzfile
-
- call Print_msg( NVERB_DEBUG, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': creating ' // Trim( tpfield%cmnhname ) )
-
- if ( Any (tpfield%ndimlist(:) == NMNHDIM_UNKNOWN ) ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': ' // Trim( tpfield%cmnhname ) &
- // ' : ndimlist must be populated' )
- return
- end if
-
- !Not very useful: call IO_Field_metadata_check( tpfield, tpfield%ntype, tpfield%ndims, 'IO_Field_create' )
-
- call IO_File_write_check( tpfile, 'IO_Field_create', iresp )
-
- call IO_Format_write_select( tpfile, glfi, gnc4 )
-
- if ( glfi ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': LFI format not supported' )
- glfi = .false.
- end if
-
- if ( iresp == 0 ) then
- Allocate( tzfield, source = tpfield )
-
- if ( All( tzfield%ntype /= [ TYPEINT, TYPELOG, TYPEREAL, TYPECHAR, TYPEDATE ] ) ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': ' &
- // Trim( tzfield%cmnhname ) // ': invalid ntype' )
- return
- end if
-
- NDIMS: select case ( tzfield%ndims )
- case ( 0 ) NDIMS
-#if MNH_SCALARS_IN_SPLITFILES
- if ( Any( tzfield%ntype == [ TYPEREAL, TYPEINT, TYPELOG ] ) then
- if ( tpfile%nsubfiles_ioz > 0 ) then
- !Create the variable in all the Z-split files
- do ik_file = 1, tpfile%nsubfiles_ioz
- tzfile => tpfile%tfiles_ioz(ik_file)%tfile
- if ( isp == tzfile%nmaster_rank ) then
-#ifdef MNH_IOCDF4
- if ( gnc4 ) call IO_Field_create_nc4( tzfile, tzfield )
-#endif
- end if
- end do
- endif
- end if
-#endif
- case ( 1 ) NDIMS
- ! Nothing to do
-
- case ( 2 ) NDIMS
- if ( All( tzfield%ntype /= [ TYPEINT, TYPEREAL ] ) ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': ' &
- // Trim( tzfield%cmnhname ) // ': invalid ntype for 2D field' )
- return
- end if
-
- if ( gsmonoproc ) call IO_Ndimlist_reduce( tpfile, tzfield )
-
- case ( 3 ) NDIMS
- if ( All( tzfield%ntype /= [ TYPEINT, TYPEREAL ] ) ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': ' &
- // Trim( tzfield%cmnhname ) // ': invalid ntype for 2D field' )
- return
- end if
-
- if ( gsmonoproc .and. ( tzfield%ntype /= TYPEREAL .or. tpfile%nsubfiles_ioz == 0 ) ) &
- call IO_Ndimlist_reduce( tpfile, tzfield )
-
- if ( tzfield%ntype == TYPEREAL .and. tpfile%nsubfiles_ioz > 0 ) then
-#ifdef MNH_IOCDF4
- ! Write the variable attributes in the non-split file
- if ( tpfile%nmaster_rank==isp .and. gnc4 ) &
- call IO_Field_header_split_write_nc4( tpfile, tzfield, &
- Int( tpfile%tncdims%tdims(tzfield%ndimlist(3))%nlen, kind = Kind( 1 ) ) )
- end if
-#endif
-
- case ( 4 ) NDIMS
- if ( All( tzfield%ntype /= [ TYPEINT, TYPEREAL ] ) ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': ' &
- // Trim( tzfield%cmnhname ) // ': invalid ntype for 4D field' )
- return
- end if
-
- if ( gsmonoproc ) call IO_Ndimlist_reduce( tpfile, tzfield )
-
- case ( 5 ) NDIMS
- if ( tzfield%ntype /= TYPEREAL ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': ' &
- // Trim( tzfield%cmnhname ) // ': invalid ntype for 5D field' )
- return
- end if
-
- if ( gsmonoproc ) call IO_Ndimlist_reduce( tpfile, tzfield )
-
- case ( 6 ) NDIMS
- if ( tzfield%ntype /= TYPEREAL ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': ' &
- // Trim( tzfield%cmnhname ) // ': invalid ntype for 6D field' )
- return
- end if
-
- ! Nothing else to do
-
- case default NDIMS
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_create', Trim( tpfile%cname ) // ': ' &
- // Trim( tzfield%cmnhname ) // ': invalid ndims' )
- end select NDIMS
-
- if ( isp == tpfile%nmaster_rank ) then
-#ifdef MNH_IOCDF4
- if ( gnc4 ) call IO_Field_create_nc4( tpfile, tzfield )
-#endif
- end if
- end if
-
-end subroutine IO_Field_create
-
-
-subroutine IO_Ndimlist_reduce( tpfile, tpfield )
- use modd_field, only: NMNHDIM_ONE, NMNHDIM_UNKNOWN, NMNHDIM_UNUSED
- use modd_io, only: gsmonoproc, l1d, l2d, lpack
- use modd_parameters_ll, only: jphext
-
- type(tfiledata), intent(in) :: tpfile
- class(tfieldmetadata_base), intent(inout) :: tpfield
-
- integer :: ihextot
- integer :: ji
-
- if ( .not. gsmonoproc ) return
-
- ihextot = 2*jphext+1
-
- ! sequential execution and non Z-split field
- if ( tpfield%ndims /= 3 .or. tpfield%ntype /= TYPEREAL .or. tpfile%nsubfiles_ioz == 0 ) then
- if ( lpack .and. l1d .and. tpfile%tncdims%tdims(tpfield%ndimlist(1))%nlen == ihextot &
- .and. tpfile%tncdims%tdims(tpfield%ndimlist(2))%nlen == ihextot ) then
- if ( tpfile%ldimreduced ) then
- tpfield%ndims = tpfield%ndims - 2
- if ( tpfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- ! Last iteration necessary if time dimension
- do ji = 1, tpfield%ndims + 1
- tpfield%ndimlist(ji) = tpfield%ndimlist(ji + 2)
- end do
- tpfield%ndimlist(tpfield%ndims + 2 : ) = NMNHDIM_UNUSED
- end if
- else
- if ( tpfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tpfield%ndimlist(1:2) = NMNHDIM_ONE
- end if
- endif
- else if ( lpack .and. l2d .and. tpfile%tncdims%tdims(tpfield%ndimlist(2))%nlen == ihextot ) then
- if ( tpfile%ldimreduced ) then
- tpfield%ndims = tpfield%ndims - 1
- if ( tpfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- ! Last iteration necessary if time dimension
- do ji = 2, tpfield%ndims + 1
- tpfield%ndimlist(ji) = tpfield%ndimlist(ji + 1)
- end do
- tpfield%ndimlist(tpfield%ndims + 2 : ) = NMNHDIM_UNUSED
- end if
- else
- if ( tpfield%ndimlist(2) /= NMNHDIM_UNKNOWN ) then
- tpfield%ndimlist(2) = NMNHDIM_ONE
- end if
- endif
- else
- !Nothing to do
- end if
- end if
-
-end subroutine IO_Ndimlist_reduce
-
-
- SUBROUTINE IO_Field_write_byname_X0(TPFILE,HNAME,PFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- REAL, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_X0',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),PFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_X0
-
-
- SUBROUTINE IO_Field_write_byfield_X0(TPFILE,TPFIELD,PFIELD,KRESP)
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- USE MODE_IO_MANAGE_STRUCT, ONLY: IO_File_find_byname
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- !
- INTEGER :: IK_FILE
- TYPE(TFILEDATA),POINTER :: TZFILE
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- TZFILE => NULL()
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_X0',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEREAL,0,'IO_Field_write_byfield_X0')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_X0',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,PFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,PFIELD,iresp_nc4)
- ELSE ! multiprocesses execution
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,PFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,PFIELD,iresp_nc4)
- END IF
- !
- END IF ! multiprocesses execution
-#if MNH_SCALARS_IN_SPLITFILES
- IF (TPFILE%NSUBFILES_IOZ>0) THEN
- ! write the data in all Z files
- DO IK_FILE=1,TPFILE%NSUBFILES_IOZ
- TZFILE => TPFILE%TFILES_IOZ(IK_FILE)%TFILE
- IF ( ISP == TZFILE%NMASTER_RANK ) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TZFILE,TPFIELD,PFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TZFILE,TPFIELD,PFIELD,iresp_nc4)
- END IF
- END DO
- ENDIF
-#endif
- END IF
- !
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_X0', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
- END SUBROUTINE IO_Field_write_byfield_X0
-
-
- SUBROUTINE IO_Field_write_byname_X1( TPFILE, HNAME, PFIELD, KRESP, koffset )
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- REAL,DIMENSION(:), INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(1), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return-code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_X1',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- if( iresp == 0 ) call IO_Field_write( tpfile, tfieldlist(id), pfield, iresp, koffset )
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_X1
-
- SUBROUTINE IO_Field_write_phy_byfield_X2(D, TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(DIMPHYEX_t), INTENT(IN) :: D
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- TYPE(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(D%NIJT,D%NKT),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(3), optional, intent(in) :: koffset
- !
- CALL IO_Field_write_phy_unpack2D(D,TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- !
- END SUBROUTINE IO_Field_write_phy_byfield_X2
-!
- SUBROUTINE IO_Field_write_phy_unpack2D(D, TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(DIMPHYEX_t), INTENT(IN) :: D
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- TYPE(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(D%NIT,D%NJT,D%NKT),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(3), optional, intent(in) :: koffset
- !
- CALL IO_Field_write_byfield_X3(TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- !
- END SUBROUTINE IO_Field_write_phy_unpack2D
-!
- SUBROUTINE IO_Field_write_phy_byfield_X1(D, TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(DIMPHYEX_t), INTENT(IN) :: D
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- TYPE(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(D%NIJT),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(3), optional, intent(in) :: koffset
- !
- CALL IO_Field_write_phy_unpack1D(D,TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- !
- END SUBROUTINE IO_Field_write_phy_byfield_X1
-!
- SUBROUTINE IO_Field_write_phy_unpack1D(D, TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(DIMPHYEX_t), INTENT(IN) :: D
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- TYPE(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(D%NIT,D%NJT),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(3), optional, intent(in) :: koffset
- !
- CALL IO_Field_write_byfield_X2(TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- !
- END SUBROUTINE IO_Field_write_phy_unpack1D
-
- SUBROUTINE IO_Field_write_byfield_X1( TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- USE MODE_IO_MANAGE_STRUCT, ONLY: IO_File_find_byname
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(:),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(1), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- INTEGER :: ISIZEMAX
- REAL,DIMENSION(:),POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- !
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_X1',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEREAL,1,'IO_Field_write_byfield_X1')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_X1',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, pfield, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, pfield, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, pfield, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, pfield, iresp_nc4 )
- end if
- ELSE ! multiprocesses execution
- CALL MPI_ALLREDUCE(SIZE(PFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_X1','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- CALL ALLOCBUFFER_ll(ZFIELDP,PFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(ZFIELDP(0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- END IF
- !
- CALL MPI_BCAST(IRESP,1,MNHINT_MPI,TPFILE%NMASTER_RANK-1,TPFILE%NMPICOMM,IERR)
- END IF ! multiprocesses execution
- END IF
- !
- IF (GALLOC) DEALLOCATE(ZFIELDP)
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_X1', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
- END SUBROUTINE IO_Field_write_byfield_X1
-
-
- SUBROUTINE IO_Field_write_byname_X2( TPFILE, HNAME, PFIELD, KRESP, koffset )
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- REAL,DIMENSION(:,:), INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(2), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return-code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_X2',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- if( iresp == 0 ) call IO_Field_write( tpfile, tfieldlist(id), pfield, iresp, koffset )
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_X2
-
-
- SUBROUTINE IO_Field_write_byfield_X2( TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- use modd_field, only: NMNHDIM_UNKNOWN, NMNHDIM_ONE, NMNHDIM_UNUSED
- USE MODD_IO, ONLY: GSMONOPROC,ISP,L1D,L2D,LPACK
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_TIMEZ, ONLY: TIMEZ
- !
- USE MODE_ALLOCBUFFER_ll
-#ifdef MNH_GA
- USE MODE_GA
-#endif
- USE MODE_GATHER_ll
- USE MODE_MNH_TIMING, ONLY: SECOND_MNH2
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(:,:),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(2), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- integer, dimension(1) :: ioffset1d
- real :: zfieldp0d
- real, dimension(:), pointer :: zfieldp1d
- REAL, DIMENSION(:,:), POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- !
- REAL(kind=MNHTIME), DIMENSION(2) :: ZT0, ZT1, ZT2
- REAL(kind=MNHTIME), DIMENSION(2) :: ZT11, ZT22
-#ifdef MNH_GA
- REAL,DIMENSION(:,:),POINTER :: ZFIELD_GA
-#endif
- INTEGER :: IHEXTOT
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- class(tfieldmetadata), allocatable :: tzfield
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- IHEXTOT = 2*JPHEXT+1
- !
- CALL SECOND_MNH2(ZT11)
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_X2',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEREAL,2,'IO_Field_write_byfield_X2')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_X2',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- if ( Present( koffset ) .and. glfi ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_partial_write_byfield_X2', Trim( tpfile%cname ) // ': LFI format not supported' )
- glfi = .false.
- end if
-
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- ! IF (LPACK .AND. L1D .AND. YDIR=='XY') THEN
- IF (LPACK .AND. L1D .AND. SIZE(PFIELD,1)==IHEXTOT .AND. SIZE(PFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 2
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1) = tzfield%ndimlist(3) !Necessary if time dimension
- tzfield%ndimlist(2:) = NMNHDIM_UNUSED
- end if
- zfieldp0d = pfield(jphext + 1, jphext + 1)
- if ( Present( koffset ) ) then
- call Print_msg( NVERB_FATAL, 'IO', 'IO_Field_partial_write_byfield_X2', Trim( tpfile%cname ) &
- // ': impossible situation/not implemented' )
- !!if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp0d, ioffset0d, iresp_lfi )
- !if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp0d, ioffset0d, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp0d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp0d, iresp_nc4 )
- end if
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1:2) = NMNHDIM_ONE
- end if
- zfieldp => pfield(jphext + 1 : jphext + 1, jphext + 1 : jphext + 1)
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp, iresp_nc4 )
- end if
- endif
- ! ELSE IF (LPACK .AND. L2D .AND. YDIR=='XY') THEN
- ELSEIF (LPACK .AND. L2D .AND. SIZE(PFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 1
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = tzfield%ndimlist(3) !Necessary if time dimension
- tzfield%ndimlist(3:) = NMNHDIM_UNUSED
- end if
- zfieldp1d => pfield(:, jphext + 1)
- if ( Present( koffset ) ) then
- ioffset1d(1) = koffset(1)
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp1d, ioffset1d, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp1d, ioffset1d, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp1d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp1d, iresp_nc4 )
- end if
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = NMNHDIM_ONE
- end if
- zfieldp => pfield(:, jphext + 1 : jphext + 1)
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp, iresp_nc4 )
- end if
- endif
- ELSE
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, pfield, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, pfield, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, pfield, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, pfield, iresp_nc4 )
- end if
- END IF
- ELSE ! multiprocesses execution
- CALL SECOND_MNH2(ZT0)
- CALL MPI_ALLREDUCE(SIZE(PFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_X2','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- ! I/O process case
- CALL ALLOCBUFFER_ll(ZFIELDP,PFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(ZFIELDP(0,0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'IO', 'IO_Field_write_byfield_X2', '2D not (yet) allowed for parallel execution' )
- CALL GATHER_XXFIELD('XX',PFIELD(:,JPHEXT+1),ZFIELDP(:,1),TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSE
-#ifdef MNH_GA
- !
- ! init/create the ga , dim3 = 1
- !
- CALL MNH_INIT_GA(SIZE(PFIELD,1),SIZE(PFIELD,2),1,YRECFM,"WRITE")
- !
- ! copy columun data to global arrays g_a
- !
- ALLOCATE (ZFIELD_GA (SIZE(PFIELD,1),SIZE(PFIELD,2)))
- ZFIELD_GA = PFIELD
- !print*,"IO_WRITE_FIELD_BYFIELD_X2::nga_put=",g_a, lo_col, hi_col,NIXO_L,NIYO_L , ld_col, YRECFM ; call flush(6)
- call ga_sync()
- call nga_put(g_a, lo_col, hi_col,ZFIELD_GA(NIXO_L,NIYO_L) , ld_col)
- call ga_sync()
- DEALLOCATE (ZFIELD_GA)
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- !
- ! this proc get the Z slide to write
- !
- lo_zplan(JPIZ) = 1
- hi_zplan(JPIZ) = 1
- !print*,"IO_WRITE_FIELD_BYFIELD_X2::nga_get=",g_a, lo_zplan, hi_zplan, ld_zplan, YRECFM ; call flush(6)
- call nga_get(g_a, lo_zplan, hi_zplan,ZFIELDP, ld_zplan)
- END IF
-#else
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
-#endif
- END IF
- END IF
- CALL SECOND_MNH2(ZT1)
- TIMEZ%T_WRIT2D_GATH=TIMEZ%T_WRIT2D_GATH + ZT1 - ZT0
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- END IF
-#ifdef MNH_GA
- call ga_sync
-#endif
- CALL SECOND_MNH2(ZT2)
- TIMEZ%T_WRIT2D_WRIT=TIMEZ%T_WRIT2D_WRIT + ZT2 - ZT1
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_X2', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- CALL SECOND_MNH2(ZT22)
- TIMEZ%T_WRIT2D_ALL=TIMEZ%T_WRIT2D_ALL + ZT22 - ZT11
- END SUBROUTINE IO_Field_write_byfield_X2
-
-
- SUBROUTINE IO_Field_write_byname_X3( TPFILE, HNAME, PFIELD, KRESP, koffset )
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- REAL, DIMENSION(:,:,:), INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(3), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_X3',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- if( iresp == 0 ) call IO_Field_write( tpfile, tfieldlist(id), pfield, iresp, koffset )
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_X3
-
-
- SUBROUTINE IO_Field_write_byfield_X3( TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- use modd_field, only: NMNHDIM_UNKNOWN, NMNHDIM_ONE, NMNHDIM_UNUSED
- USE MODD_IO, ONLY: GSMONOPROC, ISNPROC, ISP, L1D, L2D, LPACK
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_TIMEZ, ONLY: TIMEZ
- USE MODD_VAR_ll, ONLY: MNH_STATUSES_IGNORE
- !
- USE MODE_ALLOCBUFFER_ll
-#ifdef MNH_GA
- USE MODE_GA
-#endif
- USE MODE_GATHER_ll
- USE MODE_IO_MANAGE_STRUCT, ONLY: IO_File_find_byname
- USE MODE_IO_TOOLS, ONLY: IO_Level2filenumber_get
- USE MODE_MNH_TIMING, ONLY: SECOND_MNH2
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA),TARGET, INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(:,:,:),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(3), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- integer :: iresp_tmp_lfi, iresp_tmp_nc4
- integer, dimension(2) :: iresps
- integer, dimension(1) :: ioffset1d
- integer, dimension(2) :: ioffset2d
- real,dimension(:), pointer :: zfieldp1d
- real,dimension(:,:), pointer :: zfieldp2d
- REAL,DIMENSION(:,:,:),POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- INTEGER :: JK,JKK
- REAL,DIMENSION(:,:),POINTER :: ZSLICE_ll,ZSLICE
- INTEGER :: IK_FILE,IK_RANK,INB_PROC_REAL,JK_MAX
- INTEGER :: JI,IXO,IXE,IYO,IYE
- REAL,DIMENSION(:,:),POINTER :: ZTX2DP
- INTEGER, DIMENSION(MPI_STATUS_SIZE) :: ISTATUS
- LOGICAL :: GALLOC_ll
- INTEGER,ALLOCATABLE,DIMENSION(:) :: IREQ_TAB
- INTEGER :: INB_REQ
- TYPE TX_2DP
- REAL, DIMENSION(:,:), POINTER :: X
- END TYPE TX_2DP
- TYPE(TX_2DP),ALLOCATABLE,DIMENSION(:) :: T_TX2DP
- REAL(kind=MNHTIME), DIMENSION(2) :: ZT0, ZT1, ZT2
- REAL(kind=MNHTIME), DIMENSION(2) :: ZT11, ZT22
-#ifdef MNH_GA
- REAL,DIMENSION(:,:,:),POINTER :: ZFIELD_GA
-#endif
- INTEGER :: IHEXTOT
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- class(tfieldmetadata), allocatable :: tzfield
- TYPE(TFILEDATA),POINTER :: TZFILE
- !
- TZFILE => NULL()
- !
- ZSLICE => NULL()
- ZSLICE_ll => NULL()
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- GALLOC_ll = .FALSE.
- IHEXTOT = 2*JPHEXT+1
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_X3',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL SECOND_MNH2(ZT11)
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEREAL,3,'IO_Field_write_byfield_X3')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_X3',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- if ( Present( koffset ) .and. glfi ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_partial_write_byfield_X3', Trim( tpfile%cname ) // ': LFI format not supported' )
- glfi = .false.
- end if
-
- IF (IRESP==0) THEN
- IF (GSMONOPROC .AND. TPFILE%NSUBFILES_IOZ==0 ) THEN ! sequential execution
- ! IF (LPACK .AND. L1D .AND. YDIR=='XY') THEN
- IF (LPACK .AND. L1D .AND. SIZE(PFIELD,1)==IHEXTOT .AND. SIZE(PFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 2
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1) = tzfield%ndimlist(3)
- tzfield%ndimlist(2) = tzfield%ndimlist(4) !Necessary if time dimension
- tzfield%ndimlist(3:) = NMNHDIM_UNUSED
- end if
- zfieldp1d => pfield(jphext + 1, jphext + 1, :)
- if ( Present ( koffset ) ) then
- ioffset1d(1) = koffset(3)
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp1d, ioffset1d, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp1d, ioffset1d, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp1d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp1d, iresp_nc4 )
- end if
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1:2) = NMNHDIM_ONE
- end if
- zfieldp => pfield(jphext + 1 : jphext + 1, jphext + 1 : jphext + 1, :)
- if ( Present ( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp, iresp_nc4 )
- end if
- endif
- ELSEIF (LPACK .AND. L2D .AND. SIZE(PFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 1
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = tzfield%ndimlist(3)
- tzfield%ndimlist(3) = tzfield%ndimlist(4) !Necessary if time dimension
- tzfield%ndimlist(4:) = NMNHDIM_UNUSED
- end if
- zfieldp2d => pfield(:, jphext + 1, :)
- if ( Present ( koffset ) ) then
- ioffset2d(1) = koffset(1)
- ioffset2d(2) = koffset(3)
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp2d, ioffset2d, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp2d, ioffset2d, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp2d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp2d, iresp_nc4 )
- end if
- else
- if ( tzfield%ndimlist(2) /= NMNHDIM_UNKNOWN ) tzfield%ndimlist(2) = NMNHDIM_ONE
- zfieldp => pfield(:, jphext + 1 : jphext + 1, :)
- if ( Present ( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp, iresp_nc4 )
- end if
- endif
- ELSE
- if ( Present ( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, pfield, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, pfield, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, pfield, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, pfield, iresp_nc4 )
- end if
- END IF
- ELSEIF ( TPFILE%NSUBFILES_IOZ==0 .OR. YDIR=='--' ) THEN ! multiprocesses execution & 1 proc IO
- CALL MPI_ALLREDUCE(SIZE(PFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_X3','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- ! write 3D field in 1 time = output for graphique
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- CALL ALLOCBUFFER_ll(ZFIELDP,PFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(ZFIELDP(0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'IO', 'IO_Field_write_byfield_X3', '2D not (yet) allowed for parallel execution' )
- CALL GATHER_XXFIELD('XX',PFIELD(:,JPHEXT+1,:),ZFIELDP(:,1,:),TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSE
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- if ( Present ( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if (glfi) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if (gnc4) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- END IF
- !
- ELSE ! multiprocesses execution & // IO
- CALL MPI_ALLREDUCE(SIZE(PFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_X3','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
- !
- ! Write the variable attributes in the non-split file
- !
- if ( tpfile%nmaster_rank==isp .and. gnc4 ) &
- call IO_Field_header_split_write_nc4( tpfile, tpfield, size( pfield, 3 ) )
- !
- !JUAN BG Z SLICE
- !
- !
-#ifdef MNH_GA
- !
- ! init/create the ga
- !
- CALL SECOND_MNH2(ZT0)
- CALL MNH_INIT_GA(SIZE(PFIELD,1),SIZE(PFIELD,2),SIZE(PFIELD,3),YRECFM,"WRITE")
- !
- ! copy columun data to global arrays g_a
- !
- ALLOCATE (ZFIELD_GA (SIZE(PFIELD,1),SIZE(PFIELD,2),SIZE(PFIELD,3)))
- ZFIELD_GA = PFIELD
- !print*,"IO_WRITE_FIELD_BYFIELD_X3::nga_put=",g_a, lo_col, hi_col,NIXO_L,NIYO_L , ld_col, YRECFM ; call flush(6)
- call ga_sync()
- call nga_put(g_a, lo_col, hi_col,ZFIELD_GA(NIXO_L,NIYO_L,1) , ld_col)
- call ga_sync()
- DEALLOCATE(ZFIELD_GA)
- CALL SECOND_MNH2(ZT1)
- TIMEZ%T_WRIT3D_SEND=TIMEZ%T_WRIT3D_SEND + ZT1 - ZT0
- !
- ! write the data
- !
- ALLOCATE(ZSLICE_ll(0,0)) ! to avoid bug on test of size
- GALLOC_ll = .TRUE.
- !
- DO JKK=1,SIZE(PFIELD,3) ! IKU_ll
- !
- IK_FILE = IO_Level2filenumber_get(JKK,TPFILE%NSUBFILES_IOZ)
- TZFILE => TPFILE%TFILES_IOZ(IK_FILE+1)%TFILE
- !
- IK_RANK = TZFILE%NMASTER_RANK
- !
- IF (ISP == IK_RANK ) THEN
- CALL SECOND_MNH2(ZT0)
- !
- IF ( SIZE(ZSLICE_ll) .EQ. 0 ) THEN
- DEALLOCATE(ZSLICE_ll)
- CALL ALLOCBUFFER_ll(ZSLICE_ll,ZSLICE,YDIR,GALLOC_ll)
- END IF
- !
- ! this proc get this JKK slide
- !
- lo_zplan(JPIZ) = JKK
- hi_zplan(JPIZ) = JKK
- !print*,"IO_WRITE_FIELD_BYFIELD_X3::nga_get=",g_a, lo_zplan, hi_zplan, ld_zplan, YRECFM,JKK ; call flush(6)
- call nga_get(g_a, lo_zplan, hi_zplan,ZSLICE_ll, ld_zplan)
- CALL SECOND_MNH2(ZT1)
- TIMEZ%T_WRIT3D_RECV=TIMEZ%T_WRIT3D_RECV + ZT1 - ZT0
- !
- if ( Present ( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zslice_ll, koffset, iresp_tmp_lfi, &
- ! kvertlevel = jkk, kzfile = ik_file + 1 )
- !if ( iresp_tmp_lfi /= 0 ) iresp_lfi = iresp_tmp_lfi
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zslice_ll, koffset, iresp_tmp_nc4, &
- kvertlevel = jkk, kzfile = ik_file + 1 )
- if ( iresp_tmp_nc4 /= 0 ) iresp_nc4 = iresp_tmp_nc4
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zslice_ll, iresp_tmp_lfi, &
- kvertlevel = jkk, kzfile = ik_file + 1 )
- if ( iresp_tmp_lfi /= 0 ) iresp_lfi = iresp_tmp_lfi
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zslice_ll, iresp_tmp_nc4, &
- kvertlevel = jkk, kzfile = ik_file + 1 )
- if ( iresp_tmp_nc4 /= 0 ) iresp_nc4 = iresp_tmp_nc4
- end if
- CALL SECOND_MNH2(ZT2)
- TIMEZ%T_WRIT3D_WRIT=TIMEZ%T_WRIT3D_WRIT + ZT2 - ZT1
- END IF
- END DO
- !
- CALL SECOND_MNH2(ZT0)
- call ga_sync
- CALL SECOND_MNH2(ZT1)
- TIMEZ%T_WRIT3D_WAIT=TIMEZ%T_WRIT3D_WAIT + ZT1 - ZT0
-#else
- !
- ALLOCATE(ZSLICE_ll(0,0))
- GALLOC_ll = .TRUE.
- INB_PROC_REAL = MIN(TPFILE%NSUBFILES_IOZ,ISNPROC)
- Z_SLICE: DO JK=1,SIZE(PFIELD,3),INB_PROC_REAL
- !
- ! collect the data
- !
- JK_MAX=MIN(SIZE(PFIELD,3),JK+INB_PROC_REAL-1)
- !
- INB_REQ=0
- ALLOCATE(IREQ_TAB(INB_PROC_REAL))
- ALLOCATE(T_TX2DP(INB_PROC_REAL))
- DO JKK=JK,JK_MAX
- !
- ! get the file & rank to write this level
- !
- IF (TPFILE%NSUBFILES_IOZ .GT. 1 ) THEN
- IK_FILE = IO_Level2filenumber_get(JKK,TPFILE%NSUBFILES_IOZ)
- TZFILE => TPFILE%TFILES_IOZ(IK_FILE+1)%TFILE
- ELSE
- TZFILE => TPFILE
- END IF
- !
- IK_RANK = TZFILE%NMASTER_RANK
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- call Print_msg( NVERB_FATAL, 'GEN', 'IO_Field_write_byfield_X3', 'XX/YY not (yet) allowed for parallel I/O' )
- CALL GATHER_XXFIELD(YDIR,PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'GEN', 'IO_Field_write_byfield_X3', &
- '2D not (yet) allowed for parallel execution' )
- CALL GATHER_XXFIELD('XX',PFIELD(:,JPHEXT+1,:),ZFIELDP(:,1,:),TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSE
- CALL SECOND_MNH2(ZT0)
- IF ( ISP /= IK_RANK ) THEN
- ! Other processes
- CALL GET_DOMWRITE_ll(ISP,'local',IXO,IXE,IYO,IYE)
- IF (IXO /= 0) THEN ! intersection is not empty
- INB_REQ = INB_REQ + 1
- ALLOCATE(T_TX2DP(INB_REQ)%X(IXO:IXE,IYO:IYE))
- ZSLICE => PFIELD(:,:,JKK)
- ZTX2DP=>ZSLICE(IXO:IXE,IYO:IYE)
- T_TX2DP(INB_REQ)%X=ZSLICE(IXO:IXE,IYO:IYE)
- CALL MPI_ISEND(T_TX2DP(INB_REQ)%X,SIZE(ZTX2DP),MNHREAL_MPI,IK_RANK-1,99+IK_RANK &
- & ,TZFILE%NMPICOMM,IREQ_TAB(INB_REQ),IERR)
- !CALL MPI_BSEND(ZTX2DP,SIZE(ZTX2DP),MNHREAL_MPI,IK_RANK-1,99+IK_RANK,TZFILE%NMPICOMM,IERR)
- END IF
- END IF
- CALL SECOND_MNH2(ZT1)
- TIMEZ%T_WRIT3D_SEND=TIMEZ%T_WRIT3D_SEND + ZT1 - ZT0
- END IF
- END IF
- END DO
- !
- ! Write the variable attributes in the non-split file
- !
- if ( tpfile%nmaster_rank == isp .and. gnc4 ) &
- call IO_Field_header_split_write_nc4( tpfile, tpfield, size( pfield, 3 ) )
- !
- ! write the data
- !
- DO JKK=JK,JK_MAX
- IF (TPFILE%NSUBFILES_IOZ .GT. 1 ) THEN
- IK_FILE = IO_Level2filenumber_get(JKK,TPFILE%NSUBFILES_IOZ)
- TZFILE => TPFILE%TFILES_IOZ(IK_FILE+1)%TFILE
- ELSE
- TZFILE => TPFILE
- ENDIF
- IK_RANK = TZFILE%NMASTER_RANK
- !
- IF (ISP == IK_RANK ) THEN
- CALL SECOND_MNH2(ZT0)
- ! I/O proc case
- IF ( SIZE(ZSLICE_ll) .EQ. 0 ) THEN
- DEALLOCATE(ZSLICE_ll)
- CALL ALLOCBUFFER_ll(ZSLICE_ll,ZSLICE,YDIR,GALLOC_ll)
- END IF
- DO JI=1,ISNPROC
- CALL GET_DOMWRITE_ll(JI,'global',IXO,IXE,IYO,IYE)
- IF (IXO /= 0) THEN ! intersection is not empty
- ZTX2DP=>ZSLICE_ll(IXO:IXE,IYO:IYE)
- IF (ISP == JI) THEN
- CALL GET_DOMWRITE_ll(JI,'local',IXO,IXE,IYO,IYE)
- ZSLICE => PFIELD(:,:,JKK)
- ZTX2DP = ZSLICE(IXO:IXE,IYO:IYE)
- ELSE
- CALL MPI_RECV(ZTX2DP,SIZE(ZTX2DP),MNHREAL_MPI,JI-1,99+IK_RANK,TZFILE%NMPICOMM,ISTATUS,IERR)
- END IF
- END IF
- END DO
- CALL SECOND_MNH2(ZT1)
- TIMEZ%T_WRIT3D_RECV=TIMEZ%T_WRIT3D_RECV + ZT1 - ZT0
- if ( Present ( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zslice_ll, koffset, iresp_tmp_lfi, &
- ! kvertlevel = jkk, kzfile = ik_file + 1 )
- !if ( iresp_tmp_lfi /= 0 ) iresp_lfi = iresp_tmp_lfi
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zslice_ll, koffset, iresp_tmp_nc4, &
- kvertlevel = jkk, kzfile = ik_file + 1 )
- if ( iresp_tmp_nc4 /= 0 ) iresp_nc4 = iresp_tmp_nc4
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zslice_ll, iresp_tmp_lfi, &
- kvertlevel = jkk, kzfile = ik_file + 1 )
- if ( iresp_tmp_lfi /= 0 ) iresp_lfi = iresp_tmp_lfi
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zslice_ll, iresp_tmp_nc4, &
- kvertlevel = jkk, kzfile = ik_file + 1 )
- if ( iresp_tmp_nc4 /= 0 ) iresp_nc4 = iresp_tmp_nc4
- end if
- CALL SECOND_MNH2(ZT2)
- TIMEZ%T_WRIT3D_WRIT=TIMEZ%T_WRIT3D_WRIT + ZT2 - ZT1
- END IF
- END DO
- !
- CALL SECOND_MNH2(ZT0)
- IF (INB_REQ .GT.0 ) THEN
- CALL MPI_WAITALL(INB_REQ,IREQ_TAB,MNH_STATUSES_IGNORE,IERR)
- DO JI=1,INB_REQ ; DEALLOCATE(T_TX2DP(JI)%X) ; ENDDO
- END IF
- DEALLOCATE(T_TX2DP)
- DEALLOCATE(IREQ_TAB)
- CALL SECOND_MNH2(ZT1)
- TIMEZ%T_WRIT3D_WAIT=TIMEZ%T_WRIT3D_WAIT + ZT1 - ZT0
- END DO Z_SLICE
- !JUAN BG Z SLICE
-! end of MNH_GA
-#endif
- !Not global reduction because a broadcast is done in IO_Field_write_error_check
- call MPI_REDUCE( -Abs( [ iresp_lfi, iresp_nc4 ] ), iresps(:), 2, MNHINT_MPI, MPI_MIN, &
- tpfile%nmaster_rank - 1, tpfile%nmpicomm, ierr )
- iresp_lfi = iresps(1)
- iresp_nc4 = iresps(2)
- END IF ! multiprocesses execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_X3', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- IF (GALLOC_ll) DEALLOCATE(ZSLICE_ll)
- CALL SECOND_MNH2(ZT22)
- TIMEZ%T_WRIT3D_ALL=TIMEZ%T_WRIT3D_ALL + ZT22 - ZT11
- END SUBROUTINE IO_Field_write_byfield_X3
-
-
- SUBROUTINE IO_Field_write_byname_X4( TPFILE, HNAME, PFIELD, KRESP, koffset )
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- REAL,DIMENSION(:,:,:,:), INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(4), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_X4',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- if( iresp == 0 ) call IO_Field_write( tpfile, tfieldlist(id), pfield, iresp, koffset )
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_X4
-
-
- SUBROUTINE IO_Field_write_byfield_X4( TPFILE, TPFIELD, PFIELD, KRESP, koffset )
- use modd_field, only: NMNHDIM_UNKNOWN, NMNHDIM_ONE, NMNHDIM_UNUSED
- USE MODD_IO, ONLY: GSMONOPROC, ISP, L1D, L2D, LPACK
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_TIMEZ, ONLY: TIMEZ
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- USE MODE_IO_TOOLS, ONLY: IO_Level2filenumber_get
- USE MODE_MNH_TIMING, ONLY: SECOND_MNH2
- USE MODD_VAR_ll, ONLY: MNH_STATUSES_IGNORE
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(:,:,:,:),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(4), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- integer, dimension(2) :: ioffset2d
- integer, dimension(3) :: ioffset3d
- real,dimension(:,:), pointer :: zfieldp2d
- real,dimension(:,:,:), pointer :: zfieldp3d
- REAL,DIMENSION(:,:,:,:),POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- INTEGER :: IHEXTOT
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- class(tfieldmetadata), allocatable :: tzfield
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- IHEXTOT = 2*JPHEXT+1
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_X4',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEREAL,4,'IO_Field_write_byfield_X4')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_X4',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- if ( Present( koffset ) .and. glfi ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_partial_write_byfield_X4', Trim( tpfile%cname ) // ': LFI format not supported' )
- glfi = .false.
- end if
-
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- ! IF (LPACK .AND. L1D .AND. YDIR=='XY') THEN
- IF (LPACK .AND. L1D .AND. SIZE(PFIELD,1)==IHEXTOT .AND. SIZE(PFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 2
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1) = tzfield%ndimlist(3)
- tzfield%ndimlist(2) = tzfield%ndimlist(4)
- tzfield%ndimlist(3) = tzfield%ndimlist(5) !Necessary if time dimension
- tzfield%ndimlist(4:) = NMNHDIM_UNUSED
- end if
- zfieldp2d => pfield(jphext + 1, jphext + 1, :, :)
- if ( Present( koffset ) ) then
- ioffset2d(1) = koffset(3)
- ioffset2d(2) = koffset(4)
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp2d, ioffset2d, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp2d, ioffset2d, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp2d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp2d, iresp_nc4 )
- end if
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1:2) = NMNHDIM_ONE
- end if
- zfieldp => pfield(jphext + 1 : jphext + 1, jphext + 1 : jphext + 1, :, :)
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp, iresp_nc4 )
- end if
- endif
- ! ELSE IF (LPACK .AND. L2D .AND. YDIR=='XY') THEN
- ELSEIF (LPACK .AND. L2D .AND. SIZE(PFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 1
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = tzfield%ndimlist(3)
- tzfield%ndimlist(3) = tzfield%ndimlist(4)
- tzfield%ndimlist(4) = tzfield%ndimlist(5) !Necessary if time dimension
- tzfield%ndimlist(5:) = NMNHDIM_UNUSED
- end if
- zfieldp3d => pfield(:, jphext + 1, :, :)
- if ( Present( koffset ) ) then
- ioffset3d(1) = koffset(1)
- ioffset3d(2) = koffset(3)
- ioffset3d(3) = koffset(4)
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp3d, ioffset3d, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp3d, ioffset3d, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp3d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp3d, iresp_nc4 )
- end if
- else
- if ( tzfield%ndimlist(2) /= NMNHDIM_UNKNOWN ) tzfield%ndimlist(2) = NMNHDIM_ONE
- zfieldp => pfield(:, jphext + 1 : jphext + 1, :, :)
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp, iresp_nc4 )
- end if
- endif
- ELSE
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, pfield, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, pfield, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, pfield, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, pfield, iresp_nc4 )
- end if
- END IF
- ELSE
- CALL MPI_ALLREDUCE(SIZE(PFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_X4','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- CALL ALLOCBUFFER_ll(ZFIELDP,PFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(ZFIELDP(0,0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'GEN', 'IO_Field_write_byfield_X4', '2D not (yet) allowed for parallel execution' )
- CALL GATHER_XXFIELD('XX',PFIELD(:,JPHEXT+1,:,:),ZFIELDP(:,1,:,:),TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSE
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- END IF
- END IF ! multiprocess execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_X4', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- END SUBROUTINE IO_Field_write_byfield_X4
-
-
- SUBROUTINE IO_Field_write_byname_X5(TPFILE,HNAME,PFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- REAL,DIMENSION(:,:,:,:,:), INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_X5',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),PFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_X5
-
-
- SUBROUTINE IO_Field_write_byfield_X5(TPFILE,TPFIELD,PFIELD,KRESP)
- use modd_field, only: NMNHDIM_UNKNOWN, NMNHDIM_ONE, NMNHDIM_UNUSED
- USE MODD_IO, ONLY: GSMONOPROC, ISP, L1D, L2D, LPACK
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_TIMEZ, ONLY: TIMEZ
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- USE MODE_IO_TOOLS, ONLY: IO_Level2filenumber_get
- USE MODE_MNH_TIMING, ONLY: SECOND_MNH2
- USE MODD_VAR_ll, ONLY: MNH_STATUSES_IGNORE
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(:,:,:,:,:),TARGET,INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- real,dimension(:,:,:), pointer :: zfieldp3d
- real,dimension(:,:,:,:), pointer :: zfieldp4d
- REAL,DIMENSION(:,:,:,:,:),POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- INTEGER :: IHEXTOT
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- class(tfieldmetadata), allocatable :: tzfield
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- IHEXTOT = 2*JPHEXT+1
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_X5',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEREAL,5,'IO_Field_write_byfield_X5')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_X5',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- ! IF (LPACK .AND. L1D .AND. YDIR=='XY') THEN
- IF (LPACK .AND. L1D .AND. SIZE(PFIELD,1)==IHEXTOT .AND. SIZE(PFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 2
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1) = tzfield%ndimlist(3)
- tzfield%ndimlist(2) = tzfield%ndimlist(4)
- tzfield%ndimlist(3) = tzfield%ndimlist(5)
- tzfield%ndimlist(4) = tzfield%ndimlist(6) !Necessary if time dimension
- tzfield%ndimlist(5:) = NMNHDIM_UNUSED
- end if
- zfieldp3d => pfield(jphext + 1, jphext + 1, :, :, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp3d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp3d, iresp_nc4 )
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1:2) = NMNHDIM_ONE
- end if
- zfieldp => pfield(jphext + 1 : jphext + 1, jphext + 1 : jphext + 1, :, :, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp, iresp_nc4 )
- endif
- ! ELSE IF (LPACK .AND. L2D .AND. YDIR=='XY') THEN
- ELSEIF (LPACK .AND. L2D .AND. SIZE(PFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 1
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = tzfield%ndimlist(3)
- tzfield%ndimlist(3) = tzfield%ndimlist(4)
- tzfield%ndimlist(4) = tzfield%ndimlist(5)
- tzfield%ndimlist(5) = tzfield%ndimlist(6) !Necessary if time dimension
- tzfield%ndimlist(6:) = NMNHDIM_UNUSED
- end if
- zfieldp4d => pfield(:, jphext + 1, :, :, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp4d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp4d, iresp_nc4 )
- else
- if ( tzfield%ndimlist(2) /= NMNHDIM_UNKNOWN ) tzfield%ndimlist(2) = NMNHDIM_ONE
- zfieldp => pfield(:, jphext + 1 : jphext + 1, :, :, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, zfieldp, iresp_nc4 )
- endif
- ELSE
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,PFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,PFIELD,iresp_nc4)
- END IF
- ELSE
- CALL MPI_ALLREDUCE(SIZE(PFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_X5','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- CALL ALLOCBUFFER_ll(ZFIELDP,PFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(ZFIELDP(0,0,0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'GEN', 'IO_Field_write_byfield_X5', '2D not (yet) allowed for parallel execution' )
- CALL GATHER_XXFIELD('XX',PFIELD(:,JPHEXT+1,:,:,:),ZFIELDP(:,1,:,:,:),&
- & TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSE
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,ZFIELDP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,ZFIELDP,iresp_nc4)
- END IF
- END IF ! multiprocess execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_X5', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- END SUBROUTINE IO_Field_write_byfield_X5
-
-
- SUBROUTINE IO_Field_write_byname_X6(TPFILE,HNAME,PFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- REAL,DIMENSION(:,:,:,:,:,:), INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_X6',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),PFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_X6
-
- SUBROUTINE IO_Field_write_byfield_X6(TPFILE,TPFIELD,PFIELD,KRESP)
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_TIMEZ, ONLY: TIMEZ
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- USE MODE_IO_TOOLS, ONLY: IO_Level2filenumber_get
- USE MODE_MNH_TIMING, ONLY: SECOND_MNH2
- USE MODD_VAR_ll, ONLY: MNH_STATUSES_IGNORE
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- REAL,DIMENSION(:,:,:,:,:,:),TARGET,INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- REAL,DIMENSION(:,:,:,:,:,:),POINTER :: ZFIELDP
- LOGICAL :: GLFI, GNC4
- LOGICAL :: GALLOC
- INTEGER :: IHEXTOT
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- IHEXTOT = 2*JPHEXT+1
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_X6',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEREAL,6,'IO_Field_write_byfield_X6')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_X6',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,PFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,PFIELD,iresp_nc4)
- ELSE
- CALL MPI_ALLREDUCE(SIZE(PFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_X6','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- CALL ALLOCBUFFER_ll(ZFIELDP,PFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(ZFIELDP(0,0,0,0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,ZFIELDP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,ZFIELDP,iresp_nc4)
- END IF
- END IF ! multiprocess execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_X6', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- END SUBROUTINE IO_Field_write_byfield_X6
-
-
- SUBROUTINE IO_Field_write_byname_N0(TPFILE,HNAME,KFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- INTEGER, INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_N0',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),KFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_N0
-
-
- SUBROUTINE IO_Field_write_byfield_N0(TPFILE,TPFIELD,KFIELD,KRESP)
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !* 0. DECLARATIONS
- ! ------------
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- INTEGER, INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- INTEGER :: IK_FILE
- TYPE(TFILEDATA),POINTER :: TZFILE
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- TZFILE => NULL()
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_N0',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEINT,0,'IO_Field_write_byfield_N0')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_N0',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,KFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,KFIELD,iresp_nc4)
- ELSE
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,KFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,KFIELD,iresp_nc4)
- END IF
- END IF ! multiprocess execution
-#if MNH_SCALARS_IN_SPLITFILES
- IF (TPFILE%NSUBFILES_IOZ>0) THEN
- ! write the data in all Z files
- DO IK_FILE=1,TPFILE%NSUBFILES_IOZ
- TZFILE => TPFILE%TFILES_IOZ(IK_FILE)%TFILE
- IF ( ISP == TZFILE%NMASTER_RANK ) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TZFILE,TPFIELD,KFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TZFILE,TPFIELD,KFIELD,iresp_nc4)
- END IF
- END DO
- ENDIF
-#endif
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_N0', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
- END SUBROUTINE IO_Field_write_byfield_N0
-
-
- SUBROUTINE IO_Field_write_byname_N1(TPFILE,HNAME,KFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- INTEGER,DIMENSION(:), INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_N1',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),KFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_N1
-
-
- SUBROUTINE IO_Field_write_byfield_N1(TPFILE,TPFIELD,KFIELD,KRESP)
- !
- USE MODD_IO, ONLY: ISP,GSMONOPROC
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- INTEGER,DIMENSION(:),TARGET, INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- INTEGER,DIMENSION(:),POINTER :: IFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_N1',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEINT,1,'IO_Field_write_byfield_N1')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_N1',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,KFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,KFIELD,iresp_nc4)
- ELSE ! multiprocesses execution
- CALL MPI_ALLREDUCE(SIZE(KFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_N1','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- CALL ALLOCBUFFER_ll(IFIELDP,KFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(IFIELDP(0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,KFIELD,IFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,IFIELDP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,IFIELDP,iresp_nc4)
- END IF
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_N1', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(IFIELDP)
- END SUBROUTINE IO_Field_write_byfield_N1
-
-
- SUBROUTINE IO_Field_write_byname_N2(TPFILE,HNAME,KFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- INTEGER,DIMENSION(:,:), INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_N2',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),KFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_N2
-
-
- SUBROUTINE IO_Field_write_byfield_N2(TPFILE,TPFIELD,KFIELD,KRESP)
- use modd_field, only: NMNHDIM_UNKNOWN, NMNHDIM_ONE, NMNHDIM_UNUSED
- USE MODD_IO, ONLY: GSMONOPROC, ISP, L1D, L2D, LPACK
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_TIMEZ, ONLY: TIMEZ
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- USE MODE_MNH_TIMING, ONLY: SECOND_MNH2
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- INTEGER,DIMENSION(:,:),TARGET,INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- integer :: ifieldp0d
- integer,dimension(:), pointer :: ifieldp1d
- INTEGER,DIMENSION(:,:),POINTER :: IFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- !
- REAL(kind=MNHTIME), DIMENSION(2) :: ZT0, ZT1, ZT2
- REAL(kind=MNHTIME), DIMENSION(2) :: ZT11, ZT22
- INTEGER :: IHEXTOT
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- class(tfieldmetadata), allocatable :: tzfield
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- IHEXTOT = 2*JPHEXT+1
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_N2',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL SECOND_MNH2(ZT11)
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEINT,2,'IO_Field_write_byfield_N2')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_N2',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (LPACK .AND. L1D .AND. SIZE(KFIELD,1)==IHEXTOT .AND. SIZE(KFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 2
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1) = tzfield%ndimlist(3) !Necessary if time dimension
- tzfield%ndimlist(2:) = NMNHDIM_UNUSED
- end if
- ifieldp0d = kfield(jphext + 1, jphext + 1)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp0d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp0d, iresp_nc4 )
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1:2) = NMNHDIM_ONE
- end if
- ifieldp => kfield(jphext + 1 : jphext + 1, jphext + 1 : jphext + 1)
- if ( glfi) call IO_Field_write_lfi( tpfile, tzfield, ifieldp, iresp_lfi )
- if ( gnc4) call IO_Field_write_nc4( tpfile, tzfield, ifieldp, iresp_nc4 )
- endif
- ! ELSE IF (LPACK .AND. L2D .AND. YDIR=='XY') THEN
- ELSEIF (LPACK .AND. L2D .AND. SIZE(KFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 1
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = tzfield%ndimlist(3) !Necessary if time dimension
- tzfield%ndimlist(3:) = NMNHDIM_UNUSED
- end if
- ifieldp1d => kfield(:, jphext + 1)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp1d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp1d, iresp_nc4 )
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = NMNHDIM_ONE
- end if
- ifieldp => kfield(:, jphext + 1 : jphext + 1)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp, iresp_nc4 )
- endif
- ELSE
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,KFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,KFIELD,iresp_nc4)
- END IF
- ELSE ! multiprocesses execution
- CALL MPI_ALLREDUCE(SIZE(KFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_N2','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- CALL SECOND_MNH2(ZT0)
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- ! I/O process case
- CALL ALLOCBUFFER_ll(IFIELDP,KFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(IFIELDP(0,0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,KFIELD,IFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'GEN', 'IO_Field_write_byfield_N2', '2D not (yet) allowed for parallel execution' )
- CALL GATHER_XXFIELD('XX',KFIELD(:,JPHEXT+1),IFIELDP(:,1),TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSE
- CALL GATHER_XYFIELD(KFIELD,IFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- END IF
- CALL SECOND_MNH2(ZT1)
- TIMEZ%T_WRIT2D_GATH=TIMEZ%T_WRIT2D_GATH + ZT1 - ZT0
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,IFIELDP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,IFIELDP,iresp_nc4)
- END IF
- CALL SECOND_MNH2(ZT2)
- TIMEZ%T_WRIT2D_WRIT=TIMEZ%T_WRIT2D_WRIT + ZT2 - ZT1
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_N2', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(IFIELDP)
- CALL SECOND_MNH2(ZT22)
- TIMEZ%T_WRIT2D_ALL=TIMEZ%T_WRIT2D_ALL + ZT22 - ZT11
- !
- END SUBROUTINE IO_Field_write_byfield_N2
-
-
- SUBROUTINE IO_Field_write_byname_N3(TPFILE,HNAME,KFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- INTEGER,DIMENSION(:,:,:), INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_N3',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),KFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_N3
-
- SUBROUTINE IO_Field_write_byfield_N3(TPFILE,TPFIELD,KFIELD,KRESP)
- use modd_field, only: NMNHDIM_UNKNOWN, NMNHDIM_ONE, NMNHDIM_UNUSED
- USE MODD_IO, ONLY: GSMONOPROC, ISP, L1D, L2D, LPACK
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_TIMEZ, ONLY: TIMEZ
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- USE MODE_MNH_TIMING, ONLY: SECOND_MNH2
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- INTEGER,DIMENSION(:,:,:),TARGET,INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- integer, dimension(:), pointer :: ifieldp1d
- integer, dimension(:,:), pointer :: ifieldp2d
- INTEGER, DIMENSION(:,:,:), POINTER :: IFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- !
- REAL(kind=MNHTIME), DIMENSION(2) :: ZT11, ZT22
- INTEGER :: IHEXTOT
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- class(tfieldmetadata), allocatable :: tzfield
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- IHEXTOT = 2*JPHEXT+1
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_N3',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL SECOND_MNH2(ZT11)
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEINT,3,'IO_Field_write_byfield_N3')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_N3',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (LPACK .AND. L1D .AND. SIZE(KFIELD,1)==IHEXTOT .AND. SIZE(KFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 2
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1) = tzfield%ndimlist(3)
- tzfield%ndimlist(2) = tzfield%ndimlist(4) !Necessary if time dimension
- tzfield%ndimlist(3:) = NMNHDIM_UNUSED
- end if
- ifieldp1d => kfield(jphext + 1, jphext + 1, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp1d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp1d, iresp_nc4 )
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1:2) = NMNHDIM_ONE
- end if
- ifieldp => kfield(jphext + 1 : jphext + 1, jphext + 1 : jphext + 1, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp, iresp_nc4 )
- endif
- ! ELSE IF (LPACK .AND. L2D .AND. YDIR=='XY') THEN
- ELSEIF (LPACK .AND. L2D .AND. SIZE(KFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 1
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = tzfield%ndimlist(3)
- tzfield%ndimlist(3) = tzfield%ndimlist(4) !Necessary if time dimension
- tzfield%ndimlist(4:) = NMNHDIM_UNUSED
- end if
- ifieldp2d => kfield(:, jphext + 1, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp2d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp2d, iresp_nc4 )
- else
- if ( tzfield%ndimlist(2) /= NMNHDIM_UNKNOWN ) tzfield%ndimlist(2) = NMNHDIM_ONE
- ifieldp => kfield(:, jphext + 1 : jphext + 1, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp, iresp_nc4 )
- endif
- ELSE
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,KFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,KFIELD,iresp_nc4)
- END IF
- ELSE ! multiprocesses execution
- CALL MPI_ALLREDUCE(SIZE(KFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_N3','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- ! I/O process case
- CALL ALLOCBUFFER_ll(IFIELDP,KFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(IFIELDP(0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,KFIELD,IFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'GEN', 'IO_Field_write_byfield_N3', '2D not (yet) allowed for parallel execution' )
- CALL GATHER_XXFIELD('XX',KFIELD(:,JPHEXT+1,:),IFIELDP(:,1,:),TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSE
- CALL GATHER_XYFIELD(KFIELD,IFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,IFIELDP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,IFIELDP,iresp_nc4)
- END IF
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_N3', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(IFIELDP)
- CALL SECOND_MNH2(ZT22)
- TIMEZ%T_WRIT3D_ALL=TIMEZ%T_WRIT3D_ALL + ZT22 - ZT11
- !
- END SUBROUTINE IO_Field_write_byfield_N3
-
-
- SUBROUTINE IO_Field_write_byname_N4( TPFILE, HNAME, KFIELD, KRESP, koffset )
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- INTEGER,DIMENSION(:,:,:,:), INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(4), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_N4',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- if( iresp == 0 ) call IO_Field_write( tpfile, tfieldlist(id), kfield, iresp, koffset )
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_N4
-
-
- SUBROUTINE IO_Field_write_byfield_N4( TPFILE, TPFIELD, KFIELD, KRESP, koffset )
- use modd_field, only: NMNHDIM_UNKNOWN, NMNHDIM_ONE, NMNHDIM_UNUSED
- USE MODD_IO, ONLY: GSMONOPROC, ISP, L1D, L2D, LPACK
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_TIMEZ, ONLY: TIMEZ
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- USE MODE_IO_TOOLS, ONLY: IO_Level2filenumber_get
- USE MODE_MNH_TIMING, ONLY: SECOND_MNH2
- USE MODD_VAR_ll, ONLY: MNH_STATUSES_IGNORE
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- INTEGER,DIMENSION(:,:,:,:),TARGET, INTENT(IN) :: KFIELD ! array containing the data field
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(4), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- integer, dimension(2) :: ioffset2d
- integer, dimension(3) :: ioffset3d
- integer, dimension(:,:), pointer :: ifieldp2d
- integer, dimension(:,:,:), pointer :: ifieldp3d
- integer, dimension(:,:,:,:), pointer :: ifieldp
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- INTEGER :: IHEXTOT
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- class(tfieldmetadata), allocatable :: tzfield
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- IHEXTOT = 2*JPHEXT+1
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_N4',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEINT,4,'IO_Field_write_byfield_N4')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_N4',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- if ( Present( koffset ) .and. glfi ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_partial_write_byfield_N4', Trim( tpfile%cname ) // ': LFI format not supported' )
- glfi = .false.
- end if
-
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- ! IF (LPACK .AND. L1D .AND. YDIR=='XY') THEN
- IF (LPACK .AND. L1D .AND. SIZE(KFIELD,1)==IHEXTOT .AND. SIZE(KFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 2
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1) = tzfield%ndimlist(3)
- tzfield%ndimlist(2) = tzfield%ndimlist(4)
- tzfield%ndimlist(3) = tzfield%ndimlist(5) !Necessary if time dimension
- tzfield%ndimlist(4:) = NMNHDIM_UNUSED
- end if
- ifieldp2d => kfield(jphext + 1, jphext + 1, :, :)
- if ( Present( koffset ) ) then
- ioffset2d(1) = koffset(3)
- ioffset2d(2) = koffset(4)
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, ifieldp2d, ioffset2d, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, ifieldp2d, ioffset2d, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp2d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp2d, iresp_nc4 )
- end if
- else
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(1:2) = NMNHDIM_ONE
- end if
- ifieldp => kfield(jphext + 1 : jphext + 1, jphext + 1 : jphext + 1, :, :)
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, ifieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, ifieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp, iresp_nc4 )
- end if
- endif
- ! ELSE IF (LPACK .AND. L2D .AND. YDIR=='XY') THEN
- ELSEIF (LPACK .AND. L2D .AND. SIZE(KFIELD,2)==IHEXTOT) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 1
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = tzfield%ndimlist(3)
- tzfield%ndimlist(3) = tzfield%ndimlist(4)
- tzfield%ndimlist(4) = tzfield%ndimlist(5) !Necessary if time dimension
- tzfield%ndimlist(5:) = NMNHDIM_UNUSED
- end if
- ifieldp3d => kfield(:, jphext + 1, :, :)
- if ( Present( koffset ) ) then
- ioffset3d(1) = koffset(1)
- ioffset3d(2) = koffset(3)
- ioffset3d(3) = koffset(4)
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, ifieldp3d, ioffset3d, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, ifieldp3d, ioffset3d, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp3d, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp3d, iresp_nc4 )
- end if
- else
- if ( tzfield%ndimlist(2) /= NMNHDIM_UNKNOWN ) tzfield%ndimlist(2) = NMNHDIM_ONE
- ifieldp => kfield(:, jphext + 1 : jphext + 1, :, :)
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tzfield, ifieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tzfield, ifieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ifieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ifieldp, iresp_nc4 )
- end if
- endif
- ELSE
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, kfield, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, kfield, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, kfield, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, kfield, iresp_nc4 )
- end if
- END IF
- ELSE
- CALL MPI_ALLREDUCE(SIZE(KFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_N4','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- CALL ALLOCBUFFER_ll(IFIELDP,KFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(IFIELDP(0,0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,KFIELD,IFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSEIF (YDIR == 'XY') THEN
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'GEN', 'IO_Field_write_byfield_N4', '2D not (yet) allowed for parallel execution' )
- CALL GATHER_XXFIELD('XX',KFIELD(:,JPHEXT+1,:,:),IFIELDP(:,1,:,:),TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- ELSE
- CALL GATHER_XYFIELD(KFIELD,IFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, ifieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, ifieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, ifieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, ifieldp, iresp_nc4 )
- end if
- END IF
- END IF ! multiprocess execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_N4', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(IFIELDP)
- END SUBROUTINE IO_Field_write_byfield_N4
-
-
- SUBROUTINE IO_Field_write_byname_L0(TPFILE,HNAME,OFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- LOGICAL, INTENT(IN) :: OFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_L0',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),OFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_L0
-
- SUBROUTINE IO_Field_write_byfield_L0(TPFILE,TPFIELD,OFIELD,KRESP)
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- USE MODE_IO_MANAGE_STRUCT, ONLY: IO_File_find_byname
- !* 0. DECLARATIONS
- ! ------------
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- LOGICAL, INTENT(IN) :: OFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- INTEGER :: IK_FILE
- LOGICAL :: GLFI, GNC4
- TYPE(TFILEDATA),POINTER :: TZFILE
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- TZFILE => NULL()
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_L0',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPELOG,0,'IO_Field_write_byfield_L0')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_L0',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,OFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,OFIELD,iresp_nc4)
- ELSE
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,OFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,OFIELD,iresp_nc4)
- END IF
- END IF ! multiprocesses execution
-#if MNH_SCALARS_IN_SPLITFILES
- IF (TPFILE%NSUBFILES_IOZ>0) THEN
- ! write the data in all Z files
- DO IK_FILE=1,TPFILE%NSUBFILES_IOZ
- TZFILE => TPFILE%TFILES_IOZ(IK_FILE)%TFILE
- IF ( ISP == TZFILE%NMASTER_RANK ) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TZFILE,TPFIELD,OFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TZFILE,TPFIELD,OFIELD,iresp_nc4)
- END IF
- END DO
- ENDIF
-#endif
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_L0', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
- END SUBROUTINE IO_Field_write_byfield_L0
-
-
- SUBROUTINE IO_Field_write_byname_L1(TPFILE,HNAME,OFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- LOGICAL,DIMENSION(:), INTENT(IN) :: OFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_L1',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),OFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_L1
-
-
- SUBROUTINE IO_Field_write_byfield_L1(TPFILE,TPFIELD,OFIELD,KRESP)
- !
- USE MODD_IO, ONLY: ISP, GSMONOPROC
- !
- USE MODE_ALLOCBUFFER_ll
- USE MODE_GATHER_ll
- !
- IMPLICIT NONE
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- LOGICAL,DIMENSION(:),TARGET, INTENT(IN) :: OFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=2) :: YDIR ! field form
- INTEGER :: IERR
- INTEGER :: ISIZEMAX
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- LOGICAL,DIMENSION(:),POINTER :: GFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YDIR = TPFIELD%CDIR
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_L1',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPELOG,1,'IO_Field_write_byfield_L1')
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_L1',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,OFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,OFIELD,iresp_nc4)
- ELSE ! multiprocesses execution
- CALL MPI_ALLREDUCE(SIZE(OFIELD),ISIZEMAX,1,MNHINT_MPI,MPI_MAX,TPFILE%NMPICOMM,IERR)
- IF (ISIZEMAX==0) THEN
- CALL PRINT_MSG(NVERB_INFO,'IO','IO_Field_write_byfield_L1','ignoring variable with a zero size ('//TRIM(YRECFM)//')')
- IF (PRESENT(KRESP)) KRESP=0
- RETURN
- END IF
-
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- CALL ALLOCBUFFER_ll(GFIELDP,OFIELD,YDIR,GALLOC)
- ELSE
- ALLOCATE(GFIELDP(0))
- GALLOC = .TRUE.
- END IF
- !
- IF (YDIR == 'XX' .OR. YDIR =='YY') THEN
- CALL GATHER_XXFIELD(YDIR,OFIELD,GFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM)
- END IF
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,GFIELDP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,GFIELDP,iresp_nc4)
- END IF
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_L1', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(GFIELDP)
- END SUBROUTINE IO_Field_write_byfield_L1
-
-
- SUBROUTINE IO_Field_write_byname_C0(TPFILE,HNAME,HFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- CHARACTER(LEN=*), INTENT(IN) :: HFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_C0',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),HFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_C0
-
-
- SUBROUTINE IO_Field_write_byfield_C0(TPFILE,TPFIELD,HFIELD,KRESP)
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- !* 0. DECLARATIONS
- ! ------------
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- CHARACTER(LEN=*), INTENT(IN) :: HFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_C0',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPECHAR,0,'IO_Field_write_byfield_C0')
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (LEN(HFIELD)==0 .AND. GLFI) THEN
- CALL PRINT_MSG(NVERB_ERROR,'IO','IO_Field_write_byfield_C0',&
- 'zero-size string not allowed if LFI output for '//TRIM(TPFIELD%CMNHNAME))
- END IF
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_C0',IRESP)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,HFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,HFIELD,iresp_nc4)
- ELSE
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,HFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,HFIELD,iresp_nc4)
- END IF
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_C0', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
- END SUBROUTINE IO_Field_write_byfield_C0
-
-
- SUBROUTINE IO_Field_write_byname_C1(TPFILE,HNAME,HFIELD,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- CHARACTER(LEN=*),DIMENSION(:),INTENT(IN) :: HFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_C1',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),HFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_C1
-
-
- SUBROUTINE IO_Field_write_byfield_C1(TPFILE,TPFIELD,HFIELD,KRESP)
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- !* 0. DECLARATIONS
- ! ------------
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- CHARACTER(LEN=*),DIMENSION(:),INTENT(IN) :: HFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- INTEGER :: J,JJ
- INTEGER :: ILE, IP
- INTEGER,DIMENSION(:),ALLOCATABLE :: IFIELD
- INTEGER :: ILENG
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_C1',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPECHAR,1,'IO_Field_write_byfield_C1')
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF(GLFI) THEN
- ILE=LEN(HFIELD)
- IP=SIZE(HFIELD)
- ILENG=ILE*IP
- !
- IF (ILENG==0) THEN
- IP=1
- ILE=1
- ILENG=1
- ALLOCATE(IFIELD(1))
- IFIELD(1)=IACHAR(' ')
- ELSE
- ALLOCATE(IFIELD(ILENG))
- DO JJ=1,IP
- DO J=1,ILE
- IFIELD(ILE*(JJ-1)+J)=IACHAR(HFIELD(JJ)(J:J))
- END DO
- END DO
- END IF
- END IF
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_C1',IRESP)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,IFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,HFIELD,iresp_nc4)
- ELSE
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,IFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,HFIELD,iresp_nc4)
- END IF
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_C1', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (ALLOCATED(IFIELD)) DEALLOCATE(IFIELD)
- END SUBROUTINE IO_Field_write_byfield_C1
-
-
- SUBROUTINE IO_Field_write_byname_T0(TPFILE,HNAME,TFIELD,KRESP)
- USE MODD_TYPE_DATE, only: DATE_TIME
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- TYPE (DATE_TIME), INTENT(IN) :: TFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_T0',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),TFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_T0
-
-
- SUBROUTINE IO_Field_write_byfield_T0(TPFILE,TPFIELD,TFIELD,KRESP)
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- USE MODD_TYPE_DATE, only: DATE_TIME
- !
- !* 0. DECLARATIONS
- ! ------------
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- TYPE (DATE_TIME), INTENT(IN) :: TFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_T0',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEDATE,0,'IO_Field_write_byfield_T0')
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_T0',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,TFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,TFIELD,iresp_nc4)
- ELSE
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,TFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,TFIELD,iresp_nc4)
- END IF
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_T0', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
- END SUBROUTINE IO_Field_write_byfield_T0
-
-
- SUBROUTINE IO_Field_write_byname_T1(TPFILE,HNAME,TFIELD,KRESP)
- USE MODD_TYPE_DATE, only: DATE_TIME
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- TYPE (DATE_TIME),DIMENSION(:), INTENT(IN) :: TFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_T1',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write(TPFILE,TFIELDLIST(ID),TFIELD,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_T1
-
-
- SUBROUTINE IO_Field_write_byfield_T1(TPFILE,TPFIELD,TFIELD,KRESP)
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- USE MODD_TYPE_DATE, only: DATE_TIME
- !
- !* 0. DECLARATIONS
- ! ------------
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- TYPE (DATE_TIME),DIMENSION(:), INTENT(IN) :: TFIELD ! array containing the data field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_T1',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- CALL IO_Field_metadata_check(TPFIELD,TYPEDATE,1,'IO_Field_write_byfield_T1')
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_T1',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,TFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,TFIELD,iresp_nc4)
- ELSE
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,TFIELD,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,TFIELD,iresp_nc4)
- END IF
- END IF
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_T1', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
- END SUBROUTINE IO_Field_write_byfield_T1
-
-
- SUBROUTINE IO_Field_write_byname_lb(TPFILE,HNAME,KL3D,PLB,KRESP)
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CHARACTER(LEN=*), INTENT(IN) :: HNAME ! name of the field to write
- INTEGER, INTENT(IN) :: KL3D ! size of the LB array in FM
- REAL,DIMENSION(:,:,:), INTENT(IN) :: PLB ! array containing the LB field
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- INTEGER :: ID ! Index of the field
- INTEGER :: IRESP ! return_code
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byname_lb',TRIM(TPFILE%CNAME)//': writing '//TRIM(HNAME))
- !
- CALL FIND_FIELD_ID_FROM_MNHNAME(HNAME,ID,IRESP)
- !
- IF(IRESP==0) CALL IO_Field_write_lb(TPFILE,TFIELDLIST(ID),KL3D,PLB,IRESP)
- !
- IF (PRESENT(KRESP)) KRESP = IRESP
- !
- END SUBROUTINE IO_Field_write_byname_lb
-
-
- SUBROUTINE IO_Field_write_byfield_lb(TPFILE,TPFIELD,KL3D,PLB,KRESP)
- !
- use modd_field, only: NMNHDIM_UNKNOWN, NMNHDIM_ONE, NMNHDIM_UNUSED
- USE MODD_IO, ONLY: GSMONOPROC, ISNPROC, ISP, L2D, LPACK
- USE MODD_PARAMETERS_ll, ONLY: JPHEXT
- USE MODD_VAR_ll, ONLY: MNH_STATUSES_IGNORE
- !
- USE MODE_DISTRIB_lb
- USE MODE_TOOLS_ll, ONLY: GET_GLOBALDIMS_ll
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(INOUT) :: TPFIELD
- INTEGER, INTENT(IN) :: KL3D ! size of the LB array in FM
- REAL,DIMENSION(:,:,:),TARGET,INTENT(IN) :: PLB ! array containing the LB field
- INTEGER,OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- CHARACTER(LEN=28) :: YFILEM ! FM-file name
- CHARACTER(LEN=NMNHNAMELGTMAX) :: YRECFM ! name of the article to write
- CHARACTER(LEN=4) :: YLBTYPE ! 'LBX','LBXU','LBY' or 'LBYV'
- INTEGER :: IRIM ! size of the LB area
- INTEGER :: IERR
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- REAL,DIMENSION(:,:,:),ALLOCATABLE,TARGET :: Z3D
- real,dimension(:,:), pointer :: ztx2dp
- REAL,DIMENSION(:,:,:), POINTER :: TX3DP
- INTEGER :: IIMAX_ll,IJMAX_ll
- INTEGER :: JI
- INTEGER :: IIB,IIE,IJB,IJE
- INTEGER, DIMENSION(MPI_STATUS_SIZE) :: ISTATUS
- INTEGER,ALLOCATABLE,DIMENSION(:) :: IREQ_TAB
- INTEGER :: INB_REQ,IKU
- LOGICAL :: GLFI, GNC4
- TYPE TX_3DP
- REAL,DIMENSION(:,:,:), POINTER :: X
- END TYPE TX_3DP
- TYPE(TX_3DP),ALLOCATABLE,DIMENSION(:) :: T_TX3DP
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- class(tfieldmetadata), allocatable :: tzfield
- !
- YFILEM = TPFILE%CNAME
- YRECFM = TPFIELD%CMNHNAME
- YLBTYPE = TPFIELD%CLBTYPE
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_byfield_lb',TRIM(YFILEM)//': writing '//TRIM(YRECFM))
- !
- IF (YLBTYPE/='LBX' .AND. YLBTYPE/='LBXU' .AND. YLBTYPE/='LBY' .AND. YLBTYPE/='LBYV') THEN
- CALL PRINT_MSG(NVERB_ERROR,'IO','IO_Field_write_byfield_lb','unknown LBTYPE ('//YLBTYPE//')')
- RETURN
- END IF
- !
- IF (TPFIELD%CDIR/='') THEN
- CALL PRINT_MSG(NVERB_WARNING,'IO','IO_Field_write_byfield_lb','CDIR was set for '//TRIM(YRECFM))
- TPFIELD%CDIR=''
- END IF
- !
- IRIM = (KL3D-2*JPHEXT)/2
- IF (KL3D /= 2*(IRIM+JPHEXT)) THEN
- IRESP = -30
- GOTO 1000
- END IF
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_byfield_lb',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (LPACK .AND. L2D) THEN
- Allocate( tzfield, source = tpfield )
- if ( tpfile%ldimreduced ) then
- tzfield%ndims = tzfield%ndims - 1
- if ( tzfield%ndimlist(1) /= NMNHDIM_UNKNOWN ) then
- tzfield%ndimlist(2) = tzfield%ndimlist(3)
- tzfield%ndimlist(3) = tzfield%ndimlist(4) !Necessary if time dimension
- tzfield%ndimlist(4:) = NMNHDIM_UNUSED
- end if
- ztx2dp => plb(:, jphext + 1, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, ztx2dp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, ztx2dp, iresp_nc4 )
- else
- if ( tzfield%ndimlist(2) /= NMNHDIM_UNKNOWN ) tzfield%ndimlist(2) = NMNHDIM_ONE
- tx3dp => plb(:, jphext + 1 : jphext + 1, :)
- if ( glfi ) call IO_Field_write_lfi( tpfile, tzfield, tx3dp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tzfield, tx3dp, iresp_nc4 )
- endif
- ELSE
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,PLB,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,PLB,iresp_nc4)
- END IF
- ELSE
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- ! I/O proc case
- CALL GET_GLOBALDIMS_ll(IIMAX_ll,IJMAX_ll)
- IF (YLBTYPE == 'LBX' .OR. YLBTYPE == 'LBXU') THEN
- ALLOCATE(Z3D((IRIM+JPHEXT)*2,IJMAX_ll+2*JPHEXT,SIZE(PLB,3)))
- ELSE ! YLBTYPE == 'LBY' .OR. YLBTYPE == 'LBYV'
- ALLOCATE(Z3D(IIMAX_ll+2*JPHEXT,(IRIM+JPHEXT)*2,SIZE(PLB,3)))
- END IF
- DO JI = 1,ISNPROC
- CALL GET_DISTRIB_lb(YLBTYPE,JI,'FM','WRITE',IRIM,IIB,IIE,IJB,IJE)
- IF (IIB /= 0) THEN
- TX3DP=>Z3D(IIB:IIE,IJB:IJE,:)
- IF (ISP /= JI) THEN
- CALL MPI_RECV(TX3DP,SIZE(TX3DP),MNHREAL_MPI,JI-1,99,TPFILE%NMPICOMM,ISTATUS,IERR)
- ELSE
- CALL GET_DISTRIB_lb(YLBTYPE,JI,'LOC','WRITE',IRIM,IIB,IIE,IJB,IJE)
- TX3DP = PLB(IIB:IIE,IJB:IJE,:)
- END IF
- END IF
- END DO
- IF (LPACK .AND. L2D) THEN
- call Print_msg( NVERB_FATAL, 'GEN', 'IO_Field_write_byfield_lb', '2D not (yet) allowed for parallel execution' )
- TX3DP=>Z3D(:,JPHEXT+1:JPHEXT+1,:)
- ELSE
- TX3DP=>Z3D
- END IF
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,TX3DP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,TX3DP,iresp_nc4)
- ELSE
- INB_REQ=0
- ALLOCATE(IREQ_TAB(1))
- ALLOCATE(T_TX3DP(1))
- IKU = SIZE(PLB,3)
- ! Other processes
- CALL GET_DISTRIB_lb(YLBTYPE,ISP,'LOC','WRITE',IRIM,IIB,IIE,IJB,IJE)
- IF (IIB /= 0) THEN
- TX3DP=>PLB(IIB:IIE,IJB:IJE,:)
- INB_REQ = INB_REQ + 1
- ALLOCATE(T_TX3DP(INB_REQ)%X(IIB:IIE,IJB:IJE,IKU))
- T_TX3DP(INB_REQ)%X=PLB(IIB:IIE,IJB:IJE,:)
- CALL MPI_ISEND(T_TX3DP(INB_REQ)%X,SIZE(TX3DP),MNHREAL_MPI,TPFILE%NMASTER_RANK-1,99, &
- TPFILE%NMPICOMM,IREQ_TAB(INB_REQ),IERR)
- !CALL MPI_BSEND(TX3DP,SIZE(TX3DP),MNHREAL_MPI,TPFILE%NMASTER_RANK-1,99,TPFILE%NMPICOMM,IERR)
- END IF
- IF (INB_REQ .GT.0 ) THEN
- CALL MPI_WAITALL(INB_REQ,IREQ_TAB,MNH_STATUSES_IGNORE,IERR)
- DEALLOCATE(T_TX3DP(1)%X)
- END IF
- DEALLOCATE(T_TX3DP,IREQ_TAB)
- END IF
- END IF
- END IF
- !
-1000 CONTINUE
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_byfield_lb', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (ALLOCATED(Z3D)) DEALLOCATE(Z3D)
- IF (PRESENT(KRESP)) KRESP = IRESP
- END SUBROUTINE IO_Field_write_byfield_lb
-
-
- SUBROUTINE IO_Field_write_box_byfield_X2( TPFILE, TPFIELD, HBUDGET, PFIELD, KXOBOX, KXEBOX, KYOBOX, KYEBOX, KRESP, koffset )
- !
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- USE MODE_GATHER_ll
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- CHARACTER(LEN=*), INTENT(IN) :: HBUDGET ! 'BUDGET' (budget) or 'OTHER' (MesoNH field)
- REAL, DIMENSION(:,:), TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, INTENT(IN) :: KXOBOX !
- INTEGER, INTENT(IN) :: KXEBOX ! Global coordinates of the box
- INTEGER, INTENT(IN) :: KYOBOX !
- INTEGER, INTENT(IN) :: KYEBOX !
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(2), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- REAL, DIMENSION(:,:), POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_box_byfield_X2',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_box_byfield_X2',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- if ( Present( koffset ) .and. glfi ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_write_box_byfield_X3', Trim( tpfile%cname ) &
- // ': LFI format not supported' )
- glfi = .false.
- end if
-
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (HBUDGET /= 'BUDGET') THEN
- ! take the sub-section of PFIELD defined by the box
- ZFIELDP=>PFIELD(KXOBOX:KXEBOX,KYOBOX:KYEBOX)
- ELSE
- ! take the field as a budget
- ZFIELDP=>PFIELD
- END IF
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- ELSE ! multiprocesses execution
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- ! Allocate the box
- ALLOCATE(ZFIELDP(KXEBOX-KXOBOX+1,KYEBOX-KYOBOX+1))
- GALLOC = .TRUE.
- ELSE
- ALLOCATE(ZFIELDP(0,0))
- GALLOC = .TRUE.
- END IF
- !
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM,&
- & KXOBOX,KXEBOX,KYOBOX,KYEBOX,HBUDGET)
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- END IF
- END IF ! multiprocesses execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_box_byfield_X2', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- END SUBROUTINE IO_Field_write_box_byfield_X2
-
-
- SUBROUTINE IO_Field_write_box_byfield_X3( TPFILE, TPFIELD, HBUDGET, PFIELD, KXOBOX, KXEBOX, KYOBOX, KYEBOX, KRESP, koffset )
- !
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- USE MODE_GATHER_ll
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- CHARACTER(LEN=*), INTENT(IN) :: HBUDGET ! 'BUDGET' (budget) or 'OTHER' (MesoNH field)
- REAL, DIMENSION(:,:,:), TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, INTENT(IN) :: KXOBOX !
- INTEGER, INTENT(IN) :: KXEBOX ! Global coordinates of the box
- INTEGER, INTENT(IN) :: KYOBOX !
- INTEGER, INTENT(IN) :: KYEBOX !
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(3), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- REAL, DIMENSION(:,:,:), POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_box_byfield_X3',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_box_byfield_X3',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- if ( Present( koffset ) .and. glfi ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_write_box_byfield_X3', Trim( tpfile%cname ) &
- // ': LFI format not supported' )
- glfi = .false.
- end if
-
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (HBUDGET /= 'BUDGET') THEN
- ! take the sub-section of PFIELD defined by the box
- ZFIELDP=>PFIELD(KXOBOX:KXEBOX,KYOBOX:KYEBOX,:)
- ELSE
- ! take the field as a budget
- ZFIELDP=>PFIELD
- END IF
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- ELSE ! multiprocesses execution
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- ! Allocate the box
- ALLOCATE(ZFIELDP(KXEBOX-KXOBOX+1,KYEBOX-KYOBOX+1,SIZE(PFIELD,3)))
- GALLOC = .TRUE.
- ELSE
- ALLOCATE(ZFIELDP(0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM,&
- & KXOBOX,KXEBOX,KYOBOX,KYEBOX,HBUDGET)
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- END IF
- END IF ! multiprocesses execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_box_byfield_X3', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- END SUBROUTINE IO_Field_write_box_byfield_X3
-
-
- SUBROUTINE IO_Field_write_box_byfield_X4( TPFILE, TPFIELD, HBUDGET, PFIELD, KXOBOX, KXEBOX, KYOBOX, KYEBOX, KRESP, koffset )
- !
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- USE MODE_GATHER_ll
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- CHARACTER(LEN=*), INTENT(IN) :: HBUDGET ! 'BUDGET' (budget) or 'OTHER' (MesoNH field)
- REAL, DIMENSION(:,:,:,:), TARGET, INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, INTENT(IN) :: KXOBOX !
- INTEGER, INTENT(IN) :: KXEBOX ! Global coordinates of the box
- INTEGER, INTENT(IN) :: KYOBOX !
- INTEGER, INTENT(IN) :: KYEBOX !
- INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
- integer, dimension(4), optional, intent(in) :: koffset
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- REAL, DIMENSION(:,:,:,:), POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_box_byfield_X4',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_box_byfield_X4',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- if ( Present( koffset ) .and. glfi ) then
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Field_partial_write_box_byfield_X4', Trim( tpfile%cname ) &
- // ': LFI format not supported' )
- glfi = .false.
- end if
-
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (HBUDGET /= 'BUDGET') THEN
- ! take the sub-section of PFIELD defined by the box
- ZFIELDP=>PFIELD(KXOBOX:KXEBOX,KYOBOX:KYEBOX,:,:)
- ELSE
- ! take the field as a budget
- ZFIELDP=>PFIELD
- END IF
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- ELSE ! multiprocesses execution
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- ! Allocate the box
- ALLOCATE(ZFIELDP(KXEBOX-KXOBOX+1,KYEBOX-KYOBOX+1,SIZE(PFIELD,3),SIZE(PFIELD,4)))
- GALLOC = .TRUE.
- ELSE
- ALLOCATE(ZFIELDP(0,0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM,&
- & KXOBOX,KXEBOX,KYOBOX,KYEBOX,HBUDGET)
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- if ( Present( koffset ) ) then
- !if ( glfi ) call IO_Field_partial_write_lfi( tpfile, tpfield, zfieldp, koffset, iresp_lfi )
- if ( gnc4 ) call IO_Field_partial_write_nc4( tpfile, tpfield, zfieldp, koffset, iresp_nc4 )
- else
- if ( glfi ) call IO_Field_write_lfi( tpfile, tpfield, zfieldp, iresp_lfi )
- if ( gnc4 ) call IO_Field_write_nc4( tpfile, tpfield, zfieldp, iresp_nc4 )
- end if
- END IF
- END IF ! multiprocesses execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_box_byfield_X4', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- END SUBROUTINE IO_Field_write_box_byfield_X4
-
-
- SUBROUTINE IO_Field_write_box_byfield_X5(TPFILE,TPFIELD,HBUDGET,PFIELD,KXOBOX,KXEBOX,KYOBOX,KYEBOX,KRESP)
- !
- USE MODD_IO, ONLY: GSMONOPROC, ISP
- !
- USE MODE_GATHER_ll
- !
- !
- !* 0.1 Declarations of arguments
- !
- TYPE(TFILEDATA), INTENT(IN) :: TPFILE
- CLASS(TFIELDMETADATA), INTENT(IN) :: TPFIELD
- CHARACTER(LEN=*), INTENT(IN) :: HBUDGET ! 'BUDGET' (budget) or 'OTHER' (MesoNH field)
- REAL,DIMENSION(:,:,:,:,:),TARGET,INTENT(IN) :: PFIELD ! array containing the data field
- INTEGER, INTENT(IN) :: KXOBOX !
- INTEGER, INTENT(IN) :: KXEBOX ! Global coordinates of the box
- INTEGER, INTENT(IN) :: KYOBOX !
- INTEGER, INTENT(IN) :: KYEBOX !
- INTEGER,OPTIONAL, INTENT(OUT):: KRESP ! return-code
- !
- !* 0.2 Declarations of local variables
- !
- integer :: iresp, iresp_lfi, iresp_nc4, iresp_glob
- REAL,DIMENSION(:,:,:,:,:),POINTER :: ZFIELDP
- LOGICAL :: GALLOC
- LOGICAL :: GLFI, GNC4
- CHARACTER(LEN=:),ALLOCATABLE :: YMSG
- CHARACTER(LEN=6) :: YRESP
- !
- CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_write_box_byfield_X5',TRIM(TPFILE%CNAME)//': writing '//TRIM(TPFIELD%CMNHNAME))
- !
- iresp = 0
- iresp_lfi = 0
- iresp_nc4 = 0
- GALLOC = .FALSE.
- !
- CALL IO_File_write_check(TPFILE,'IO_Field_write_box_byfield_X5',IRESP)
- !
- CALL IO_Format_write_select(TPFILE,GLFI,GNC4)
- !
- IF (IRESP==0) THEN
- IF (GSMONOPROC) THEN ! sequential execution
- IF (HBUDGET /= 'BUDGET') THEN
- ! take the sub-section of PFIELD defined by the box
- ZFIELDP=>PFIELD(KXOBOX:KXEBOX,KYOBOX:KYEBOX,:,:,:)
- ELSE
- ! take the field as a budget
- ZFIELDP=>PFIELD
- END IF
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,ZFIELDP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,ZFIELDP,iresp_nc4)
- ELSE ! multiprocesses execution
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- ! Allocate the box
- ALLOCATE(ZFIELDP(KXEBOX-KXOBOX+1,KYEBOX-KYOBOX+1,SIZE(PFIELD,3),&
- & SIZE(PFIELD,4),SIZE(PFIELD,5)))
- GALLOC = .TRUE.
- ELSE
- ALLOCATE(ZFIELDP(0,0,0,0,0))
- GALLOC = .TRUE.
- END IF
- !
- CALL GATHER_XYFIELD(PFIELD,ZFIELDP,TPFILE%NMASTER_RANK,TPFILE%NMPICOMM,&
- & KXOBOX,KXEBOX,KYOBOX,KYEBOX,HBUDGET)
- !
- IF (ISP == TPFILE%NMASTER_RANK) THEN
- IF (GLFI) CALL IO_Field_write_lfi(TPFILE,TPFIELD,ZFIELDP,iresp_lfi)
- IF (GNC4) CALL IO_Field_write_nc4(TPFILE,TPFIELD,ZFIELDP,iresp_nc4)
- END IF
- END IF ! multiprocesses execution
- END IF
-
- call IO_Field_write_error_check( tpfile, tpfield, 'IO_Field_write_box_byfield_X5', iresp, iresp_lfi, iresp_nc4, iresp_glob )
- if ( Present( kresp ) ) kresp = iresp_glob
-
- IF (GALLOC) DEALLOCATE(ZFIELDP)
- END SUBROUTINE IO_Field_write_box_byfield_X5
-
-
-SUBROUTINE IO_Fieldlist_write(TPOUTPUT)
-!
-USE MODE_MODELN_HANDLER, ONLY: GET_CURRENT_MODEL_INDEX
-!
-IMPLICIT NONE
-!
-TYPE(TOUTBAK), INTENT(IN) :: TPOUTPUT !Output structure
-!
-INTEGER :: IDX
-INTEGER :: IMI
-INTEGER :: JI
-!
-IMI = GET_CURRENT_MODEL_INDEX()
-!
-DO JI = 1,SIZE(TPOUTPUT%NFIELDLIST)
- IDX = TPOUTPUT%NFIELDLIST(JI)
- NDIMS: SELECT CASE (TFIELDLIST(IDX)%NDIMS)
- !
- !0D output
- !
- CASE (0)
- NTYPE0D: SELECT CASE (TFIELDLIST(IDX)%NTYPE)
- !
- !0D real
- !
- CASE (TYPEREAL)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_X0D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X0D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_X0D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X0D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_X0D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 0D real fields' )
- END IF
- !
- !0D integer
- !
- CASE (TYPEINT)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_N0D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_N0D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_N0D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_N0D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_N0D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 0D integer fields' )
- END IF
- !
- !0D logical
- !
- CASE (TYPELOG)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_L0D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_L0D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_L0D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_L0D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_L0D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 0D logical fields' )
- END IF
- !
- !0D string
- !
- CASE (TYPECHAR)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_C0D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_C0D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_C0D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_C0D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_C0D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 0D character fields' )
- END IF
- !
- !0D date/time
- !
- CASE (TYPEDATE)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_T0D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_T0D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_T0D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_T0D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_T0D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 0D date/time fields' )
- END IF
- !
- !0D other types
- !
- CASE DEFAULT
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': type not yet supported for 0D output' )
- END SELECT NTYPE0D
- !
- !1D output
- !
- CASE (1)
- NTYPE1D: SELECT CASE (TFIELDLIST(IDX)%NTYPE)
- !
- !1D real
- !
- CASE (TYPEREAL)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_X1D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X1D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_X1D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X1D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_X1D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 1D real fields' )
- END IF
- !
- !1D integer
- !
- CASE (TYPEINT)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_N1D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_N1D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_N1D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_N1D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_N1D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 1D integer fields' )
- END IF
- !
- !1D logical
- !
- CASE (TYPELOG)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_L1D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_L1D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_L1D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_L1D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_L1D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 1D logical fields' )
- END IF
- !
- !1D string
- !
- CASE (TYPECHAR)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_C1D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_C1D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_C1D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_C1D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_C1D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 1D character fields' )
- END IF
- !
- !1D date/time
- !
- CASE (TYPEDATE)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_T1D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_T1D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_T1D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_T1D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_T1D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 1D date/time fields' )
- END IF
- !
- !1D other types
- !
- CASE DEFAULT
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': type not yet supported for 1D output' )
- END SELECT NTYPE1D
- !
- !2D output
- !
- CASE (2)
- NTYPE2D: SELECT CASE (TFIELDLIST(IDX)%NTYPE)
- !
- !2D real
- !
- CASE (TYPEREAL)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_X2D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X2D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_X2D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X2D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_X2D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 2D real fields' )
- END IF
- !
- !2D integer
- !
- CASE (TYPEINT)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_N2D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_N2D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_N2D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_N2D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_N2D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not allowed for 2D integer fields' )
- END IF
- !
- !2D other types
- !
- CASE DEFAULT
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': type not yet supported for 2D output' )
- END SELECT NTYPE2D
- !
- !3D output
- !
- CASE (3)
- NTYPE3D: SELECT CASE (TFIELDLIST(IDX)%NTYPE)
- !
- !3D real
- !
- CASE (TYPEREAL)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_X3D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X3D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_X3D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X3D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_X3D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not (yet) allowed for 3D real fields' )
- !PW: TODO?: add missing field in TFIELDLIST?
- !CALL IO_Field_write_lb(TPOUTPUT%TFILE,TFIELDLIST(IDX),***,TFIELDLIST(IDX)%TFIELD_X3D(IMI)%DATA)
- END IF
- !
- !3D integer
- !
- CASE (TYPEINT)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_N3D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_N3D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_N3D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_N3D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_N3D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not (yet) allowed for 3D integer fields' )
- !PW: TODO?: add missing field in TFIELDLIST?
- !CALL IO_Field_write_lb(TPOUTPUT%TFILE,TFIELDLIST(IDX),***,TFIELDLIST(IDX)%TFIELD_N3D(IMI)%DATA)
- END IF
- !
- !3D other types
- !
- CASE DEFAULT
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': type not yet supported for 3D output' )
- END SELECT NTYPE3D
- !
- !4D output
- !
- CASE (4)
- NTYPE4D: SELECT CASE (TFIELDLIST(IDX)%NTYPE)
- !
- !4D real
- !
- CASE (TYPEREAL)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_X4D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X4D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_X4D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X4D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_X4D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not (yet) allowed for 4D real fields' )
- !PW: TODO?: add missing field in TFIELDLIST?
- !CALL IO_Field_write_lb(TPOUTPUT%TFILE,TFIELDLIST(IDX),***,TFIELDLIST(IDX)%TFIELD_X4D(IMI)%DATA)
- END IF
- !
- !4D other types
- !
- CASE DEFAULT
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': type not yet supported for 4D output' )
- END SELECT NTYPE4D
- !
- !5D output
- !
- CASE (5)
- NTYPE5D: SELECT CASE (TFIELDLIST(IDX)%NTYPE)
- !
- !5D real
- !
- CASE (TYPEREAL)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_X5D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X5D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_X5D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X5D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_X5D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not (yet) allowed for 5D real fields' )
- !PW: TODO?: add missing field in TFIELDLIST?
- !CALL IO_Field_write_lb(TPOUTPUT%TFILE,TFIELDLIST(IDX),***,TFIELDLIST(IDX)%TFIELD_X5D(IMI)%DATA)
- END IF
- !
- !5D other types
- !
- CASE DEFAULT
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': type not yet supported for 5D output' )
- END SELECT NTYPE5D
- !
- !6D output
- !
- CASE (6)
- NTYPE6D: SELECT CASE (TFIELDLIST(IDX)%NTYPE)
- !
- !6D real
- !
- CASE (TYPEREAL)
- IF ( .NOT.ALLOCATED(TFIELDLIST(IDX)%TFIELD_X6D) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X6D is NOT allocated ' )
- END IF
- IF ( .NOT.ASSOCIATED(TFIELDLIST(IDX)%TFIELD_X6D(IMI)%DATA) ) THEN
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': TFIELD_X6D%DATA is NOT associated' )
- END IF
- IF ( TFIELDLIST(IDX)%CLBTYPE == 'NONE' ) THEN
- CALL IO_Field_write(TPOUTPUT%TFILE,TFIELDLIST(IDX),TFIELDLIST(IDX)%TFIELD_X6D(IMI)%DATA)
- ELSE
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': CLBTYPE/=NONE not (yet) allowed for 6D real fields' )
- !PW: TODO?: add missing field in TFIELDLIST?
- !CALL IO_Field_write_lb(TPOUTPUT%TFILE,TFIELDLIST(IDX),***,TFIELDLIST(IDX)%TFIELD_X6D(IMI)%DATA)
- END IF
- !
- !6D other types
- !
- CASE DEFAULT
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': type not yet supported for 6D output' )
- END SELECT NTYPE6D
- !
- !Other number of dimensions
- !
- CASE DEFAULT
- call Print_msg( NVERB_ERROR, 'IO', 'IO_Fieldlist_write', trim(tfieldlist(idx)%cmnhname)// &
- ': number of dimensions not yet supported' )
- END SELECT NDIMS
-END DO
-!
-END SUBROUTINE IO_Fieldlist_write
-
-
-SUBROUTINE IO_Field_user_write(TPOUTPUT)
-!
-#if 0
-USE MODD_DYN_n, ONLY: XTSTEP
-USE MODD_FIELD_n, ONLY: XUT, XVT, XRT, XTHT, XSVT
-USE MODD_PARAMETERS, ONLY: JPVEXT
-USE MODD_PRECIP_n, ONLY: XINPRR
-#endif
-!
-IMPLICIT NONE
-!
-TYPE(TOUTBAK), INTENT(IN) :: TPOUTPUT !Output structure
-!
-TYPE(TFIELDMETADATA) :: TZFIELD
-!
-#if 0
-INTEGER :: IKB
-!
-IKB=JPVEXT+1
-!
-TZFIELD%CMNHNAME = 'UTLOW'
-TZFIELD%CSTDNAME = 'x_wind'
-TZFIELD%CLONGNAME = ''
-TZFIELD%CUNITS = 'm s-1'
-TZFIELD%CDIR = 'XY'
-TZFIELD%CCOMMENT = 'X_Y_Z_U component of wind at lowest physical level'
-TZFIELD%NGRID = 2
-TZFIELD%NTYPE = TYPEREAL
-TZFIELD%NDIMS = 2
-TZFIELD%LTIMEDEP = .TRUE.
-CALL IO_Field_write(TPOUTPUT%TFILE,TZFIELD,XUT(:,:,IKB))
-!
-TZFIELD%CMNHNAME = 'VTLOW'
-TZFIELD%CSTDNAME = 'y_wind'
-TZFIELD%CLONGNAME = ''
-TZFIELD%CUNITS = 'm s-1'
-TZFIELD%CDIR = 'XY'
-TZFIELD%CCOMMENT = 'X_Y_Z_V component of wind at lowest physical level'
-TZFIELD%NGRID = 3
-TZFIELD%NTYPE = TYPEREAL
-TZFIELD%NDIMS = 2
-TZFIELD%LTIMEDEP = .TRUE.
-CALL IO_Field_write(TPOUTPUT%TFILE,TZFIELD,XVT(:,:,IKB))
-!
-TZFIELD%CMNHNAME = 'THTLOW'
-TZFIELD%CSTDNAME = 'air_potential_temperature'
-TZFIELD%CLONGNAME = ''
-TZFIELD%CUNITS = 'K'
-TZFIELD%CDIR = 'XY'
-TZFIELD%CCOMMENT = 'X_Y_Z_potential temperature at lowest physical level'
-TZFIELD%NGRID = 1
-TZFIELD%NTYPE = TYPEREAL
-TZFIELD%NDIMS = 2
-TZFIELD%LTIMEDEP = .TRUE.
-CALL IO_Field_write(TPOUTPUT%TFILE,TZFIELD,XTHT(:,:,IKB))
-!
-TZFIELD%CMNHNAME = 'RVTLOW'
-!TZFIELD%CSTDNAME = 'humidity_mixing_ratio' !ratio of the mass of water vapor to the mass of dry air
-TZFIELD%CSTDNAME = 'specific_humidity' !mass fraction of water vapor in (moist) air
-TZFIELD%CLONGNAME = ''
-TZFIELD%CUNITS = 'kg kg-1'
-TZFIELD%CDIR = 'XY'
-TZFIELD%CCOMMENT = 'X_Y_Z_Vapor mixing Ratio at lowest physical level'
-TZFIELD%NGRID = 1
-TZFIELD%NTYPE = TYPEREAL
-TZFIELD%NDIMS = 2
-TZFIELD%LTIMEDEP = .TRUE.
-CALL IO_Field_write(TPOUTPUT%TFILE,TZFIELD,XRT(:,:,IKB,1))
-!
-TZFIELD%CMNHNAME = 'ACPRRSTEP'
-TZFIELD%CSTDNAME = 'rainfall_amount'
-TZFIELD%CLONGNAME = ''
-TZFIELD%CUNITS = 'kg m-2'
-TZFIELD%CDIR = ''
-TZFIELD%CCOMMENT = 'X_Y_ACcumulated Precipitation Rain Rate during timestep'
-TZFIELD%NGRID = 1
-TZFIELD%NTYPE = TYPEREAL
-TZFIELD%NDIMS = 2
-TZFIELD%LTIMEDEP = .TRUE.
-!XACPRR is multiplied by 1000. to convert from m to kg m-2 (water density is assumed to be 1000 kg m-3)
-CALL IO_Field_write(TPOUTPUT%TFILE,TZFIELD,XINPRR*XTSTEP*1.0E3)
-!
-TZFIELD%CMNHNAME = 'SVT001'
-TZFIELD%CSTDNAME = 'concentration in scalar variable'
-TZFIELD%CLONGNAME = ''
-TZFIELD%CUNITS = 'kg kg-1'
-TZFIELD%CDIR = 'XY'
-TZFIELD%CCOMMENT = 'X_Y_Z_concentration in scalar variable'
-TZFIELD%NGRID = 1
-TZFIELD%NTYPE = TYPEREAL
-TZFIELD%NDIMS = 3
-TZFIELD%LTIMEDEP = .TRUE.
-CALL IO_WRITE_FIELD(TPOUTPUT%TFILE,TZFIELD,XSVT(:,:,:,1))
-!
-#endif
-!
-END SUBROUTINE IO_Field_user_write
-
-END MODULE MODE_IO_FIELD_WRITE
-
diff --git a/src/mesonh/aux/mode_io_field_write_phy.f90 b/src/mesonh/aux/mode_io_field_write_phy.f90
new file mode 100644
index 0000000000000000000000000000000000000000..785677a4a967a515e81bfa2ff8360deed6c0f115
--- /dev/null
+++ b/src/mesonh/aux/mode_io_field_write_phy.f90
@@ -0,0 +1,93 @@
+!MNH_LIC Copyright 2022-2023 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
+!MNH_LIC for details. version 1.
+!-----------------------------------------------------------------
+! Modifications:
+! Q.Rodier 02/2023 Creation call to mode_io_field_write inside PHYEX
+!-----------------------------------------------------------------
+!
+MODULE MODE_IO_FIELD_WRITE_PHY
+ USE MODD_IO, ONLY: TFILEDATA
+ USE MODD_FIELD, ONLY: TFIELDMETADATA
+ IMPLICIT NONE
+ INTERFACE IO_Field_write_phy
+ MODULE PROCEDURE IO_Field_write_phy_byfield_X2, IO_Field_write_phy_byfield_X1
+ END INTERFACE
+CONTAINS
+ SUBROUTINE IO_Field_write_phy_byfield_X2(D, TPFILE, TPFIELD, PFIELD, KRESP, koffset )
+ USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
+ !
+ IMPLICIT NONE
+ !
+ !* 0.1 Declarations of arguments
+ !
+ TYPE(DIMPHYEX_t), INTENT(IN) :: D
+ TYPE(TFILEDATA), INTENT(IN) :: TPFILE
+ TYPE(TFIELDMETADATA), INTENT(IN) :: TPFIELD
+ REAL,DIMENSION(D%NIJT,D%NKT),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
+ INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
+ integer, dimension(3), optional, intent(in) :: koffset
+ !
+ CALL IO_Field_write_phy_unpack2D(D,TPFILE, TPFIELD, PFIELD, KRESP, koffset )
+ !
+ END SUBROUTINE IO_Field_write_phy_byfield_X2
+!
+ SUBROUTINE IO_Field_write_phy_unpack2D(D, TPFILE, TPFIELD, PFIELD, KRESP, koffset )
+ USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
+ USE MODE_IO_FIELD_WRITE, ONLY: IO_Field_write
+ !
+ IMPLICIT NONE
+ !
+ !* 0.1 Declarations of arguments
+ !
+ TYPE(DIMPHYEX_t), INTENT(IN) :: D
+ TYPE(TFILEDATA), INTENT(IN) :: TPFILE
+ TYPE(TFIELDMETADATA), INTENT(IN) :: TPFIELD
+ REAL,DIMENSION(D%NIT,D%NJT,D%NKT),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
+ INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
+ integer, dimension(3), optional, intent(in) :: koffset
+ !
+ CALL IO_Field_write(TPFILE, TPFIELD, PFIELD, KRESP, koffset )
+ !
+ END SUBROUTINE IO_Field_write_phy_unpack2D
+!
+ SUBROUTINE IO_Field_write_phy_byfield_X1(D, TPFILE, TPFIELD, PFIELD, KRESP, koffset )
+ USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
+ !
+ IMPLICIT NONE
+ !
+ !* 0.1 Declarations of arguments
+ !
+ TYPE(DIMPHYEX_t), INTENT(IN) :: D
+ TYPE(TFILEDATA), INTENT(IN) :: TPFILE
+ TYPE(TFIELDMETADATA), INTENT(IN) :: TPFIELD
+ REAL,DIMENSION(D%NIJT),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
+ INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
+ integer, dimension(3), optional, intent(in) :: koffset
+ !
+ CALL IO_Field_write_phy_unpack1D(D,TPFILE, TPFIELD, PFIELD, KRESP, koffset )
+ !
+ END SUBROUTINE IO_Field_write_phy_byfield_X1
+!
+ SUBROUTINE IO_Field_write_phy_unpack1D(D, TPFILE, TPFIELD, PFIELD, KRESP, koffset )
+ USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
+ USE MODE_IO_FIELD_WRITE, ONLY: IO_Field_write
+ !
+ IMPLICIT NONE
+ !
+ !* 0.1 Declarations of arguments
+ !
+ TYPE(DIMPHYEX_t), INTENT(IN) :: D
+ TYPE(TFILEDATA), INTENT(IN) :: TPFILE
+ TYPE(TFIELDMETADATA), INTENT(IN) :: TPFIELD
+ REAL,DIMENSION(D%NIT,D%NJT),TARGET, INTENT(IN) :: PFIELD ! array containing the data field
+ INTEGER, OPTIONAL, INTENT(OUT) :: KRESP ! return-code
+ integer, dimension(3), optional, intent(in) :: koffset
+ !
+ CALL IO_Field_write(TPFILE, TPFIELD, PFIELD, KRESP, koffset )
+ !
+ END SUBROUTINE IO_Field_write_phy_unpack1D
+!
+!
+END MODULE MODE_IO_FIELD_WRITE_PHY
diff --git a/src/mesonh/aux/shuman_phy.f90 b/src/mesonh/aux/shuman_phy.f90
index 98ab1fd072d37acd923bc6c7a417c848c9934999..e3e493afe37c2a8414312116367dfad2c3c0bbc1 100644
--- a/src/mesonh/aux/shuman_phy.f90
+++ b/src/mesonh/aux/shuman_phy.f90
@@ -1,4 +1,4 @@
-MODULE SHUMAN_PHY
+MODULE MODE_SHUMAN_PHY
IMPLICIT NONE
CONTAINS
! ###############################
@@ -1605,4 +1605,4 @@ END DO
!
END SUBROUTINE DYF_PHY
!
-END MODULE SHUMAN_PHY
+END MODULE MODE_SHUMAN_PHY
diff --git a/src/mesonh/filesToSuppress.txt b/src/mesonh/filesToSuppress.txt
index 6ed551e32c24e786a0e3607c71a6b12c9d3fe49d..7a2f0c47d76783189ca323990cece0a931eae1b6 100644
--- a/src/mesonh/filesToSuppress.txt
+++ b/src/mesonh/filesToSuppress.txt
@@ -15,6 +15,7 @@ aux/modi_shuman.f90
aux/modi_second_mnh.f90
aux/second_mnh.f90
aux/modi_gradient_m.f90
+aux/mode_io_field_write.f90
#To be un-commented when the PHYEX version will be merged into the Meso-NH code
#aux/mode_thermo.f90