File indexing completed on 2024-12-17 18:38:01 UTC
view on githubraw file Latest commit b55e95f1 on 2018-09-19 15:37:37 UTC
b55e95f1ff Oliv*0001 #include "RADTRANS_OPTIONS.h"
0002
0003
0004
0005
0006
0007 SUBROUTINE RADTRANS_DIAGNOSTICS_INIT( myThid )
0008
0009
0010
0011
0012
0013 IMPLICIT NONE
0014 #include "SIZE.h"
0015 #include "EEPARAMS.h"
0016 #include "PARAMS.h"
0017 #include "RADTRANS_SIZE.h"
0018
0019
0020
0021 INTEGER myThid
0022
0023
0024 #ifdef ALLOW_RADTRANS
0025 #ifdef ALLOW_DIAGNOSTICS
0026
0027
0028 INTEGER ILNBLNK
0029 EXTERNAL ILNBLNK
0030 CHARACTER*16 DIAGS_MK_UNITS
0031 EXTERNAL DIAGS_MK_UNITS
0032
0033
0034
0035
0036 INTEGER j
0037 INTEGER diagNum
0038 CHARACTER*8 diagName
0039 CHARACTER*16 diagCode
0040 CHARACTER*16 diagUnits
0041 CHARACTER*(80) diagTitle
0042 INTEGER ilnb
0043 CHARACTER*7 trcUnits
0044 CHARACTER*30 locName
0045 CHARACTER*4 diagSufx
0046
0047
0048
0049 WRITE(diagName,'(A)')'rmud'
0050 WRITE(diagTitle,'(A)')'inverse cosine of solar zenith angle'
0051 diagUnits = '1 '
0052 diagCode = 'SM P L1 '
0053 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0054 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0055
0056 DO j = 1, nlam
0057 WRITE(diagName,'(A,I3.3)')'Rirr', j
0058 WRITE(diagTitle,'(A, I3)')
0059 & 'irradiance reflectance for waveband ', j
0060 diagUnits = '1 '
0061 diagCode = 'SM P L1 '
0062 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0063 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0064
0065 WRITE(diagName,'(A,I3.3)')'Ed',j
0066 WRITE(diagTitle,'(A,I3)')'direct irradiance for waveband ',j
0067 diagUnits = 'W/m^2 '
0068 diagCode = 'SMRP LR '
0069 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0070 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0071
0072 WRITE(diagName,'(A,I3.3)')'Es',j
0073 WRITE(diagTitle,'(A,I3)')
0074 & 'diffuse downwelling irradiance for waveband ',j
0075 diagUnits = 'W/m^2 '
0076 diagCode = 'SMRP LR '
0077 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0078 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0079
0080 WRITE(diagName,'(A,I3.3)')'Eu',j
0081 WRITE(diagTitle,'(A,I3)')
0082 & 'diffuse upwelling irradiance for waveband ',j
0083 diagUnits = 'W/m^2 '
0084 diagCode = 'SMRP LR '
0085 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0086 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0087
0088 WRITE(diagName,'(A,I3.3)')'E',j
0089 WRITE(diagTitle,'(A,I3)')'vector irradiance waveband ',j
0090 diagUnits = 'W/m^2 '
0091 diagCode = 'SMRP MR '
0092 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0093 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0094
0095 WRITE(diagName,'(A,I3.3)')'E0F',j
0096 WRITE(diagTitle,'(A,I3)')
0097 & 'scalar irradiance at W point waveband ',j
0098 diagUnits = 'W/m^2 '
0099 diagCode = 'SMRP LR '
0100 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0101 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0102
0103 #ifdef RADTRANS_DIAG_SOLUTION
0104 WRITE(diagName,'(A,I3.3)')'Estop',j
0105 WRITE(diagTitle,'(A,I3)')
0106 & 'diffuse downwelling irradiance for waveband ',j
0107 diagUnits = 'W/m^2 '
0108 diagCode = 'SMRP LR '
0109 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0110 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0111
0112 WRITE(diagName,'(A,I3.3)')'Eubot',j
0113 WRITE(diagTitle,'(A,I3)')
0114 & 'diffuse upwelling irradiance for waveband ',j
0115 diagUnits = 'W/m^2 '
0116 diagCode = 'SMRP LR '
0117 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0118 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0119
0120 WRITE(diagName,'(A,I3.3)')'amp1_',j
0121 WRITE(diagTitle,'(A,I3)')
0122 & 'amplitude of d-w increasing mode for waveband ',j
0123 diagUnits = 'W/m^2 '
0124 diagCode = 'SMR MR '
0125 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0126 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0127
0128 WRITE(diagName,'(A,I3.3)')'amp2_',j
0129 WRITE(diagTitle,'(A,I3)')
0130 & 'amplitude of d-w decreasing mode for waveband ',j
0131 diagUnits = 'W/m^2 '
0132 diagCode = 'SMR MR '
0133 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0134 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0135
0136 WRITE(diagName,'(A,I3.3)')'x_',j
0137 WRITE(diagTitle,'(A,I3)')
0138 & 'amplitude of direct in diffuse for waveband ',j
0139 diagUnits = 'W/m^2 '
0140 diagCode = 'SMR MR '
0141 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0142 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0143
0144 WRITE(diagName,'(A,I3.3)')'y_',j
0145 WRITE(diagTitle,'(A,I3)')
0146 & 'amplitude of direct in upward for waveband ',j
0147 diagUnits = 'W/m^2 '
0148 diagCode = 'SMR MR '
0149 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0150 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0151
0152 WRITE(diagName,'(A,I3.3)')'r1_',j
0153 WRITE(diagTitle,'(A,I3)')
0154 & 'R coefficient of d-w increasing mode for waveband ',j
0155 diagUnits = '1 '
0156 diagCode = 'SMRP MR '
0157 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0158 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0159
0160 WRITE(diagName,'(A,I3.3)')'r2_',j
0161 WRITE(diagTitle,'(A,I3)')
0162 & 'R coefficient of d-w decreasing mode for waveband ',j
0163 diagUnits = '1 '
0164 diagCode = 'SMRP MR '
0165 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0166 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0167
0168 WRITE(diagName,'(A,I3.3)')'att1_',j
0169 WRITE(diagTitle,'(A,I3)')
0170 & 'attenuation of d-w increasing mode for waveband ',j
0171 diagUnits = '1/m '
0172 diagCode = 'SMRP MR '
0173 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0174 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0175
0176 WRITE(diagName,'(A,I3.3)')'att2_',j
0177 WRITE(diagTitle,'(A,I3)')
0178 & 'attenuation of d-w decreasing mode for waveband ',j
0179 diagUnits = '1/m '
0180 diagCode = 'SMRP MR '
0181 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0182 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0183 #endif
0184 ENDDO
0185
0186 #ifdef ALLOW_EXF
0187 WRITE(diagName,'(A)')'RTice'
0188 WRITE(diagTitle,'(A)')'ice area fraction in radtrans package'
0189 diagUnits = 'm^2/m^2 '
0190 diagCode = 'SM P MR '
0191 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0192 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0193 #endif
0194
0195 #endif /* ALLOW_DIAGNOSTICS */
0196 #endif /* ALLOW_RADTRANS */
0197
0198 RETURN
0199 END
0200
