- Timestamp:
- Aug 2, 2024, 9:58:25 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/writedynav_loc.f90
r5158 r5159 16 16 USE lmdz_comgeom 17 17 18 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 19 USE lmdz_paramet 18 20 IMPLICIT NONE 19 21 20 ! 22 21 23 ! Ecriture du fichier histoire au format IOIPSL 22 ! 24 23 25 ! Appels succesifs des routines: histwrite 24 ! 26 25 27 ! Entree: 26 28 ! histid: ID du fichier histoire … … 34 36 ! ps :pression au sol 35 37 ! phis : geopotentiel au sol 36 ! 37 ! 38 39 38 40 ! Sortie: 39 41 ! fileid: ID du fichier netcdf cree 40 ! 42 41 43 ! L. Fairhead, LMD, 03/99 42 ! 44 43 45 ! ===================================================================== 44 ! 46 45 47 ! Declarations 46 INCLUDE "dimensions.h" 47 INCLUDE "paramet.h" 48 49 ! 48 49 50 51 50 52 ! Arguments 51 53 ! … … 62 64 ! This routine needs IOIPSL 63 65 ! Variables locales 64 ! 66 65 67 INTEGER, SAVE, ALLOCATABLE :: ndex2d(:), ndexu(:), ndexv(:) 66 68 INTEGER :: iq, ii, ll … … 73 75 !$OMP THREADPRIVATE(first) 74 76 75 ! 77 76 78 ! Initialisations 77 ! 79 78 80 IF (adjust) return 79 81 … … 101 103 CALL covnat_loc(llm, ucov, vcov, unat, vnat) 102 104 103 ! 105 104 106 ! Appels a histwrite pour l'ecriture des variables a sauvegarder 105 ! 107 106 108 ! Vents U 107 109 ! … … 117 119 !$OMP END MASTER 118 120 119 ! 121 120 122 ! Vents V 121 ! 123 122 124 ije = ij_end 123 125 IF (pole_sud) jjn = jj_nb - 1 … … 130 132 131 133 132 ! 134 133 135 ! Temperature potentielle moyennee 134 ! 136 135 137 ijb = ij_begin 136 138 ije = ij_end … … 141 143 !$OMP END MASTER 142 144 143 ! 145 144 146 ! Temperature moyennee 145 147 ! … … 159 161 160 162 161 ! 163 162 164 ! Geopotentiel 163 ! 165 164 166 !$OMP MASTER 165 167 CALL histwrite(histaveid, 'phi', itau_w, phi(ijb:ije, :), & … … 168 170 169 171 170 ! 172 171 173 ! Traceurs 172 ! 174 173 175 !!$OMP MASTER 174 176 ! DO iq=1,nqtot … … 179 181 180 182 181 ! 183 182 184 ! Masse 183 ! 185 184 186 !$OMP MASTER 185 187 CALL histwrite(histaveid, 'masse', itau_w, masse(ijb:ije, :), & … … 188 190 189 191 190 ! 192 191 193 ! Pression au sol 192 ! 194 193 195 !$OMP MASTER 194 196 … … 197 199 !$OMP END MASTER 198 200 199 ! 201 200 202 ! Geopotentiel au sol 201 ! 203 202 204 !$OMP MASTER 203 205 ! CALL histwrite(histaveid, 'phis', itau_w, phis(ijb:ije), … … 205 207 !$OMP END MASTER 206 208 207 ! 209 208 210 ! Fin 209 ! 211 210 212 !$OMP MASTER 211 213 IF (ok_sync) THEN
Note: See TracChangeset
for help on using the changeset viewer.