Changeset 3531 for LMDZ6/trunk/libf/phylmd/readaerosolstrato.F90
- Timestamp:
- Jun 6, 2019, 5:08:45 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/readaerosolstrato.F90
r3436 r3531 15 15 USE aero_mod 16 16 USE dimphy 17 USE print_control_mod, ONLY: prt_level,lunout 17 18 #ifdef CPP_XIOS 18 19 USE xios … … 61 62 data alpha_strat_wave/3.36780953,3.34667683,3.20444202,3.0293026,2.82108808/ 62 63 64 CHARACTER (len = 20) :: modname = 'readaerosolstrato' 65 CHARACTER (len = 80) :: abort_message 66 63 67 !-------------------------------------------------------- 64 68 … … 72 76 73 77 IF (nbands.NE.2) THEN 74 print *,'nbands doit etre egal a 2 dans readaerosolstrat'75 STOP78 abort_message='nbands doit etre egal a 2 dans readaerosolstrat' 79 CALL abort_physic(modname,abort_message,1) 76 80 ENDIF 77 81 … … 82 86 n_lev = size(lev) 83 87 IF (n_lev.NE.klev) THEN 84 print *,'Le nombre de niveaux n est pas egal a klev'85 STOP88 abort_message='Le nombre de niveaux n est pas egal a klev' 89 CALL abort_physic(modname,abort_message,1) 86 90 ENDIF 87 91 … … 89 93 CALL nf95_gw_var(ncid_in, varid, latitude) 90 94 n_lat = size(latitude) 91 print *,'LAT aerosol strato=', n_lat, latitude95 WRITE(lunout,*) 'LAT aerosol strato=', n_lat, latitude 92 96 IF (grid_type/=unstructured) THEN 93 97 IF (n_lat.NE.nbp_lat) THEN 94 print *,'Le nombre de lat n est pas egal a nbp_lat'95 STOP98 abort_message='Le nombre de lat n est pas egal a nbp_lat' 99 CALL abort_physic(modname,abort_message,1) 96 100 ENDIF 97 101 ENDIF … … 101 105 n_lon = size(longitude) 102 106 IF (grid_type/=unstructured) THEN 103 print *,'LON aerosol strato=', n_lon, longitude107 WRITE(lunout,*) 'LON aerosol strato=', n_lon, longitude 104 108 IF (n_lon.NE.nbp_lon) THEN 105 print *,'Le nombre de lon n est pas egal a nbp_lon'106 STOP109 abort_message='Le nombre de lon n est pas egal a nbp_lon' 110 CALL abort_physic(modname,abort_message,1) 107 111 ENDIF 108 112 ENDIF … … 111 115 CALL nf95_gw_var(ncid_in, varid, time) 112 116 n_month = size(time) 113 print *,'TIME aerosol strato=', n_month, time117 WRITE(lunout,*) 'TIME aerosol strato=', n_month, time 114 118 IF (n_month.NE.12) THEN 115 print *,'Le nombre de month n est pas egal a 12'116 STOP119 abort_message='Le nombre de month n est pas egal a 12' 120 CALL abort_physic(modname,abort_message,1) 117 121 ENDIF 118 122 … … 124 128 CALL nf95_inq_varid(ncid_in, "TAUSTRAT", varid) 125 129 ncerr = nf90_get_var(ncid_in, varid, tauaerstrat) 126 print *,'code erreur readaerosolstrato=', ncerr, varid130 WRITE(lunout,*) 'code erreur readaerosolstrato=', ncerr, varid 127 131 128 132 CALL nf95_close(ncid_in) … … 130 134 !---select the correct month 131 135 IF (mth_cur.LT.1.OR.mth_cur.GT.12) THEN 132 print *,'probleme avec le mois dans readaerosolstrat =', mth_cur136 WRITE(lunout,*) 'probleme avec le mois dans readaerosolstrat =', mth_cur 133 137 ENDIF 134 138 tauaerstrat_mois(:,:,:) = tauaerstrat(:,:,:,mth_cur)
Note: See TracChangeset
for help on using the changeset viewer.